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PREFACE 


\ 

Fliis  First  Teclinical  Report  describes  the  work  performed  on  Contract 
DAI  1C15-73-C-0127  on  Theoretical  Studies  of  Iligh-Power  Infrared  Window 
Materials  during  the  period  from  December  7,  1972  through  June  30,  1973. 

'Die  work  on  the  present  contract  is  a  continuation  of  that  of  the  previous 
Contract  DA1IC15-72-C-0129. 

The  following  investigators  contributed  to  this  report: 

Mr.  11.  C.  Chow,  research  associate 

Dr.  C.  j.  Duthler,  principal  research  scientist 

Dr.  A.  M.  Karo,  consultant,  Lawrence  Livermore  Laboratory, 

Livermore,  California 

Dr.  A.  A.  Maradudin,  consultant,  University  of  California,  Irvine,  California 

Dr.  D.  L.  Mills,  consultant,  University  of  California,  Irvine,  California 

Mr.  A.  Mcreira,  research  associate 

Dr.  L.  J.  Sham,  consultant.  University  of  California,  San  Diego,  California 

Dr.  M.  Sparks,  principal  investigator 

The  material  in  this  report  constitutes  the  final  results  on  the  subjects  cov¬ 
ered.  The  preliminary  discussions  of  the  theory  of  infrared  absorption  and  material 
failure  in  crystals  containing  inclusions  and  of  the  theory  of  multiphonon  infrared 
absorption  presented  in  the  reports  of  the  preceding  contract  DAHC15-72-C-0129 
are  superseded  by  the  results  of  die  present  report.  In  particular,  the  section  on 
the  effect  of  inclusion  of  optical  absorption  in  the  previous  Final  Report  contained 
an  error,  which  fortunately  is  not  of  practical  consequence.  The  appropriate  value 
of  reflection  is  the  average  value  (R>  rather  than  the  perpendicular  incidence 
value  R^  . 
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SUMMARY 


Theory  of  Infrared  Absorption  and  Material  Failure  in  Crystals  Containing 

Inclusions.  Two  effects  of  inclusions  in  or  on  the  surface  of  infrared-transmitting 

materials  are  to  increase  the  average  value  of  the  optical  absorption  coefficient  0 

and  to  cause  localized  heating  that  could  lead  to  material  failure  at  high -power 

-7  ~8 

levels.  Volume  fractions  as  low  as  10"  to  10  of  such  inclusions  can  give  rise 
to  a  value  of  the  optical  absorption  coefficient  /3  of  10  4  cm  \  a  typical  value  of 
current  interest.  For  various  types  of  inclusions,  the  frequen  cy  dependence  of  fi 
ranges  from  increasing  as  Oi  »  to  independent  of  co,  to  exponentially  decreasing 
with  a).  The  temperature  dependence  ranges  from  independent  of  T ,  to  increasing 
as  Tp  in  the  high-temperature  limit,  where  p  -  2-4  typically.  Simple  expres¬ 
sions  for  the  absorption  cross  section  are  derived  for  various  cases  of  practical 
interest.  The  cross  sections  are  used  to  derive  expressions  for  0  for  the  four 
cases  of  large  inclusions  of  strong  and  weak  absorbers  and  of  small  inclusions  of 
dielectric  and  metallic  particles.  The  material  failure  resulting  from  local  heating 
of  inclusions  is  a  far  greater  problem  in  high -intensity  short -pulse  systems  than  in 
low-intensity  long-pulse  or  cw  systems  having  the  same  average  intensity.  Micro¬ 
second  pulses  with  energy  densities  as  low  as  a  few  joules  per  square  centimeter 
can  cause  material  failure. 


Theory  of  Multiphon-  l  Absorption  in  Insulating  Crystals.  The  nearly  expo¬ 
nential  frequency  dependence  of  the  infrared  absorption  coefficient  0  recently 
observed  in  fifteen  crystals  up  to  several  Umes  the  reststrahl  frequency  is  ex¬ 
plained  in  terms  of  multiphonon  absorption  processes.  The  central-limit  theorem 
is  used  to  reduce  the  multiphonon  contribution  to  a  simple  closed  form.  The  theore¬ 
tical  estimates  for  the  magnitude  of  the  absorption  coefficient,  with  no  adjustable 
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parameters,  are  also  in  good  agreen.ent  with  experiment.  The  temperature 
dependence  of  0  at  a  fixed  frequency  is  shown  to  be  considerably  weaker  than 
0  „  Tn_1,  where  n  is  the  number  of  created  phonons.  Higher-order  processes 
m  the  perturbation  expansion  are  shown  to  be  negligible  for  small  n ,  to  be  com¬ 
parable  to  that  of  the  lowest-order,  single-vertex  terms  for  n  2:  S  ,  and  to  domi¬ 
nate  for  large  n  in  a  typical  case.  Difference  processes,  in  which  some  thermally 
excited  phonons  are  annihilated,  are  shown  to  be  negligible  with  respect  to  the 
summation  processes  in  the  nearly  exponential  region.  An  explanation  involving 
finite  phonon  lifetimes  is  proposed  to  explain  the  fact  chat  the  alkali  halides  show 
less  structure  in  the  0-  CO  curves  than  do  the  semiconductor  crystals. 

Temperature  Dependence  of  Multiphonon  Infrared  Absorption.  Measurements 
of  Harrington  and  Hass  and  of  Barker  indicate  that  the  temperature  dependence  of 
the  infrared  absorption  coefficient  0  in  the  n -phonon  region  is  considerably  weaker 
than  0  ~  Tn_1,  which  had  been  predicted  for  the  high-temperature  limit  of  multi¬ 
phonon  absorption.  This  discrepancy  is  resolved  by  taking  into  account  the  tempera¬ 
ture  dependence  of  the  phonon  frequencies  and  the  lattice  constant.  The  agreement 
between  the  experimental  and  theoretical  results  with  no  adjustable  parameters  is 
good.  A  new  evaluation  of  the  multiphonon  sums  yields  0  ~  exp  (-to  T  )  directly, 
rather  than  as  a  sum  on  n . 

Theory  of  Infrared  Absorption  by  Crystals  in  the  High  Frequency  Wing  of 
Their  Fundamental  Lattice  Absorption.  We  have  calculated  the  frequency  depend¬ 
ence  of  infrared  absorption  in  the  classical  limit  for  an  exactly  soluble  model  of  a 
lattice  of  noninteracting  diatomic  molecules,  each  bound  internally  by  a  potential 

for  which  the  classical  equation  of  motion  can  be  solved  in  closed  form.  Four  poten 

2 

tials  have  been  used:  a  Morse  potential,  a  potential  of  the  form  V  (x )  =  ( a/x  )  +  bx 
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an  infinite  square  well  potential,  and  a  triangular  well  potential.  The  analytic 
results  we  obtain  show  that  the  absorption  coefficient  for  large  frequencies  associ¬ 
ated  with  potentials  which  admit  an  harmonic  approximation  decreases  nearly  ex¬ 
ponentially  over  the  frequency  region  covered  by  recent  experiments,  with  signifi¬ 
cant  deviations  from  exponential  behavior  at  higher  frequencies.  For  the  square 

-2 

and  triangular  well  potentials,  the  absorption  decreases  like  oo  for  frequencies 
large  compared  to  a  characteristic  frequency. 

Temperature  Dependence  of  the  Absorption  Coefficient  of  Alkali  Halides  in 
the  Multiphonon  Regime.  The  theory  of  infrared  absorption  by  an  array  of  inde¬ 
pendent,  anharmonic  oscillators  is  discussed.  When  the  oscillator  potential  is  the 
Morse  potential,  the  theory  provides  an  excellent  description  of  the  temperature 
dependence  of  the  absorption  coefficient  at  10.  6jU  in  NaCl  and  NaF  reported  by 
Harrington  and  Hass. 

Temperature  and  Frequency  Dependence  of  Infrared  Absorption  as  a 
Diagnostic  Tool.  Recent  developments  render  untenable  a  proposed  method  of 
distinguishing  between  intrinsic  and  extrinsic  infrared  absorption  on  the  basis 
of  the  proposed  temperature  dependence.  However,  when  the  proper  temperature 
dependence  of  multiphonon  absorption  is  accounted  for  and  the  possibility  of  other 
intrinsic  processes  is  taken  into  account,  the  temperature  and  frequency  depend¬ 
ence  of  the  absorption  of  both  the  best  available  and  intentionally  imperfected 
crystals  should  be  useful  in  studying  extrinsic  processes. 

Short-Pulse  Operation  of  Infrared  Windows  without  Thermal  Defocusing. 

The  possibility  of  transmitting  short  infrared  pulses  through  materials  with  little 

thermally  induced  optical  distortion  is  shown  to  exist.  For  sufficiently  short  pulses, 
-8  -9 

of  the  order  of  10  -  10  sec,  the  absorbed  energy  does  not  have  time  to  thermalize. 


Thus,  the  thermally  induced  optical  distortion  is  greatly  reduced. 


Sec.  A 


A.  INTRODUCTION 

The  motivation  for  this  program  on  theoretical  studies  of  high-power  infrared 
window  materials,  which  is  a  continuation  of  a  previous  contract  DAHC15-72-C- 
0129,  was  the  availability  of  high-power  infrared  lasers  for  current  Department  of 
Defense  programs  and  the  realization  that  lack  of  transparent  materials  for  windows 
may  limit  the  usefulness  of  many  laser  systems.  Values  of  the  optical  absorption 
coefficient  0  of  candidate  window  materials  were  needed  in  order  to  evaluate  the 
potential  performance  of  the  materials.  There  had  been  no  previous  calculations 
of  the  numerical  values  of  j3  in  the  highly  transparent  regions  for  materials  of  in¬ 
terest  (such  as  KBr  and  ZnSe  at  10.6f.tm),  and  the  currently  available  correspond¬ 
ing  experimental  values  were  of  questionable  efficacy  since  they  were  believed  to  be 
extrinsic  (i.  e. ,  caused  by  imperfections  that  can  be  removed  in  principle. ). 

The  paucity  of  experimental  and  theoretical  information  on  the  values  of  /3 
was  one  of  the  most  pressing  problems  in  the  present  Department  of  Defense  high- 
power-window  programs.  It  was  especially  important  to  know  if  the  values  of  )3 
were  intrinsic  or  extrinsic  and  to  have  reliable  estimates  of  the  intrinsic  value  of  jS 
before  undertaking  imperfection -identification  and  sample-purification  programs 
since  th'ire  were  many  candidate  materials  and  these  expensive  programs  should  be 
undertaken  only  if  there  were  a  good  chance  of  reducing  |3  to  the  required  value. 

During  the  early  stages  of  the  previous  contract  it  became  increasingly  ap¬ 
parent  that  in  order  to  obtain  this  information  on  j3,  theoretical  and  experimental 
values  of  0  were  needed  not  only  at  10. ^m,  but  also  over  a  large  range  of  values 
of  frequencies  and  temperature.  In  the  intervening  eighteen  nonths  there  has  been 
considerable  progress  in  our  theoretical  program  and  in  experimerLdl  and  theoretical 
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Sec.  A 


programs  at  other  laboratories.  We  have  explained  quantitatively  the  nearly 
exponential  frequency  dependence  of  the  optical  absorption  frequency  /3  observed 
by  Rupprecht  and  by  Deutsch  and  the  substantial  deviations  from  the  expected 
temperature  dependence  observed  by  Harrington  and  Hass.  The  calculations  are 
based  on  a  reasonable  model  of  the  lattice  with  the  Born-Mayer  interaction  poten¬ 
tial.  They  include  the  dispersion  relations  of  the  phonons,  and  the  approximations 
made  were  shown  to  be  reasonable.  The  theory  of  intrinsic  multiphcnon  absorption 
now  appears  to  be  complete,  and  the  emphasis  of  the  program  has  shifted  to  ex¬ 
trinsic  and  nonlinear  and  other  high-power  absorption  mechanisms.  Although  no 
attempt  will  be  made  to  review  the  progress  of  other  laboratories,  it  should  be 
mentioned  that  the  Raytheon  measurements  of  0(u>)  and  our  theoretical  prediction 
that  there  should  be  no  drastic  deviations  from  the  extrapolations  of  the  measured 
/3 (co)  have  settled  the  question  of  whether  the  values  of  f}  measured  at  10.  6|im  are 
intrinsic  or  extrinsic  for  most  materials  of  interest  and  have  afforded  estimates  of 
the  intrinsic  values.  With  the  exception  of  KC1,  the  estimated  intrinsic  values  of 
^10  6  ^or  canc^date  10.  6fim  window  materials  are  well  below  the  lowest  measured 
values. 


A  study  (Sec.  C  )  of  tile  effects  of  macroscopic  inclusions  in  crystals,  including 
the  increase  in  absorption  and  damage  thresholds,  has  been  completed,  A  study  of 
a  proposed  quasi-selection  rule  for  absorption  and  the  first  phase  of  a  study  of  the 
effects  on  absorption  of  parametric  processes  are  nearing  completion,  and  a  number 
of  other  problems,  listed  below,  are  under  investigation.  The  following  publications 
and  reports  have  been  prepared  under  this  and  the  previous  contract: 
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M.  Sparks  and  T.  Azzarelli,  "Theoretical  Studies  of  High-Power  Infrared  Window 
Materials,"  Xonics  Quarterly  Technical  Progress  Report  No.  1,  Contract  DAHC15- 
72-C-0129,  March  1972. 

M.  Sparks,  "Recent  Developments  in  High-Power  Infrared  Window  Research,"  Invited 
Talk,  4th  ASTM  Damage  in  Laser  Materials  Symposium,  Boulder,  Colorado,  June 
14-15,  1972. 

M.  Sparks  and  T.  Azzarelli,  "Theoretical  Studies  of  High-Power  Infrared  Window 
Materials,"  Xonics  Quarterly  Technical  Progress  Report  No. 2,  Contract  DAHC15- 
72-C-0129,  June  1972. 

M.  Sparks  and  L.  J.  Sham,  "Theory  of  Multiphonon  Infrared  Absorption,"  AFCRL 
Conference  on  High-Power  Infrared  Laser  Window  Materials,  Hyannis,  Massachusetts, 
Oct.  30-Nov.  1,  1972. 

M.  Sparks  and  M.  Cottis,  "Pressure-Induced  Optical  Distortion  in  Infrared  Windows,  " 
AFCRL  Conference  on  High  Power  Infrared  Laser  Window  Materials,  Hyannis,  Massa¬ 
chusetts,  Oct.  30-Nov.  1,  1972. 

M.  Sparks  and  L.  J.  Sham,  "Exponential  Frequency  Dependence  of  Multiphonon  Sum¬ 
mation  Infrared  Absorption,"  Solid  State  Commum  JO,  1451  (1972). 

M.  Sparks,  "Theoretical  Studies  of  High-Power  Infrared  Window  Materials,"  Xonics 
Final  Report,  Contract  DAHC15-72-C-0129,  December  1972. 

M.  Sparks  and  M.  Cottis,  "Pressure-Induced  Optical  Distortion  in  Laser  Windows," 

J.  Appl.  Phys.  44 ,  787  (1973). 

M.  Sparks,  "Stress  "'id  Temperature  Analysis  for  Surface  Cooling  or  Heating  of  Laser 
Window  Materials,"  J.  Appl.  Phys.,  in  press,  September  1973. 

M.  Sparks  and  C.  J.  Duthler,  "Theory  of  Infrared  Absorption  and  Material  Failure  in 
Crystals  Containing  Inclusions, ”  J.  Appl.  Phys.,  in  press,  July  1973. 

M.  Sparks  and  L.  J.  Sham,  "Theory  of  Multiphonon  Absorption  in  Insulating  Crystals,” 
Phys.  Rev. ,  in  press. 

M.  Sparks,  "Short-Pulse  Operation  of  Infrared  Windows  without  Thermal  Defocusing, " 
AodL  Opt. ,  in  press. 
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M.  Sparks  and  L.  j.  Sham,  "Temperature  Dependence  of  Infrared  Absorption," 
submitted  to  Phys.  Rev.  Letters. 


C.  J.  Duthler  and  M.  Sparks,  "Theory  of  Material  Failure  in  Crystals  Containing 
Infrared  Absorbing  Inclusions,"  ASTM  1973  Symposium  on  Damage  in  Laser  Ma¬ 
terials,  Boulder,  Colorado,  May  15-16,  1973. 


M.  Sparks,  "Temperature  and  Frequency  Dependence  of  Infrared  Absorption  as  a 
Diagnostic  Tool,"  submitted  to  Appl,  Phys.  Letters. 


D.  L.  Mills  and  A.  A.  Maradudin,  "Theory  of  Infrared  Absorption  by  Crystals  in  the 
High  Frequency  Wing  of  Their  Fundamental  Lattice  Absorption,  "  Phys.  Rev. ,  in  press. 

A.  A.  Maradudin  and  D.  L.  Mills,  "Temperature  Dependence  of  the  Absorption  Coef¬ 
ficient  of  Alkali  Halides  in  the  Multiphonon  Regime, "  submitted  to  Phys.  Rev.  Letters. 


C.  J.  Duthler  and  M.  Sparks,  "Quasi -Selection  Rule  for  Infrared  Absorption  by  NaCl- 
Structure  Crystals, "  to  be  published. 


C.  J.  Duthler  and  R.  Hellwarth,  "Mechanism  for  Surface  Damage  in  Laser  Window 
Materials,  "  to  be  published. 


M.  Sparks  and  H.  C.  Chow,  "Nonlinear  Infrared  Absorption:  Parametric  Instabilities 
of  Phonons, "  to  be  published. 


L.  J.  Sham  and  M.  Sparks,  "Explicit  Exponential  Frequency  Dependence  of  Multi - 
phonon  Infrared  Absorption,  "  to  be  published. 


The  following  topics  will  be  covered  in  the  final  report: 


Parametric  processes  in  infrared  absorption 

Infrared  absorption  by  imperfections  in  crystals;  ionic  impurities,  disloca¬ 
tions,  band-mode  plus  impurity  mode  absorption,  and  surface  imperfections 

Quasi-selection  rule  for  infrared  absorption 


Explicit  exponential  frequency  dependence  of  multiphonon  infrared 
absorption 


Numerical  evaluation  of  multiphonon  absorption  coefficients 
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Explanation  of  surface  damage  cones  observed  in  high-power  laser 
experiments 

Explanation  of  well  known  anomalies  in  stimulated  Raman  and 
Brilloain  scattering  and  moving  focus  filaments 

Relative  strengths  of  anharmonic  interaction  and  higher-order-dipole 
interaction  in  infrared  absorption 

Ultraviolet- induced  infrared  absorption. 
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B.  THEORY  OF  INFRARED  ABSORPTION  AND  MATERIAL 
FAILURE  IN  CR  YSTALS  CONTAINING  INCLUSIONS 


M.  Sparks  and  C.  J.  Duthler 
Xonics,  Incorporated,  Van  Nuys,  California  91406 


Two  effects  of  inclusions  in  or  on  the  surface  of  infrared- 


transmitting  materials  are  to  increase  the  average  value  of  the  optical 


absorption  coefficient  (3  and  to  cause  localized  heating  that  could 


lead  to  material  failure  at  high-power  levels.  Volume  fractions 


-7  -  8 

as  low  as  10  to  10  of  such  inclusions  can  give  rise  to  a  value 


of  the  optical  absorption  coefficient  j3  of  10"4cm~\  a  typical 


value  of  current  interest.  For  various  types  of  inclusions,  the 


frequency  dependence  of  /3  ranges  from  increasing  as  00  ,  to 


independent  of  to,  to  exponentially  decreasing  with  to.  The  tem¬ 


perature  dependence  ranges  from  independent  of  T,  to  increasing 
as  tP  in  the  high-temperaturt  limit,  where  p  -  2  -  4  typically. 


Simple  expressions  for  the  absorption  cross  section  are  derived 


for  various  cases  of  practical  interest.  The  cross  sections  are 


used  to  derive  expressions  for  /3  for  the  four  cases  of  large  in¬ 


clusions  of  strong  and  weak  absorbers  and  of  small  inclusions  of 


dielectric  and  metallic  particles.  The  material  failure  resulting 


from  local  heating  of  inclusions  is  a  far  greater  problem  in  high- 


intensity  short-pulse  systems  than  in  low-intensity  long-pulse  or 


cw  systems  having  t!  •  same  averse  intensity.  Microsecond  pulses 


with  energy  densities  as  low  as  a  few  joules  per  square  centimeter 


can  cause  material  failure. 


1 
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I.  INTRODUCTION 


The  problem  of  obtaining  highly  transparent  window  materials  for 

1  2 

high-power  infrared  laser  systems  is  of  considerable  interest.  In  par¬ 
ticular,  there  is  great  interest  in  lowering  the  value  of  |3  for  candidate 

materials  such  as  7inSe,  CdTe,  KC1,  KBr,  and  TI  1173  glass  (Ge2g  Sb^  Se^  )• 

-4  -1 

Materials  having  values  of  fi  at  least  as  low  as  10  cm  are  needed.  It  has 

3 

been  shown  that  the  absorption  with  /3  decaying  exponentially  with  frequency  to, 

4 

observed  in  a  number  of  materials,  is  intrinsic  and  results  from  multiphonon 
absorption. 

The  present  investigation  is  concerned  with  another  aspect  of  the  problem  - 

that  of  extrinsic  absorption  by  macroscopic  inclusions  either  in  the  bulk  of  the 

crystals  or  on  their  surfaces.  The  results  of  this  investigation  are  of  practical 

interest  since  it  is  believed  the  current  experimental  values  of  /3  for  all  candi- 

5-7 

date  window  materials  for  high-power  10.  6fim  systems  are  extrinsic.  The 
temperature  and  frequency  dej  .’ndence  of  0  derived  in  Sec.  Ill  should  be  useful 
in  experiments  to  determine  whether  j3  is  intrinsic  or  extrinsic,  especially  as 
improved  materials  become  available. 

Two  aspects  of  optical  absorption  by  inclusions  are  considered.  First, 

the  spatial  average  of  /3  is  increased,  thus  increasing  the  overall  heating  of  the 

material.  Second,  the  local  heating  in  the  region  of  an  inclusion  can  lead  to 

material  failure.  The  overall  increase  in  the  value  of  /3  can  cause  either 

irreversible  system  failure,  such  as  thermally  induced  fracture,  or  reversible 

5-7 

failure,  such  as  thermal  defocusing  of  the  laser  beam  by  the  heated  window. 
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tor  sufficiently  large  concentrations  of  inclusions,  the  localized  heating  can  also 
cause  considerable  optical  degradation  of  the  beam.  In  the  present  study  it  is  as 
sumed  that  the  concentration  of  inclusions  is  so  small  that  this  localized- heating 
type  of  optical  distortion  is  negligible.  However,  it  should  be  mentioned  that  scat¬ 
tering  may  be  considerably  greater  at  high-power  levels  than  at  low  levels  as  a 
result  of  the  scattering  by  the  heated  host  material  near  the  inclusion.  This  effect 
she  vJ  be  greatest  for  scattering  near  the  forward  angle. 

Local  material  failure  at  discrete  inclusion  sites  has  been  observed  and  studied 

by  others.  ^  In  these  treatments,  which  were  concerned  with  metallic  inclusions 

2 

in  glass  hosts,  the  absorption  cross  section  aabg  =  e^tra  was  used,  where  is 
the  bulk  emissivity  of  the  inclusion  material  and  a  is  the  inclusion  radius.  The 
resulting  errors  in  <j^  can  be  large,  especially  for  inclusions  with  diameters 
less  than  the  laser  wavelength.  The  heating  in  transparent  hosts  of  both  dielectric 
and  metallic  spherical  inclusions  will  be  considered  as  a  function  of  inclusion  radius 
and  laser  pulse  length,  using  more  realistic  models  for  the  absorption  cross  section 
and  for  the  heat  transfer  from  the  inclusion  to  the  host. 

Absorption  is  of  greater  interest  than  scattering  in  the  study  of  high-power  in¬ 
frared  windows.  A  value  of  0  =  10'4  cm'1  for  the  absorption  coefficient  can  cause 

sufficient  heating  for  the  system  to  fail,  for  example,  by  thermal  defocusing  or  by 

-4-1 

window  fracture.  But  an  equal  amount  of  scattering  Pscat  =  10  cm  may  be  toler¬ 
able.  Thus,  scattering  wili  be  neglected  here,  except  to  mention  that  observation  of 
the  associated  scattering  may  help  to  identify  the  absorption  mechanism  in  some 
cases.  Winsor10  has  shown  theoretically  that  scattering,  especially  in  conjunction 
with  total  internal  reflection  at  the  host -crystal  boundaries,  may  increase  the  mea 
sured  value  of  /3  by  increasing  the  path  lengths  of  the  rays  in  the  crystal. 


11 
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II.  ANALYSIS  OF  ABSORPTION  CROSS  SECTIONS 


In  this  section,  the  absorption  cross  section  of  an  individual  spherical 

inclusion  of  radius  a  will  be  considered.  The  absorption  cross  section  is  not 

2 

generally  equal  to  the  geometrical  cross  section  ira  .  For  ka  «  1,  where 

k  =  2 it /A,. ,  with  A.,  the  wavelength  of  the  radiation  in  the  host  material,  the 
H  i  i 

2 

value  of  typically  is  small  (O^g  Fa  )•  1°  ^ie  case  of  ka  »  1  and 

|  e  |  »  1,  where  €  =  ,  with  ^  and  the  dielectric  const?nts  of 

the  inclusion  and  host,  respectively,  the  reflection  at  the  surface  of  the  sphere 

2  _ 

is  great,  which  again  makes  cr^  «  IT  a  .  E  tact  absorption  cross  sections 

11  12 

for  spheres  of  arbitrary  size  can  be  obtained  from  the  classic  result  of  Mie.  * 
The  Mie  solutions  are  complicated  in  general,  but  reduce  to  simple  results  in 
the  limits  ka  »  1  and  ka  «  1 . 

For  small  spheres  ( ka  «  1 ) ,  the  Mie  series  is  well  approximated  by  the 
first  term,  which  yields 


12  ka 
(  +  2)  + 


for  ka  «  1  , 


(2.  n 


where  €  =  c  ^  t-  i  €j  .  Two  limiting  cases  of  (2.  1)  are  of  interest.  For 
Cj  «  ,  which  is  satisfied  for  nonmetals  at  frequencies  not  too  near  the 

fundamental  resonance  frequency  or  the  high- frequency  absorption  edge,  j3j 
and  (j  are  related  by  the  expression 


PT  =  2  na  k  a  -f-  , 


(2.2) 
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where  n  =  n^  +  inj .  Using  this  result  to  eliminate  e^k  in  (2. 1)  gives 


12  n  jrt  n 

Tabs  =  - <0ia,ffa  > 

(fft  +  2T 


Cc9  <<C  € 


(2.3) 


The  second  limiting  case  of  (2. 1)  is  that  of  small  metallic  inclusions.  The 

13 

Drude  expression  for  the  dielectric  constant  is 


e  =  e. 


eH(to2  +  itoT ) 


(2.4) 


wheTe  F  is  the  electron  relaxation  frequency  (often  written  as  1  /  t),  em  is  the 

2  1/2 

contribution  to  €  from  the  core  electrons,  and  to^  =  (4  7lNe  /m)  is  the  plasma 

frequency,  which  has  a  typical  value  of  top  =  5  x  1015  sec  *.  At  10. 6fim, 

;o  2:  1.  9  X  1014  sec"1.  There  are  two  contributions  to  the  relaxation  frequency  T: 

13 

r  =  Fgu  +  Tgu  ,  where  the  bulk  contribution  has  a  typical  value  of  rBu 
^  5  X  lO1^  sec  1 .  The  value  of  the  surface  scattering  contribution14  Tgu  is  ~  Vp/a  , 

g 

where  the  Fermi  velocity  Vp  has  the  value  vp  -  10  cm/sec  for  many  metals. 
With  Tpu  -  5  x  10-  sec  1  and  Vp  2:  10^  cm/sec,  >  Tpu  for  a  <  a^,  -  200  & . 

For  typical  metals  at  10.  6}im,  to  »  |  to  TitoT  |  ,  and  (2.4)  gives 

? 


f  f  s; - 1 -  - - 1 - 

H  w2+itor  to2  +  r2 


eR  +  2  "  e(R  *  (2,5) 


Both  the  real  and  imaginary  parts  of  the  dielectric  constant  are  large  in  magni¬ 
tude  for  small  particles  of  typical  metals  at  10.  6  Jim  .  Substituting  (2.  5)  into 
(2.1)  gives 
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°abs  =  12  fH 


3/2  {jy  r  a 


2  c 


— —  7f  a.  ;  CO  »  co -  +  ico  £  , 


(2.6) 


for  metals  with  ka  «  1 . 

Next  consider  the  case  of  large  spheres  (ka  »  1 ).  Using  geometrical 
optics  and  the  identity  l+x  +  x2  +  *"  =  (1-x)"1,  where  x  =  R.exp(-^d), 
with  R.  the  internal  reflection  coefficient  and  d  the  distance  the  ray  travels 


in  traversing  the  sphere  once,  gives 


-Mr1 


°ab«  =  *a 


f  9  (  "MW  ""ld\ 

d(cos26)(l-R)(  1- e  j  ^ i -Rie  j 

n 


(2.7) 


where  R  =  ( |  rp  |2  +  |  rn  |2  )/2 ,  with  rp  and  rn  the  Fresnel  reflection  coeffi¬ 
cients  for  the  two  polarizations,  d  the  distance  that  the  refracted  ray  travels 
through  the  sphere,  and  6  the  angle  of  incidence. 

There  are  two  limiting  cases  of  (2.7)  of  interest.  First,  for  »  1 , 
which  is  typically  satisfied  for  metals  and  strongly  absorbing  dielectrics,  (2.7) 


yields 


°ab>,  “  *a  <‘-<R>> 


(2,8) 


where  the  axerage  reflection  coefficient  (R)  is  defined  as 


<  R)  = 


R  d(cos^8 )  . 


For  j3ja  «  1 ,  in  the  limit  of  small  index  of  refraction,  using  |  n^  -  1  |  «  1 , 


R  2*  R  2:  o ,  and  d  S  a  cos6  in  (2. 7)  yields 
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abs 


4  r  —  n  3 
3  ^ 


(2.9) 


As  n^  departs  from  1 ,  a^g  first  increases  slightly  for  small  |  n  -  1  |  and 
eventually  goes  to  zero  when  R  goes  to  one.  Notice  that  (2.9)  has  the  same 
functional  dependence  on  a  as  (2.3),  but  with  a  slightly  different  coefficient. 


For  metallic  inclusions,  a  schematic  illustration  of  the  absorption  efficiency 
a^bg/^3^*  obtained  by  sketching  the  results  (2.  6)  and  (2.  8),  is  shown  in  Fig.  1. 
The  dashed  line  represents  the  asymptotic  values  obtained  from  (2.  6)  with 
r  ~  a_1  ( for  a  «  ap  ),  from  (2.  6)  with  T  ~  a°  (  for  ap  «  a  «  k  1 ),  and 
from  (2.  8)  (for  a  »  k"1 ).  The  extrapolated  dashed  curve  from  (2.  o)  with 
a  »  ap  intersects  the  dashed  curve  from  (2.  8)  at 


*k  = 


top  c(l-  a>) 

372“! - 


12  € 


H 


to 


’Bu 


(2.  10) 


which  occurs  near  k"1  for  many  metals  in  the  infrared.  The  solid  curve  sche¬ 
matically  illustrates  the  results  in  the  intermediate  regions  a  ~  ap  and  a  ~  a^ . 
The  dashed  curve  will  be  sufficient  for  the  order-of-magnitude  estimates  of 
overall  absorption  and  failure  intensity  in  Secs.  Ill  and  IV. 


In  dielectric  inclusions,  fya  «  1  typically  is  satisfied  except  for  very 

4 

strong  absorbers  with  rather  large  radii.  (For  example,  ,3ja  =  1  for  /3j  =  10 
cm"1  and  a  -  l^m.)  For  fya  «  1,  the  absorption  cross  section  is  given  by 
(2.3)  for  ka  «  1  and  by  (2.9)  for  ka  »  1.  Since  these  two  limiting  results 

3 

both  are  of  the  form  CTabg  ~  jSja  ,  with  only  slightly  different  coefficients,  the 

3 

approximation  CJ^g  =  A  fya  ,  with  A  the  average  of  the  two  coefficients,  will 
be  used  for  all  values  of  ka . 
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The  absorption  efficiency  for  dielectric  inclusions  is  qualitatively  similar 

to  that  sketched  in  Fig.  1  for  metallic  inclusions,  with  two  exceptions.  First, 

2 

for  ka  «  1,  (2.  3)  indicates  that  “  a  »  and  the  constant  region  at 

small  a  does  not  occur  for  dielectrics.  Second-,  the  extrapolated  linear  region 

crosses  the  large  a  asymptotic  region  near  rather  than  k  *.  This  can 

be  seen  from  (2.  8)  and  (2.  9),  which  are  valid  for  a  >  1  and  a  <  ^  1 , 

respectively.  In  some  cases  with  n^-1  and  n^  small,  the  exact  Mie  solution 
2 

yields  oabg  >  IT  a  near  ka  =  1 . 

For  large-bandgap  semiconductors,  such  as  A1P,  SiC,  and  ZnS,  the  dielectric- 
inclusion  results  above  can  be  applied.  The  absorption  by  small-bandgap  semi¬ 
conductors  is  more  complicated  than  that  by  dielectrics  and  metals  in  general, 
and  will  not  be  considered  explicitly.  Problems  can  occur  involving  tempera¬ 
ture  dependence  of  the  electrical  conductivity,  increased  absorption  caused  by 
free  carriers  that  are  created  in  the  absorption  process,  and  the  resulting 

4 

thermal  runaway. 
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III.  ABSORPTION  EFFICIENCIES  FOR  VARIOUS  TYPES  OF  INCLUSIONS 


The  above  cross  sections  will  now  be  used  to  calculate  properties  of  the 
absorption  coefficient  for  crystals  containing  various  types  of  inclusions.  The 
volume  fraction  f  of  inclusions  required  to  make  0  =  10  4  cm  ,  a  value  of 
current  interest,  will  be  determined. 

Consider  a  sample  consisting  of  a  nonabsorbing  host  material  of  dielectric 
constant  €„  (real)  containing  N.  inclusions  per  unit  volume,  each  inclusion 
having  absorption  cross  section  o^  .  Multiple  scattering  will  be  neglected  - 
a  reasonable  approximation  for  the  present  case  of  small  impurity  concentrations. 
The  well  known  result  for  the  absorption  coefficient  0  of  the  sample  is  then 


0  = 


CT  .  N. 
abs  I 


(3.1) 


For  small  dielectric  inclusions  ,  a^g  is  given  by  (2.1),  which  when  substi 
tuted  into  (3. 1)  yields 

9n« 

0  S!  - 2 - 2 

'V2)  +CJ 

where  the  factor  in  the  bracket  typically  has  a  value  near  unity.  The  absorption 
coefficient  in  (3.2)  is  independent  of  a  ,  but  is  strongly  temperature  and  fre¬ 
quency  dependent  in  general.  If  0j  is  controlled  by,  say,  the  n-phonon  summa¬ 
tion  process,  then  0  ~  Tn_1  in  the  high-temperature  limit,  and  0  decays 
exponentially  with  frequency.  Using  0j  =  10  to  104  cm  for  strongly  absorbing 
inclusions,  (3.2)  gives  0  =  lO^cm'1  for  volume  fractions  of  f  =  10  to  10  . 


0  f  ,  for  ka  «  1,  €jj  ^  1  »  (3.2) 


mtmmtmammmm  w 
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In  the  case  of  large  dielectric  inclusions,  there  are  two  possibilities. 
First,  for  ka  »  1  and  ^a  «  1 ,  using  the  geometrical-optics  absorption 
cross  section  (2.9)  yields  jS  «  0jf  with  a  numerical  factor  near  unity,  as  in 
(3.2).  For  ka  »  1  and  fi^a  «  1,  (2.  8)  and  (3. 1)  give 


P=  7T  a  (1-  (R) )  Nj  =  3  ( 1  -  (R)  f  /4  a  , 


ka  »  1  ,  fya  «  1  .  (3.3) 


The  absorpti  ‘'efficient  in  (3.  3)  is  proportional  to  a  for  a  given  value 

of  f  and  is  generally  temperature  independent,  except  near  the  reststrahl  region 

-2 

where  (R)  is  strongly  temperature  dependent.  For  1  -  (R)  -  1  and  a  -  10 
to  10  ^  cm ,  (3.  3)  gives  =  10  ^  cm  for  f  in  the  range  from  10  ^  to  10  . 

For  large  metallic  inclusions,  the  value  of  ( 1  -  (R>)  in  (3.  3)  is  small 

(  <  1/10),  since  the  reflectivity  of  metals  in  the  infrared  is  great.  As  a  result, 

-4 

the  volume  fraction  of  inclusions  for  jS  =  10  cm  is  increased  by  a  factor  of  at 
least  ten  over  the  corresponding  dielectric  case. 

Next  consider  the  case  of  small  metallic  inclusions.  From  (2.  6)  and  (3.  1), 


(3.4) 


This  expression  shows  that  0  increases  quadratically  with  frequency.  For 

temperatures  greater  than  room  temperature,  the  electron  relaxation  frequency 

Tgu  and  hence  /3  increases  linearly  with  temperature  having  a  typical  fractional 
-2  3/2 

increase  of  10  per  degree  Kelvin.  With  eH  =  10  ,  to  =  cc _/25  ,  and 

F  =  5  x  10^  sec  \  an  inclusion  volume  fraction  of  f  =  4  x  10  results  in 

-4  -1 

0  =  10  cm  . 
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The  absorption  by  small  metallic  inclusions  has  a  peak  that  typically  lies 

13 


in  the  ultraviolet  or  visible  region.  The  value  of  j3  at  this  peak  is  much 
greater  than  the  value  in  the  infrared.  Thus,  ultraviolet  and  visible  measure¬ 
ments  of  absorption  or  scattering  can  be  used  to  verify  the  source  of  absorption 
by  small  metallic  particles  in  the  infrared. 


Consider,  as  an  example,  small  potassium  spheres  in  KC1  or  KBr 


F  centers  can  be  transformed  to  colloidal  potassium  (small  spheres)  by  heating 


1S-17 

the  crystal.  The  transformation  is  enhanced  by  ultraviolet  radiation. 


The  wavelengths  of  the  peaks  for  the  small  potassium  spheres  in  KC1  and 


13 

KBr  are  0.730  and  0.  770  Jim,  respectively.  In  the  visible  and  ultraviolet, 


r2  «  a;2,  and  (2.  1)  and  (2.4)  give,  with  k  =  n^  to/c 


to. 


a.' 


€  -  C. 


+  i 


nl_ 


€h  u 


j3  s 


rV/U 


2  2  ^  2  2 
(U  ■  OJ0  )  +  (T  o-'q  /to) 


(3.5) 


wher.; 


to. 


to. 


o  3  4 

9nH  w0 


2  €H+€I*> 


pk 


2  T 

c  to  r 

p 


Notice  that  for  fairly  narrow  lines  (T  ^  to^/lO),  |S  2:  0  ^  at  the  peak  a 


to  -  cOq  ,  and  T  is  the  full  line  width  between  the  points  $  =  y  |3  ^  . 
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Dividing  j3^  by  the  infrared  absorption  coefficient  in  (3.  4)  gives 


ftpk 


60 


0 


P 


IR  ^IR 


2r2 


(3.6) 


14  -1 

For  XQ  =  2  irc/coQ  =  0.75fim,  XJR  =  10.6fim,  and  T  =  5  X  10  sec  , 
(3.6)  gives 


in*-  =  6.4  X  103  .  (3-7) 

*10.6 

For  |3  =  lO-4  cm"1,  (3.7)  gives  =6.4  cm-1,  which  would  produce  a 
visibly  colored  crystal.  Uncolored  KBr  or  KC1  with  fi  ,  <  10  cm  would 
have  fitn  .  <  1.6  X  10_4cm_1,  and  the  contribution  from  potassium  colloids 
to  the  infrared  absorption  would  be  small.  Colored  crystals  could  have  a 
greater  contribution  to  0  from  this  source.  It  should  also  be  mentioned  that 
impurities  in  the  form  of  F  centers,  which  give  rise  to  strong  absorption  in 
the  visible  region,  may  not  give  rise  to  detectable  absorption  in  the  infrared. 
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IV.  MATERIAL  FAILURE  FROM  LOCAL  KEATING 


The  heating  of  macroscopic  inclusions  can  give  rise  to  localized  regions 
of  high  temperature  that  can  cause  material  failure  when  the  intensity  is  great. 
This  is  not  a  nonlinear  effect,  but  is  usually  important  only  at  relatively  high 
intensities.  For  times  short  with  respect  to  a  characteristic  time  for  heat  to 
diffuse  a  distance  a,  very  roughly  speaking,  most  of  the  energy  absorbed 
by  the  inclusion  remains  in  the  inclusion.  Thus,  the  temperature  in  the  inclu¬ 
sion  increases  linearly  with  time.  For  times  large  with  respect  to  the  char¬ 
acteristic  time,  part  of  the  absorbed  energy  has  diffused  into  the  host  crystal, 
and  the  temperature  rise  in  the  inclusion  is  considerably  less  than  the  value 
obtained  by  neglecting  diffusion.  Thus,  a  given  amount  of  energy  will  cause  a 
greater  temperature  rise  if  it  is  deposited  in  a  time  that  is  short  with  respect 
to  the  characteristic  time  than  if  deposited  in  a  time  long  with  respect  to  the 
characteristic  time ;  the  local  heating  of  macroscopic  inclusions  is  a  more 
severe  problem  in  high -intensity  short-pulse  systems  than  in  low-intensity 
long-pulse  or  cw  systems  of  equal  average  intensity. 

The  criterion  for  failure  of  the  window  material  depends  on  the  details  of 
the  laser  system  and  the  type  of  material  and  inclusions.  Since  there  is  no 
universal  criterion,  it  will  be  assumed  th  :  a  temperature  rise  of  1,000K 
constitutes  failure.  This  is  a  reasonable  choice  for  the  following  reasons: 

This  temperature  is  approximately  the  correct  value  for  melting  temperatures 
and  fracture-inducing  temperatures.  The  latter  have  typical  values  of  the  order 

IQ,  £0 

of  '  Oj/ttE,  where  cr^  is  the  material  strength,  a  the  linear  thermal  expan¬ 
ds  -5 

sion  coefficient,  and  E  the  Youngs  modulus.  For  =  1CT  psi ,  a  =  10  and 
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E  =  IQ7  psi ,  the  temperature  corresponding  to  fracture  is  1.000  K.  Heats  of 
fusion  have  typical  values  corresponding  to  several  hundred  degrees  Kelvin. 

At  1, 000  K  above  ambient  temperature  in  materials  which  do  not  melt,  the 
ionic  diffusion  may  be  important.  Although  order-of-magnitude  accuracy  in 
temperature  usually  is  not  sufficient,  order-of-magnimde  accuracy  of  intensi¬ 
ties  corresponding  to  failure  is  all  that  can  be  expected  at  present,  and  this  is 
often  adequate.  Since  the  present  interest  is  in  this  failure  intensity  If  .  and 
the  failure  temperature  Tf  is  linearly  related  to  If,  the  value  Tf  =  1.000  K 
should  be  sufficient  for  present  purposes. 


Two  important  features  of  high-power  laser-window  failure  are  that  failure 
of  the  weakest  part  of  the  window  can  constitute  system  failure,  and  that  fatigue 
and  other  multiple -pulse  effects  must  be  considered  when  repeated  pulses  must 
be  withstood.  Concerning  the  former,  a  single  inclusion  in  a  window  conceivably 
could  cause  failure.  As  an  example  of  the  latter,  in  a  single-pulse  measure¬ 
ment.  a  laser  glass  conceivably  could  melt  locally  and  recrystallize  without 
leaving  detectable  damage.  For  a  window  in  a  pulse-operated  system,  the  local 
absorption  coefficient  could  be  changed  by  the  high  temperature  associated  with 
the  first  pulse  or  the  first  n  pulses,  thus  causing  increased  absorption  in  sub¬ 
sequent  pulses  with  eventual  failure. 


Bloembergen21  has  suggested  that  local  field  enhancement,  such  as  that 
occurring  at  the  edge  of  a  crack  tn  a  material,  may  give  rise  to  local  intensi¬ 
ties  up  to  100  times  greater  than  the  nominal  external  intensity.  Thus,  if  an 
inclusion  is  in  the  high-intensity  region,  this  effect  could  lower  the  failure  in¬ 
tensities  calculated  below  by  a  factor  of  the  order  of  100.  Local  field  enhance¬ 
ment  can  also  occur  when  one  inclusion  is  at  the  focal  point  of  anothe'.  Since 
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focusing  is  limited  by  diffraction  when  ka  <  1 ,  the  focusing  by  large,  weakly 
absorbing  inclusions  is  most  severe.  Focusing  by  surface  imperfections  also 
could  occur. 

The  temperature  rise  below  the  melting  point  can  be  calculated  simply  for 
the  following  model.  The  spherical  inclusion  of  radius  a  is  assumed  to  have 
temperature- independent  values  and  Kj  of  heat  capacity  per  unit  volume 

and  thermal  conductivity.  The  host  crystal  is  assumed  to  have  temperature- 
independent  values  Cj^  and  K^.  The  boundary  between  the  two  is  assumed  to 
be  thermally  perfect ;  that  is,  there  is  no  thermal  impedance.  Heat  absorption 
by  the  host  crystal  is  assumed  negligible.  The  relaxation  time  required  to 
transfer  energy  to  heat  from  the  modes  that  absorb  energy  is  assumed  to  be 
much  shorter  than  the  laser  pulse  duration. 

Simple  solutions  to  the  heat -flow  equation 

-KV2T  +  C  -  S  (4.1) 

0  t 

will  be  derived  for  a  series  of  limiting  times  for  the  cases  of  uniform  heat  gen¬ 
eration  within  the  volume  of  the  inclusion  and  of  uniform  heat  generation  over 
the  surface  of  the  inclusion.  In  (4. 1),  S  is  the  rate  at  which  heat  is  generated 
per  unit  volume. 

First  consider  the  case  of  spatially  uniform  heat  generation  within  the 
inclusion  at  the  rate 

S  =  3  aabs  I/4  tra3  (4.2) 

per  unit  volume,  where  I  is  the  incident  intensity.  This  applies  for  dielectric 
inclusions  with  ^a  «  1  or  for  metallic  inclusions  with  skin  depth  6  >  a. 
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Roughly  speaking,  the  thermal  time  constant 


T  =  C  a  /  4  K 


(4.  3) 


is  the  time  required  for  heat  to  diffuse  a  distance  a  in  either  the  inclusion  or 
the  host  when  the  appropriate  values  of  C  and  K  are  used.  Subscripts  I  and  H 
will  denote  the  values  of  T  in  the  inclusion  and  in  the  host  crystal,  respectively. 
For  short  times  t  «  ,  the  diffusion  of  heat  out  of  the  inclusion  is  negligible. 

The  term  KjV  2  T  in  (4. 1)  is  then  negligible.  The  temperature  at  the  center  of 
the  inclusion  is,  from  (4. 1)  and  (4.  2), 


Tc  =  3oabs“/4"a  CI  ' 


for  t  «  T. 


(4.4) 


for  T  defined  as  zero  at  the  time  the  laser  is  turned  on  (t  =  0 ) . 

For  long  times  t  »  Tj,  T^,  equilibrium  is  reached  with  the  host  material 

conducting  heat  away  from  the  inclusion  at  the  same  rate  that  it  is  generated 

within  the  inclusion.  In  this  case  the  temperature  is  obtained  by  solving  (4. 1) 

2  2  2 

with  dT/d  t  =  0  in  both  the  inclusion  and  the  host.  Using  V  r  =  6  ,  V  (1/r)  =  0 
for  r  >  a ,  and  the  boundary  conditions  that  the  temperature  and  the  heat  flow 
be  continuous  at  the  inclusion  surface,  yields  the  steady-state  value  at  the  in¬ 


clusion  center. 


r  -  3  gabs 1 

c®  8  IT  a  K  ,, 


(4.5) 


For  short  and  long  times,  the  temperatures  are  approximated  by  (4.4)  and 
(4.  5),  respectively.  These  are  shown  as  dashed  lines  in  Fig.  2,  while  the  actual 
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temperature  is  sketched  as  the  solid  lire,  The  extrapolated  short-time  curve 


intersects  the  steady-state  value  at  the  time  r  ^  =  (2Kj/Ke^)  Tj,  which  is 
found  by  equating  (4. 4)  and  (4.  5). 


Next  consider  the  case  in  which  the  heat  is  generated  near  the  surface  of 


the  inclusion,  rather  than  uniformly  throughout  its  volume.  Such  is  the  case 


for  metallic  inclusions  in  dielectric  hosts,  where  a  typical  value  for  the  skin 


depth  6  is  40  A  (Cu  at  X  =  10.  6  Jim),  or  for  dielectric  inclusions  with 


fi^a  »  1,  where  1/ is  of  the  order  of  10  cm  for  strong  absorption. 


For  spherical  inclusions  with  radius  a  »  6 ,  it  is  assumed  that  heat  is 


generated  uniformly  within  a  layer  of  thickness  6  over  the  entire  surface  of 


the  inclusion.  This  is  a  good  approximation  for  a  «  X  because  the  electric 


field  is  nearly  constant  over  distances  of  the  order  of  a.  In  general , 


there  will  be  local  hot  spots  over  the  surface  --  not  only  on  the  front  surface, 


but  also  on  the  rear  surface. 


There  are  three  characteristic  times  of  interest.  First, 


T6  -  4^6^/ffKj 


(4.  6) 


is  roughly  the  time  in  which  heat  diffuses  out  of  the  skin  depth,  assuming  negli¬ 


gible  diffusion  into  the  host  for  small  time,  since  (dielectric)  «  Kj  (metal). 


Second, 


Ta  =  4CIaz/9irKI 


(4.7) 


is  the  time  in  which  heat  diffuses  from  the  inclusion  surface  to  the  center, 


roughly  speaking.  Third, 


1 
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tH  "  CHa  /3KH 


(4.  8) 


is  roughly  the  time  in  which  heat  diffuses  a  distance  equal  to  the  radius  of  the 
inclusion  into  the  host.  The  values  of  the  numerical  coefficients  in  (4.6)-(4. 8) 
are  chosen  for  later  convenience. 

For  t  «  Tg  and  6  «  a ,  the  spherical  shape  of  the  inclusion  surface 

is  not  important,  and  the  solution  to  the  simpler  problem  of  heat  generation  in 

a  thin  plane  slab  can  be  used.  In  this  case,  the  temperature  is  obtained  from 
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the  Laplace -transform  solution  of  (4. 1).  The  solution  to  the  transform  equation 


at  the  inclusion  surface  is 


T(x  =  6 ,  p)  =  £ 


CtKt 


I  /C^+/CHKt 


X/ 4. 

where  q^  =  (pC^/K^)  .  Taking  the  in/erse  transform  and  keeping  only  the 

dominant  term  for  t  small  yields  the  temperature  at  the  inclusion  surface 


1  *aabs  *aabs 

T  =  -4 - 1  1  ■  —  — t  ,  t  «  Tc  .  (4.9) 

&  I  (/C^  +/ChKh  )  47ra“6  4ffa*6Cj 


In  (4.  9)  and  in  the  equations  below,  the  approximate  equalities  are  valid  for 
CjKj  »  .  For  «  T  «  Ta,  the  heat  generation  in  a  plane-slab 

problem  can  be  approximated  by  a  delta-function  source  at  x  =  0.  Using  the 
method  of  Laplace  transforms  and  keeping  the  dominant  term  for  small  t 
yields 


T  = 

1  O 


* CTabs  1/2 


S  JT  (/CjjKjj  +/C^Kj)  4  77a2 


r6  «  t  «  T  .  (4.  10) 
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When  t  »  ra,  me  temperature  Inside  the  spherical  inclusion  reaches 
spatially  uniform  value.  The  present  problem  then  is  equivalent  to  uniform 
heat  generation  within  an  inclusion  of  infinite  conductivity  and  the  previous 
result  (4.  4)  can  be  used  to  obtain 


T  = 

\s 


3I(W 


4ira3C. 


T  «  t  «  t 


H 


(4. 11) 


For  t  »  th  ,  the  temperature  inside  the  spherical  inclusion  reaches 
an  equilibrium  value  that  is  spatially  uniform.  This  case  is  again  equivalent 
to  uniform  volume  heating  within  an  inclusion  of  infinite  conductivity,  and  (4.  5) 


yields 


la 


abs 


S°°  4rraK 


H 


t  »  t 


H 


(4. 12) 


The  temperature  at  the  surface  of  the  spherical  inclusion  is  sketched  in 
Fig.  3  for  the  various  time  regimes  in  (4.  9)-(4. 12).  The  extrapolated  low- 
temperature  linear  time-dependent  section  of  the  curve  intersects  the  t  ^2 
curve  at  the  time  tg  =  Tg  ,  at  which  time  the  temperature  is 


T6  = 


I6trabS  I 

~T~r  k 


16  a 


abs 


7T"a"  *1  (/C^  +  ZC^,)  tt2  a2  Kj 


(4. 13) 


1  /9 

At  t  =  ta,  the  extrapolated  t  7  curve  intersects  the  second  linear  region. 
The  temperature  at  this  intersection  has  the  value 


T  = 
a 


la  u  CT 
abs  I 


I  a 


abs 


3tt  a(/C^K^  +*/CjKj  )  3  tt  aK. 


(4. 14) 


29 


Figure  .3.  Temperature  at  the  surface  of  a  spherical  inclusion  as  a  function 


of  time  in  the  case  of  surface  heating, 
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turn  tv 


and  the  value  of  t  is 

U 


t  = 
a 


2  2 
4aZCj 


9l7(v/CnKH  +  n/CjKj  ) 


—  4  Ta 


(4.  15) 


At 


the  time  t  =  (C^/C^,)  TH  *  the  extrapolated  linear  curve  crosses  the 
equilibrium  value  given  in  (4. 12). 


Consider  the  effect  of  inclusion  size  on  the  failure  intensity  If .  Hie  various 


time  regimes  in  Fig.  3  depend  on  inclusion  size  in  such  a  wav  that  if  the  pulse 
length  is  fixed,  long  times  in  the  figure  are  associated  with  small  inclusions. 


For  example,  the  time  t  at  which  the  temperature  reaches  equilibrium  is 

6C[ 


equivalent  to  an  inclusion  of  radius  agq  -  (  3  K^t/  Cj)  for  a  pulse  of  dura 


tion  t.  The  temperature  is  then  given  by  (4. 12)  for  a  <  agq.  Similarly,  the 


time  t  in  Fig.  3  corresponds  to  a  size  at  s:  (9tt  K1t/4Cj)  ,  and  the  tern- 

perature  is  given  by  (4. 11)  for  aeq  <  a  <  at  .  Next,  the  time  tfi  is  independent 

c«  size  and  is  determined  by  the  skin  depth  and  thermal  properties  of  the  inclusion. 

,-13 


In  the  case  of  metallic  inclusions  ,  t&  is  typically  of  the  order  of  10  sec, 
which  is  much  less  than  most  laser  pulse  durations  of  interest.  Hence,  the  first 
Unear  region  in  Fig.  3  given  by  (4.  9)  does  not  occur,  in  general,  and  for  a  > 
the  temperature  is  given  by  (4. 10).  To  determine  the  temperature  in  a  pulsed 
system  as  a  function  of  the  inclusion  radius,  the  dependence  of  aabg  on  a  in 
(4.  9)-(4.  12)  must  be  included.  In  the  case  of  metallic  inclusions,  the  dashed 
curve  of  Fig.  1  is  used  as  an  approximation  for  aabg .  For  a  given  type  of  in¬ 
clusion  in  a  particular  host,  the  radii  a^  and  a^  at  which  the  functional  de¬ 
pendence  of  the  cross  section  changes  are  independent  of  the  pulse  duration, 
while  the  radii  a,  and  apn  are  both  proportional  to  t1/2  and  decrease  with 

L  C'-j 

decreasing  pulse  length. 
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There  are  six  combinations  of  the  sequence  of  ap,  ak>  a{ ,  agfj  which  can 
occur  for  different  pulse  lengths.  The  temperature  as  a  function  of  inclusion 
size  is  sketched  for  each  of  these  cases  in  Fig.  4,  with  the  pulse  duration  de¬ 
creasing  in  going  from  A  to  F.  For  a  given  type  of  inclusion,  the  curves  C 
and  D  cannot  both  occur.  The  applicable  case  will  depend  on  the  relative  size 
of  ratios  at/a  and  a^/ap. 

Consider  the  example  of  metallic  inclusions  with  the  laser  wavelength 
equal  to  10.  6  Jim .  Using  the  typical  values  u:p  =  5  x  1015  sec  l,  T]]u  =  5  x  10 
sec'1,  (1  -  (R) )  =  0.1,  and  €  ,.3/2  =  10  in  (2.10)  yields  ak  =  4pm.  A  typical 
value  of  ap  is  200  A .  Using  CH  a  Cj  a  2J/cm3K,  Kj  =2W/cmK  and 
K  =  10"2W/cmK,  these  typical  values  indicate  that  aeq  =  ak  when  t=  10 
sec,  at  =  ak  when  t  =  2  x  10'8  sec,  aeq  =  ap  when  t  =  3  X  10  10  sec,  and 
a  =  ap  when  t  =  6  x  10'13  sec.  Hence,  curve  A  applies  for  t  <  10  sec  , 
that  is,  in  the  range  of  cw  or  millisecond  pulses.  Curve  B  applies  to  micro¬ 
second  pulses,  curve  C  or  D  to  nanosecond  pulses,  curve  E  to  picosecond 
pulses,  and  curve  F  to  subpicosecond  pulses. 

The  maximum  temperature  for  a  microsecond  duration  pulse  occurs  for 

a  <  a  <  a,  in  Fig.  4B,  where  a  =  ljm  using  the  above  parameters, 
eq  k  eq 

For  radii  in  this  range,  the  temperature  is  given  by  (4. 11),  which  is  evalu¬ 
ated  at  a  =  ak  using  the  absorption  cross  section  (2.  8).  With  the  above  values 
in  (4. 11),  the  failure  temperature  of  1,000K  occurs  for  a  pulse  energy  of 
3J  /cm2  with  micron-size  inclusions. 

With  a  nanosecond  duration  pulse,  the  maximum  temperature  again  occurs 

near  a  =  a^  in  Fig.  4D  and  is  determined  from  (4. 10).  In  this  case,  the 

2 

failure  temperature  occurs  at  a  pulse  energy  of  2J  /cm  . 
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Figure  4.  Temperature  at  the  surface  of  a  metallic  inclusion  as  a  function  of 
radius  for  various  pulse  durations. 


2 

Next  consider  a  cw  system  with  I  =  300W/cm  ,  an  intensity  which  is 
equal  to  the  average  intensity  of  a  repetitively  pulsed  system  with  a  pulse 
energy  of  3J  /cm2  and  a  100  pulse/sec  repetition  rate,  hi  this  cw  case,  the 
temperature  of  an  inclusion  with  a  =  a^  is  determined  by  (4. 12),  which  yields 
a  temperature  rise  of  less  than  1  K  above  ambient.  Hence,  the  pulsed  system 
would  fail  on  a  single  pulse,  while  the  cw  system  of  the  same  average  intensity 
would  have  a  negligible  temperature  rise. 

At  laser  wavelengths  other  than  10.  6  pm,  the  various  curves  of  Fig.  4 

correspond  to  pulse  durations  different  from  those  listed  above.  Consider  the 

case  of  platinum  inclusions  in  a  glass  host  with  A  =  l(im.  Fitting  the  Drude 

expression  (2.4)  for  the  dielectric  constant  of  platinum  to  tabulated  values  of 

the  refractive  index23  yields  the  values  for  the  parameters  and  Ij^, 

which  together  with  =  3  when  substituted  into  (2. 10)  gives  ak  =  400  A 

2 

for  platinum  inclusions  in  a  glass  host.  In  this  case,  using  Kj  =  0.7W/cm  K, 
Kh  =  1.3x  lO-2  W/cm  K ,  Cf  =2.8J/cm3K,  and  CH  =  3.  8J/cm3K,  aeq  =  a^ 
when  t  =  10_9sec.  Hence,  curve  A  applies  for  pulses  longer  than  a  nano¬ 
second  duration,  and  the  maximum  temperature  T  occurs  at  a  =  a  . 

rnax 

2 

This  example  with  a  pulse  energy  of  20  j  /  cm  and  a  pulse  duration  of  30  nsec 

8 

has  been  considered  previously  by  Hopper  and  Uhlmann.  With  this  pulse 
duration,  Tmax  occurs  at  a  =  0. 2  pm .  Using  (2.  8)  with  ( 1  -  (R> )  =  0. 2  at 
this  wavelength,  (4. 12)  yields  Tma^  =  5  X  104  K.  These  results  are  in  rea¬ 
sonable  agreement  with  results  shown  in  Hopper  and  Uhlmann' s  Fig.  3,  which 
presents  a  plot  of  T  versus  a  that  agrees  with  Fig.  4  A  only  for  the  region 
a,  <  a  <  a  =  2.3pm.  In  their  curve,  a  maximum  temperature  of  T=  2.5  x  10 

K  t 
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occurs  at  a  =  0.2|im,  the  difference  in  temperature  resulting  from  tbr 


extrapolations  made  in  curve  A.  In  the  region  a  <a^»  their  results  are  not 


valid  since  the  incorrect  cross  section  a  ^  =  c^tfa  was  used. 


In  the  case  of  dielectric  inclusions,  volume  heating  is  used  for  <  1 , 
and  surface  heating  is  used  for  /3j  a  >  1 .  The  time  T  ^  separating  the  two 


regions  of  Fig.  2  given  by  (4.  4)  and  (4.  5)  corresponds  to  a  radius  a 


,1/2 


eff 


=  (2  K„  t/Cj)  '  for  a  pulse  of  duration  t.  In  the  volume  heating  range 


with  /3j  a  <  1 ,  the  absorption  cross  section  is  proportional  to  a 6  with 


slightly  different  coefficients  for  ka  «  1  in  (2.3)  and  for  ka  »  1  in  (2.9). 


3  . 


Using  cr  ^  “  a  with  the  average  of  the  two  coefficients  for  the  entire 


volume  heating  region,  (4. 4)  gives  T£  *  a^  for  agff  <  a  <  /3j  1  and  (4.  5) 


gives  T£  «  a  for  a  <  a^  where  it  is  assumed  that  <  k .  These  results 
are  sketched  at  the  left  side  of  Fig.  5. 


For  surface  heating  of  strongly  absorbing  dielectric  inclusions  with 


j9  =  1/6  =  103  cm"1,  Cj  =  2  J/cm3  K  and  Kj  =  10_2W/cmK,  T&  is  of  the 


order  of  10  sec.  This  value  is  much  greater  than  pulse  lengths  of  interest, 
and  the  first  linear  region  of  Fig.  3  given  by  (4.  9),  which  was  only  of  academic 
interest  for  metals,  is  now  the  only  surface  heating  region  that  applies.  Be¬ 
cause  >  1  for  surface  heating  and  generally  ka  »  1,  the  absorption  cross 
section  is  =  ( 1  -  (R)  )  tra  ,  according  to  (2.  8).  Substituting  (2.  8)  for 

CTabs  *nto  ^)»  sur^ace  temperature  is  independent  of  a  .  This  is  sketched 
at  the  right  side  of  Fig.  5  for  /3j  a  >  1. 

a.  i  Fig.  5  the  maximum  temperature  occurs  for  a^  <  a  <  1  / /3j .  Using 
aabs  =  Cj  =2J/cm3K,  and  /3j  ranging  from  10  to  10^  cm  *,  the 
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failure  temperature  of  1,000K  occurs  for  pulse  energies  ranging  from  0.  3J/cm^ 


2  -4 

to  300  J  /cm  .  For  microsecond  duration  pulses  ag^  is  of  the  order  of  10  cm, 

°  -2 

and  for  nanosecond  pulses  is  of  the  order  of  300  A ,  where  =  10  W/cmK 


is  used. 


Many  other  examples,  in  nearly  every  case  of  practical  interest,  could  be 


derived  using  the  simple  results  developed  above.  The  extrapolations  from  the 


simple  limiting  cases  to  the  intermediate  regions,  such  as  using  the  dashed 


rather  than  solid  curve  in  Fig.  1 ,  tend  to  overestimate  the  temperature,  or 


underestimate  the  pulse  energy  at  the  damage  threshold,  but  only  by  factors 


typically  of  order  2  . 
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V.  CONCLUSIONS 


Simple  limiting  expressions  for  the  absorption  cross  sections  of  inclusions 
derived  in  Sec.  II  are  used  in  Sec.  Ill  to  obtain  expressions  for  the  optical  ab¬ 
sorption  coefficient  0  for  the  cases  of  large  (ka  >1)  and  small  (ka  <  1  )  dielec¬ 


tric  and  metallic  inclusions.  For  various  types  of  inclusions,  the  frequency 
dependence  of  0  ranges  from  increasing  as  U)^,  to  independent  of  to,  to  expo¬ 
nentially  decreasing  with  to.  The  temperature  dependence  ranges  from  inde¬ 
pendent  of  T  to  increasing  as  in  the  high -temperature  limit,  where  p  =;2-4 
typically.  The  examples  in  Sec.  Ill  illustrate  that  for  strongly  absorbing  dielec¬ 
tric  or  metallic  inclusions,  impurity  volume  fractions  as  small  as  f  =  10-8  can 
result  in  infrared  absorption  coefficients  of  the  order  of  10~4cm-1,  which  are 
currently  observed.  The  impurities  are  not  necessarily  limited  to  the  bulk  of 
the  crystal,  but  may  be  on  the  surface  as  would  occur  for  a  surface  contam¬ 
inated  by  polishing  compounds,  which  generally  have  large  absorption  coef¬ 
ficients. 

In  Sec.  IV  failure  due  to  local  heating  of  dielectric  and  metallic  inclusions 
in  pulsed  and  cw  systems  is  examined,  and  schematic  results  are  given  for 
many  limiting  cases.  Local  heating  is  a  iai  greater  problem  in  short-pulse 
systems  than  in  long-pulse  or  cw  systems  having  the  same  average  intensity 
as  the  short-pulse  system.  In  the  case  of  micron-size  metallic  or  dielectric 
inclusions,  pulse  energies  of  a  few  joules  per  square  centime. ter  are  sufficient 
to  cause  local  damage.  In  special  cases  where  the  inclusion  is  adjacent  to  a 
crack  or  other  imperfection  or  near  the  focal  point  of  another  inclusion  or 
other  imperfections,  the  damage  thresholds  could  be  lower  than  our  calculated 
values  by  as  much  as  two  orders  of  magnitude. 
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C.  THEORY  OF  MULTIPHONON  ABSORPTION  IN  INSULATING  CRYSTALS* 

M.  Sparks 

Xonics,  Incorporated,  Van  Nuys,  California  91406 

and 

L.  J.  Sham 

University  of  California,  San  Diego,  La  Jolla,  California  92037,  and 
Xonics,  Incorporated,  Van  Nuys,  California  91406 

The  nearly  exponential  frequency  dependence  of  the  infrared 
absorption  coefficient  0  recently  observed  in  fifteen 
crystals  up  to  several  times  the  reststrahl  frequency  is  explained  in 
terms  of  multiphonon  absorption  processes.  The  central-limit  theorem 
is  used  to  reduce  the  multiphonon  contribution  to  a  simple  closed  form. 

The  theoretical  estimates  for  the  magnitude  of  the  absorption  coefficient, 
with  no  adjustable  parameters,  are  also  in  good  agreement  with  experi¬ 
ment.  The  temperature  dependence  of  0  at  a  fixed  frequency  is  shown 

n  1 

to  be  considerably  weaker  than  0  ~  T  ,  where  n  is  the  number  of 
created  phonons.  Higher-order  processes  in  the  perturbation  expan¬ 
sion  are  shown  to  be  negligible  for  small  n  ,  to  be  comparable  to  that 
of  the  lowest-order,  single-vertex  terms  for  n^5,  and  to  dominate 
for  large  n  in  a  typical  case.  Difference  processes,  in  which  some 
thermally  excited  phonons  are  annihilated,  are  shown  to  be  negligible 
with  respect  to  the  summation  processes  in  the  nearly  exponential 
region.  An  explanation  involving  finite  phonon  lifetimes  is  proposed 
to  explain  the  fact  that  the  alkali  halides  show  less  structure  in  the 
|3-oo  curves  than  do  the  semiconductor  crystals. 
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I.  INTRODUCTION 

The  intensity  I  of  infrared  radiation  propagating  through  a  solid  typically 
decays  according  to  Beer’s  law,  I  =  IQ  exp(-0z),  where  0  is  defined  as  the 
optical  absorption  coefficient.  Extensive  experimental  and  theoretical  studies 
have  been  conducted  on  the  absorption  due  to  phonons  in  insulating  or  semi¬ 
conducting  crystals.  Refs.  1-4  represent  some  recent  reviews  on  this  topic. 

The  main  interest  has  been  focused  on  the  two -phonon  region  where  0»  1  cm”1, 
and  particularly  on  the  structure  of  the  frequency  dependence  that  determines 
the  critical  points  of  the  phonon  spectra.  5  The  availability  of  high-power 
infrared  lasers  has  shifted  attention  to  higher-order  phonon  processes,  where 
0  «  1  cm  l.  Not  only  the  positions  of  the  multiphonon  peaks  are  of  interest, 
but  also  the  magnitude  of  0  is  of  great  importance  now  that  high  intensities 
are  available. 

6  7b 

It  has  been  observe  1  ’  '  that  for  frequencies  to  greater  than  several 

times  the  reststrahl  frequency  tOp  the  optical  absorption  coefficient  varies 
nearly  exponentially  vdth  frequency, 

0  ~  exp  (-Aw)  ,  (1.  l) 

for  a  number  of  crystals  including  LiF,  NaF,  NaCl,  KC1,  KBr,  MgF2 ,  CaF2, 
Bal'2,  Srh2,  MgO,  A^O^,  Si02 ,  TiC>2  ,  BaTiO^ ,  and  SrTiO^ .  This  is  true  for 
0  £  10  cm  1  and  to  ^  2  to^- ,  roughly.  In  NaCl  at  room  temperature,  for  instance 
0  decreases  nearly  exponentially  for  over  four  orders  of  magnitude  as  the  fre¬ 
quency  increases  from  2.2  to f  to  5.  8  oof  ,  as  shown  in  Fig.  1 . 
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Figure  1.  Experimental  frequency  dependence  of  the  infrared  absorption  coefficient  j3  for 
N'aCl  after  Horrigan  and  Deutsch  (+,©  )  Ref.  7,  Smart, et  al  (— •  — )  Ref.  7a,  and  Genzel  (— ) 
Ref.  7b . 
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At  first  sight,  the  nearly  exponential  behavior  might  suggest  the  form 

0  „  exp  (-■Rw/kjjT).  However,  the  room -temperature  values  of  the  coeffi¬ 
cient  A  in  (1.1 1  differ  by  factors  of  2  -  4  from  the  value  of  ft/ kg T.  Further¬ 
more,  the  temperature  dependence  8  of  0 ,  though  not  extensively  studied  to  date, 
appears  to  be  less  strong  than  exp(-hto/kBT). 

In  this  paper,  an  investigation  of  the  optical  absorption  by  multiphonon 
processes  is  presented.  It  is  shown  that  the  sum  of  n-phonon  summation  pro¬ 
cesses  is  approximately  exponentially  decreasing  with  increasing  frequency 
over  the  frequency  range  of  interest,  i.  e.,  about  2  tOf  -  7  tof  ,  typically.  As 
illustrated  in  Fig.  2  a  ,  we  consider  the  n-phonon  summation  process  in  which 
the  photon  is  absorbed  by  the  crystal  through  the  virtual  excitation  of  the  funda¬ 
mental  reststrahlen  mode  which  finally  emits  n  phonons.  In  other  words,  the 
electromagnetic  field  drives  the  fundamental  mode  (off  resonance  since  W  >  Wf  ) , 
whose  relaxation  time  is  determined  by  the  sum  of  all  possible  processes  of 
splitting  into  n  normal  modes  of  lattice  vibrations.  The  Lax-Burstein-Born 
higher-order  dipole-moment  mechanism9  is  not  considered  explicitly,  although 
most  of  the  analysis  still  applies  to  that  case. 

By  energy  conservation,  the  energy  ftw  of  the  photon  absorbed  is  equal  to 
the  sum  of  the  energies  of  the  n  final-state  phonons.  It  follows  that  the  n-phonon 
summation  process  cannot  contribute  to  0  when  to  >  n  cogr ,  where  is  the 
greatest  frequency  of  the  phonon  spectrum.  For  to  «  n  cogr  ,  the  contribution 
0  of  the  n-phonon  summation  process  to  0  is  small  because  the  low  frequen¬ 
cies  of  the  final- state  phonons  greatly  restrict  the  amount  of  phase  space  avail¬ 
able.  Thus,  0n  must  peak  at  a  frequency  not  far  below  nto^ .  As  n  increases, 
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the  peak  shifts  to  higher  frequencies  and  decreases  in  height  since  higher-order 
phonon  processes  involve  weaker  coupling  coefficients.  The  sum  of  the  then 
has  a  frequency  dependence  nearly  exponential  in  the  experimental  frequency 
range.  This  behavior  of  the  0n  and  the  sum  of  the  0n  is  demonstrated  expli¬ 
citly  in  Sec.  VI . 

A  preliminary  account  of  these  results  lias  been  published. 10  Subsequent 
investigations  are  discussed  in  Sec.  II .  The  exponential  frequency  dependence 
of  the  absorption  was  fiist  suggested  by  Rupprecht6  to  be  due  to  n-phonon 
processes,  although  he  did  not  investigate  the  theory  in  detail. 


In  Sec.  II,  formal  expressions  for  the  contribution  to  0  due  to  multiphonon 
processes  are  given.  A  practical  approximation  for  the  anharmonic  coelficient 
is  chosen.  In  Sec.  Ill ,  an  asymptotic  approximation  for  evaluation  of  the 
n-phonon  contribution  is  developed.  In  Sec.  IV  ,  confluence  phonon  processes 
are  shown  to  be  unimportant  in  the  nearly  exponential  region.  In  Sec.  V,  all 
possible  processes  that  convert  the  fundamental  phonon  to  n  phonons  are  ex¬ 
amined,  and  the  contributions  of  vertex  corrections  are  estimated.  In  Sec.  VI, 
the  explicit  evaluation  of  0n  is  described,  and  comparison  of  theory  with  experi¬ 
ment  is  made.  In  Sec.  VII ,  a  summary  of  all  the  assumpations  and  approximations 
that  have  gone  into  the  theory  is  given,  and  the  relatnn  of  a  computer-calculation 
program  to  the  present  results  is  discussed. 
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11.  ANHARMONIC  CONTRIBUTION  TO  THE  ABSORPTION  COEFFICIENT 


The  infrared  radiation  perturbs  the  insulating  crystal  by  excitation  of  the 
dipole  moment  of  the  crystal  by  the  oscillating  electric  field.  The  absorption 
coefficient  is  simply  related  to  the  imaginary  part  of  the  electric  susceptibility 


j3(w)  =  4ttx1(^)  to/nr  c  ,  (2.1) 

where  c  is  the  speed  of  light  and  nr  is  the  refractive  index  at  frequency  a,'. 

3  11 

Tlie  susceptibility,  in  turn,  is  just  the  linear  response  of  the  dipole  moment. 

In  an  anharmonic  crystal,  the  dipole  moment  can  be  expanded  in  powers  of  the 

12 

ionic  displacement  s.  For  infrared-active  crystals,  the  leading  nonzero 
term  is  linear  in  the  ionic  displacements.  The  nonlinear  terms  (the  dominant 
mechanism  for  infrared  absorption  in  such  non-infrared-active  crystals  as 


9  13  14 

diamond  )  are  probably  small  in  polar  crystals,  especially  in  alkali  halides,  ’ 

and  shall  be  neglected  in  this  work.  However,  there  are  contrary  conclusions.1'’ 

Then,  the  absorption  coefficient  is  given  by  the  imaginary  part  of  the  Green's 

3  13 

function  of  the  fundamental  mode  ’ 


.  4  it  Me*2  O  w  Mf  r(M> 

cmrnr  a  <w2-w2)2  +  [wfr(u))]2 


(2.2) 


where  N  is  the  number  of  unit  cells,  0  the  volume  of  the  crystal,  e  is  the  Bora 

effective  charge,  mr  the  reduced  mass  of  the  two  ions  in  the  unit  cell,  and  T  is 

the  energy  relaxation  frequency  of  the  fundamental  mode  (equal  to  twice  the  T  in 

13 

R.  A.  Cowley's  notation  ).  The  real  part  of  the  phonon  self-energy  is  understood  to 
have  been  included  in  producing  the  renormalized  reststrahl  frequency  to^  ,  and  its 
frequency  dependence  is  neglected  in  Eq.  (2.2).  A  simple  classical  model  of  a 
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harmonic  oscillator  (the  fundamental  lattice  mode)  driven  by  the  applied  electric 

12 

field  gives  (2.2),  but  with  cofr  replaced  by  cor  in  the  numerator  and  denominator. 


The  contribution  Tn  from  the  n-phonon  summation  processes  to  T  can  be  calcu- 
16 


lated  by  applying1  °  the  standard  perturbation- theory  result  that  the  probability  per 

unit  time  of  a  transition  between  two  states  is  2  it  /la  times  the  product  of  the 

square  of  the  matrix  element  and  the  energy  conserving  delta  function,  giving 

o  9  on  n 

2tt  z 


r  (w)=  2j(n+l)  n!  E  |A(fQ1*-*Qn)|  A(  E  n  . )  6(co  -  £  coQ  )  n  , 
n  f/  1  ^  j-i'-J  j  =  l  n 


(2.3) 


where  Q.  is  the  phonon  mode  with  wavevector  q  .  and  branch  b. ,  A  is  the  modified 
J  J  J 

Kroriecker  delta  which  is  unity  when  the  argument  is  ze  ro  or  a  reciprocal-lattice 
vector  and  zero  otherwise,  and 


n 


nn  =  n.  ("j+D/n  +1) 

J=1  J 


(2.4) 


with 


Wq/  Wt 

n.  =  n(Q.  )  =  l/[e  ^  -1]  , 


J  '  T 

OO/ClJ'-p 

=  l/(e  1  -1)  , 


(2.5) 

(2.6) 


and 


C0T  =  kgT/'fi  . 


(2.7) 


Furthermore,  A(fQj*  •  •  Q^)  denotes  the  renormalized  n  +  1  phonon  vertex, 
represented  by  the  circle  in  Fig.  2  a ,  and  is  the  sum  of  all  possible  n+1  phonon 


vertices.  The  simplest  one  is  the  unrenormalized  vertex  V  (fQj  '  "  Qn)  from 

12, 13 


the  anharmonic  Hamiltonian  given  by 


n+1 


K 


n+1 


q  ‘q  ' 

^1  vnfl  J 


(2.8) 
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where 


A,^  =  /  ,  =  a  ,  +  a  ,  , 

Q  qb  qb  -qb 


(2.9) 


a  and  a  being  the  phonon  creation  and  annihilation  operators  normalized  to 
unit  commutators,  as  usual.  This  simple  vertex  is  represented  diagrammatically 
in  Fig.  2b.  Other  more  complicated  processes  are  examined  in  Sec.  V, 
where  we  derive  an  approximate  form  for  the  total  vertex 


A(fQ1-**Qn)  =  An  V(fQ1--*Qn) 


(2. 10) 


To  obtain  a  reasonable  approximation  for  the  anharmonic  coefficients,  let  us 

confine  our  attention  to  diatomic  polar  crystals  with  cubic  symmetry,  especially 

12  17 

NaCl- structure  crystals.  The  model  interaction  potential  ’  between  ions  is 
composed  of  a  Coulomb  potential  and  a  nearest-neighbor  overlap  exchange  re¬ 


pulsion  of  the  form 


0(r )  =  C  exp  (  -  r/pa  )  , 


(2. 11) 


where  a  is  the  equilibrium  nearest-neighbor  distance.  The  Coulomb  interaction 

is  used  only  in  determining  the  constants  C  and  p  in  Eq.  (2.  11)  from  the  equili- 

12  17 

brium  condition  and  the  value  of  the  bulk  modulus  ,  B  ,  yielding  ’ 


C  =  3  a3  B  e1/p  p2/(l  -  2 p)  . 


(2. 12) 


In  the  anharmonic  coefficients,  only  the  derivatives  of  the  repulsive  potential  (2. 11) 
are  retained.  Since  p  is  of  the  order  of  0. 1  for  NaCl ,  the  derivatives  of  the 
Coulomb  potential  are  smaller  than  the  corresponding  ones  of  the  repulsive 
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potential  for  orders  up  to  at  least  n  —  10.  Had  we  used  an  inverse  power  law 
for  the  repulsive  potential,  this  would  be  true  to  any  order.  This  model,  in¬ 
cluding  the  neglect  of  the  Coulomb  potential  in  the  anharmonic  terms ,  has  been 


13 


used  previously  with  much  success. 


The  anharmonic  coefficients  can  be  obtained  in  a  straightforward  calcula¬ 
tion  from  this  nearest-neighbor  exchange  repulsion  potential. 12  The  mth  order 
coefficient  V(Qj  •  ■  *Qm)  involves  derivatives  of  0(r)  up  to  order  m  .  From 
the  exponential  form  (2. 11),  it  is  clear  that 


|a0<m>(a)/0<ra-1>(a)|~lO  . 


(2.13) 


Thus,  it  is  a  good  approximation  to  retain  only  the  highest-order  derivative. 
Using  these  results  and  assuming  central  forces  yields 


6  m 


/  v  U  UI  1/2 

v (Qf  *  * Qm )  =  ( N/m  1)0  (a)  £  n  U  (Q  )(D/2Nm<w0  )  ,  (2.14) 

y=  l  j=i  7  1  vj 


where 


r  1/2  lq  •  x 

uy(<y  =  xy  •  [Scq -<m</ra>)  X>Q.e~  ~y 


]  ' 


(2. 15) 


and  m^j.  and  denote  the  smaller  and  larger  ionic  masses,  respectively. 


The  positions  of  nearest  neighbors  measured  from  the  lighter  ion  are  x^  ,  and 


x  is  the  unit  vector  in  the  same  direction.  The  polarization  vector 


w 


is 


"y  “  - - - - - -  *  ••  WWWAV..I  •  *»W  1UU11UII  »  CULUl  VV  Q 

defined  in  terms  of  the  ionic  displacement  u^  from  the  equilibrium  position 


x  .  bv  the  relation  ^2’  ^ 
T 
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1  I 


1/2  >3  •  J,T 

u,T  =  E(ti/2NmTup)  e  AQ  wtQ  , 


(2.  16) 


with  r  denoting  the  ion  type.  For  the  fundamental  mode, 


Uy(f)  =  xy  •  wf  (m</mr) 


(2. 17) 


I  |  *  th 

with  mf  =  (m<”J  +  m>  )  .  From  Eq.  (2.11),  we  obtain  the  rn  derivative 


0(m)(a)  =  3Ba3p2/  j  (l-2p)(-pa)m 


(2.18) 


Substituting  the  approximate  expression  (2. 14)  for  the  anhanaonic  coefficient 
into  Eq.  (2.3)  for  r  ,  we  obtain,  by  using  Eq.  (2.2),  the  contribution  of  the 
n -phonon  summation  process  to  the  optical  absorption  in  the  form 

1/2  5  -to/coT  4-1  n  2 

j3  =  (  it  /  2 )  K  to  3  D  1-e  1  )  (to  n! )  (to  D  )  Az£  . 

fn'//  mx  p  ' '  'mxe  nn 


(2. 19) 


2  2  2 

We  have  used  the  approximation  for  high  frequency  (to  »  tof  +  T  )  and  intro¬ 


duced  the  following  groups  of  constants: 


K  =  B2  e  a  to, /lie  m  2  n  to  3 
f  r  r  mx 


i 


Dp  =  ( 2  it  )  [6Trp/(l-2p)]  , 


(2.20) 


2  2 

D  =  Ti  /  2  p  a  nv  to 
e  H  <  mx 


The  frequency  to  ,  introduced  for  later  use,  cancels  out  in  Eq.  (2.  19). 


(2.21) 


1  y&  c 

For  crystal-  of  NaCl  structure,  symmetry  ensures  that  Zn  and,  therefore, 
0n  are  independent  of  the  direction  of  wf  .  Let  us  choose  wf  to  be  along  the 
positive  x-axis  .  Then,  Eq.  (2.21)  becomes 


Sn  =  2Sn++2(’1)  Sn- 


(2.22) 


where 


2  =  N  "n  Y  NA(Lq.)6(to-E  Wq  )  n  W+(Q  )(nQ  +1)/o^q 

n±  Q  — Q  j=l~J  j=l  j=l  "  J  J  J 

i  ^  /n  no 


and,  with  Re  and  Im  denoting  real  and  imaginary  parts,  respectively, 

2  2 

W  (Q)  =  [Re  Ux(Q)]  ±[ImUx(Q)]  . 


(2.23) 


(2.24) 


In  passing,  notice  that  the  evaluation  of  the  sums  in  (2.23)  is  trivial  if  the 
density  of  states  g(to)  is  approximated  by  the  Einstein  model 


g(to)  =  (to  -oje  ) 

and  the  angle  dependence  of  W± ( Q  )  is  neglected  :  W±(Q)  =  W±  .  It  will  be 
shown  later  than  W+n  »  W_n.  Then  (2.23)  and  (2. 19)  give  directly 


6  6 
=  E  Z 
y— 1  y'=l 


The  factor  Z  contains  the  dynamical  information  of  the  n -phonon 

fi- 

ab  sorption,  and  is  given  by 
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p  _  |H2  |  v^w1-1}  8(B-n 

n=l  n4n!  '  'I  WE  ) 


aJE)  ’ 


(2.25) 


where  E  =  (2 it )  /  Kcc^D^.  According  to  (2.25),  the  spectrum  is  approximated 
by  a  series  of  delta  functions,  which  is,  of  course,  not  realistic.  Even  though  such 
a  model  is  not  of  significant  practical  value,  it  does  crudely  approximate  some  of 
the  features  of  the  more  realistic  model  discussed  below.  For  example,  plotting 
the  coefficients  of  the  delta  functions  in  (2.25),  or  formally  replacing  the  delta 
functions  by  line-shape  functions  of  finite  width,  gives  a  nearly  exponential  de¬ 
crease  with  increasing  frequency.  In  Ref.  19b,  the  result  (2.25)  was  rederived 
using  the  simpler  model  of  a  one -dimensional  lattice  with  the  Einstein  approxima¬ 
tion  and  a  simpler  interaction  potential  that  neglects  the  angle  dependence  [our 

factor  W+(Q)]  from  the  outset, and  an  independent -molecule  model  was  considered. 

2 

Use  of  this  simpler  interaction  potential  gives  unreasonably  large  values  of  An  , 

19b 

which  causes  noticeable  deviation  from  an  exponential  frequency  dependence. 

19c 

Mills  and  Maradudin  independently  used  a  single -frequency  anharmonic- 

molecule  type  lattice  to  study  various  types  of  interaction  potentials,  effects  of  im- 

19d 

purities,  and  high -temperature  effects.  Bendow,  Ying,  and  Yukon  have  used  a 
differ  jnt  mathematical  method  that  starts  with  partially  summed  terms.  The  method 
is  potentially  powerful,  but  to  date  they  have  recovered  only  our  terms  without  the 
vertex  correction.  Since  the  validity  of  the  perturbation  expansion  is  justified  by 
showing  that  all  diagrams  not  included  in  the  result  are  negligible,  it  is  expected 
that  new  methods  of  calculation  should  give  equivalent  results.  Foj  example,  the 
factor  exp  (-y  T )  resulting  from  vertices  with  phonon  loops  [  AqAqS  ,  AqAqAq,Aq,S  , 
etc. ,  where  S  is  the  simple  vertex  and  the  Aq  are  defined  in  (2.9)]  is  well  approxi¬ 
mated  by  1  since  the  phonon-loop  terms  are  negligible  with  respect  to  simple  ver¬ 
tices.  See  Sec.  V . 
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HI.  ASYMPTOTIC  APPROXIMATION’  FOR  ABSORPTION  BY  A 
LARGE  NUMBER  OF  PHONONS 


Eq.  (2.  19)  gives  the  contribution  to  the  absorption  coefficient  by  the 
n-phonon  summation  process.  It  contains  the  factor  given  by  Eqs.  (2.22) 
to  (2.24)  which  involves  n-fold  Brillouin -zone  sums.  Although  these  are  not 
beyond  the  means  of  modem  computing  capabilities  for  n  in  the  experimental 
range  of  2  to  8 ,  we  are  still  interested  in  analytical  approximations  that  will 
give  us  general  properties  of  )3  which  appear  to  be  shared  by  a  rather  large 
number  of  crystals.  The  method  of  evaluation  used  in  this  section  is  correct 
in  the  large-n  limit. 

For  n^  2  we  can  neglect  the  quasi-momentum  conservation  restriction  in 

the  sums  given  by  Eq.  (2.23).  We  shall  justify  this  laer  in  the  section.  First 

notice  that  if  the  angle  dependence  of  W+((J.)  is  neglected,  the  summand  in  (2.23) 

is  a  function  of  phonon  frequencies  gOq  only.  Then,  replacing  the  sums  over 

j 

Q i  by  integrals  over  dten  g(tOn  )»  where  g(con  )  is  the  phonon  density  of 
J  Llj  ^ 

states,  reduces  (2.23)  to  the  form 


n± 


Jd“Qlf<“Pl)'"  /d“Qnf<"Q„ 


)  6  ( a)  -  L  ccn  ) 

j=l 


to  which  the  central-limit  theorem  applies  directly.  Here  f(u.'Q  ) 

.1  i 

=  6N  W+(Qj  )(nq  + 1 )  g(cOq  )/Wq  ,  where  the  normalization 


constant 


_1  j  j  j 

6  N  arises  since  g  is  normalized  to  unity. 


Eq.  (2.23)  can  be  cast  into  this  central-limit-theorem  form  without  neg¬ 
lecting  the  angle  dependence  of  W+(Qj)  as  follows:  We  introduce  two  functions 
which  are  kindred  to  the  phonon  propagators 
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0,(0  =  l^/K^tj  SW,(Q>  j<V1)/“Qj  £,<C-“q)  •  (3'" 


where 


a0±'“n«N''‘^w±<Q)<nQ  +  1,/“Q 


(3.2) 


are  constants  for  normalizing  the  integrals  of  o+(C)  over  £  to  unity. 

The  n-fold  sums  in  Eq.  (2.23)  can  be  written  as  n-dimensional  integrals, 


00  00 

S„t  =  <°0±/a;mx  "  j  d«10±«:i>-"  _[  <iCnO,(Cn)  eiw-L^j).  (3.3) 


These  convolution  integrals  are  well  known  in  statistics.  For  n  -»  ®  ,  the  integral 
tends  to  a  Gaussian  (the  central-limit  theorem), 


Sn±  =  [ao±n/(27rn)1/2o£2±Wmxn+1]  exp  [  "  (t0  ’  nai±  Wmx)  /2na2±Wmx  ] 

(3.4) 


where  aQ+  are  defined  in  Eq.  (3.2),  and 


al±  =  Winx 


-1 


00 

\  d^± 


<oe  . 


(3.5) 


(a 


2± 


V»2=  \  dt°±<?)  ?2  '  (“1±wmx  >Z  ' 


(3.6) 


For  small  n ,  it  is  possible  to  improve  Eq.  (3. 4)  with  an  asymptotic  series, 


20 


A  particular  series  in  terms  of  Hermite  polynomials  has  been  used  by  Sjolander 
to  evaluate  the  multiphonon  background  in  neutron  scatterings  in  a  harmonic 
crystal. 
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It  is  obvious  from  Eq,  (2.24)  that  G!q_  ,  as  defined  by  Eq.  (3.2),  is  less 
than  .  The  estimates  discussed  in  Sec.  VI  show  that  0!q_  is  about  one 
third  to  one  half  of  aQ+,  at  most.  Since  SQ±  ~  (aQ±  f  according  to  (3.  4),  Ln_ 
becomes  negligible  compared  with  for  large  n.  Therefore,  from  Eqs.  (3.4), 
(2.22),  and  (2. 19),  the  absorption  coefficient  has  the  explicit  form 

4 


D„  K  /  cu  \  "  ,  nn 

ferr  (-?)  3^ 


to 


^n  a 


2+ 


exp  /^'Vm'  ] 


(3.7) 


By  virtue  of  the  central-limit  theorem,  the  multiple  sum  over  Qj  ,  •  •  •  Q  has 
been  reduced  to  sums  over  a  single  phonon  coordinate  Qj  ,  as  given  by  Eqs.  (3.2), 
(3.5),  and  (3. 6). 


The  neglect  of  momentum  conservation  appears  to  be  physically  reasonable 
since,  for  larger  and  larger  n  ,  the  restriction  cn  phase  space  becomes  less 
and  less  important.  However,  if  we  wish  not  to  neglect  the  momentum  conserva¬ 
tion  in  Eq.  (2.2),  we  can  extend  the  foregoing  procedure  by  treating  the  summa¬ 
tions  over  q  in  the  same  manner  as  the  integrals  over  £. .  Thus,  we  introduce 
"*J  J 

the  functions  a,  (q,  £ )  similar  to  Eq.  (3. 1) ,  except  omitting  the  sum  over  q. 

'  r*j  ** 


Eq.  (3.  3)  becomes  not  only  multiple  integrals  over  £.  but  also  over  q  .  with 

J 


four  6  functions,  one  for  the  frequency  and  three  for  the  wavevectors.  The 
convolution  integral  is  evaluated  in  the  same  way  by  means  of  the  central-limit 


theorem.  The  integrals  over  q.  contribute  a  factor  which  is  a  lattice  sum  of 

2 


Gaussians  of  the  form  exp(-na  x  ^  )  and  is,  therefore,  approximately  unity 
for  large  n.  We  arrive  at  the  same  answer  as  Eq.  (3.  4),  thereby  justifying 
the  neglect  of  momentum  conservation. 
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IV.  THE  CONFLUENCE  PROCESSES 


In  the  preceding  calculation  of  the  multiphonon  absorption  of  light,  only  a 
particular  type  of  phonon  processes,  called  the  n-phonon  summation  processes 
and  Illustrated  in  Fig.  2a  ,  was  considered.  We  have  neglected  the  confluence 
processes  ,  illustrated  in  Fig.  2  c  .  Instead  of  creating  n  phonons  after  the 
annihilation  of  the  fundamental  phonon,  m  phonons  are  absorbed,  and  n-m 
phonons  are  created. 

A  confluence  process  involving  n  phonons  (some  of  which  are  created  and 
some  annihilated)  is  governed  by  the  same  vertex  as  the  n-phonon  summation 
process.  The  contribution  to  the  absorption  coefficient  of  all  confluence  pro¬ 
cesses  and  the  summation  process  is  easily  obtained  by  replacing  6  (£  -  Wq) 
by  6(C  -Wq)  -  6(C  +  Wq)  and  n(-u)Q)by  -[n(uiQ)+  1]  in  (3.1).  The  cross- 
product  terms  in  (3.3)  containing  m  factors  of  -6(C  +  u>Q)  and  n-m  factors 
of  6  (C  -  0Oq  )  correspond  to  the  confluence  process  with  m  thermally  excited 
phonons  absorbed,  as  shown  in  Fig.  2  c  .  Applying  the  central-limit  theorem 
to  this  term  yields  a  Gaussian  peaked  at  (n-m)ai+o)mx-m0j+Wmx  instead  of 
n  a  oo  ,  where  fl, ,  is  defined  by  Eq.  (3. 5)  with  the  new  a (C  )  •  ™s 
contribution  will  be  masked  by  the  summation  process  of  n-2  m  phonons  which 
peaks  at  about  the  same  frequency  but  has  greater  strength,  being  an  anharmonic 

process  of  lower  order. 
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V.  VERTEX  CORRECTIONS 


th 


Now  we  consider  all  possible  processes  that  contribute  to  the  (n+1)  -order 

vertex  A(fQ,  •••())  and  estimate  the  vertex  correction  factor  A  ,  defined 
in  n 

by  Eq.  (2.  10).  Standard  perturbation  theory  can  be  applied  in  a  straight¬ 
forward  way  to  all  the  higher -order  terms.  For  example,  for  n  =  3,  the 
diagram  in  Fig.  4  b  below  has  one  intermediate  state.  The  contribution  from 
this  diagram  is  easily  calculated,  but  it  must  be  remembered  that  this  dia¬ 
gram  represents  four  diagrams  when  the  arrow  is  added  to  the  intermediate - 
state  phonon.  (There  are  two  time  orderings  of  the  two  vertices,  and  the 
arrow  can  go  in  either  direction,  corresponding  to  a  a^  and  a^  a ,  in  each 


time  ordering. )  This  procedure  has  been  carried  out  for  a  number  of  low- 


10 


order  diagrams,  and  the  results  agree  with  those  presented  below. 


Since  the  number  of  diagrams  increases  rapidly  as  n  increases,  this 
method  becomes  tedious  and  time  consuming  when  applied  to  larger  values 
of  n  .  The  following  method  is  more  convenient.  First,  all  the  self-energy 
corrections,  such  as  those  illustrated  in  Fig.  3a  ,  are  taken  to  be  accounted 
for  by  using  the  measured  phonon  frequencies,  i.  e. ,  they  are  included  in  the 
corresponding  "skeleton"  diagram  (Fig.  3  b).  The  lifetime  of  the  intermediate- 
and  final-state  phonons  is  taken  to  be  infinite. 


There  are  two  types  of  vertices:  (1)  the  irreducible  ones  that  cannot  be 
rent  asunder  by  cutting  a  single  phonon  line,  such  as  those  in  Fig.  3c ,  and 
(2)  the  reducible  ones  that  can  be  separated  by  cutting  a  single  line,  such  as 
those  in  Fig.  3d. 
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In  the  sum  of  all  irreducible  vertices  of  the  same  number  of  external 


22  . 


phonon  lines,  the  simple  vertex  dominates.  We  follow  Van  Hove  in  order- 

n-2 


ing  the  anharmonic  terms  in  the  Hamiltonian,  V(Qj  •  •  *Q  ),  with  e 
where  e  is  the  small  parameter  given  by  the  ratio  of  the  root -mean -square 
displacement  of  the  ions  to  the  nearest -neighbor  distance.  The  value  of  e 
is  less  than  0.05  in  alkali  halides.  A  complex  irreducible  vertex  must  be  of 
higher  order  in  e  than  the  simple  vertex  with  the  same  number  of  external 
lines,  since  cutting  a  phonon  line  will  produce  one  vertex  with  a  larger 


number  of  external  lines.  For  example,  the  simple  vertex  in  Fig.  3  c  is 


2  4 

0  (  e  ) ,  but  all  the  other  irreducible  vertices  in  Fig.  3  c  are  0(e). 


On  the  other  hand,  a  reducible  vertex  composed  of  simple  irreducible  ver¬ 
tices  is  of  the  same  order  in  c  as  the  simple  vertex  with  the  same  number  of 


external  lines.  For  example,  the  first  diagram  in  Fig.  3d  is  0(e  ).  A  reducible 
vertex  that  contains  one  or  more  complex  irreducible  vertices  is  again  negli¬ 
gible.  Therefore,  for  the  total  vertex  contribution,  we  need  only  sum  the  simple 
vertex  and  the  reducible  vertices  which  are  composed  of  simple  vertices  only. 


To  illustrate  the  procedure  of  obtaining  the  vertex  renormalization  to  the 
n-phonon  summation  process,  the  simplest  non-trivial  vertex  correction,  namely 


A„,  is  first  calculated.  The  two  vertex  terms  that  contribute  to  the  three-phonon 
o 


absorption  are  given  in  Fig.  4  a  and  4  b .  The  ratio  of  the  latter  to  the  simple 
vertex  is 


ilji  £  V(fQ3Q4)  D(Q4,C4)  V(Q4Q1Q2)/|f  V(fQ1Q2Q3)  .  (5.1) 

*  Q4  ’ 


62 


Iinr  II  ilMifUMjaimi  1'iiri  ~  (i 


a  I  I  tr  fk'ii v .  u 


Sec.  C 


23  13 

The  two  factors  of  3!  represent  the  number  of  ways  ’  the  phonon  states  aie 
attached  o  the  limbs  of  each  vertex  in  Fig.  4b.  The  divisor  2  !  represents  the 
fact  that  interchanging  the  labels  on  phonon  lines  1  and  2  again  produces  the  same 
term.  The  factor  4!  is  the  number  of  ways  the  four-phonon  vertex  in  Fig.  4  a 
can  be  labelled,  and  the  divisor  3 !  is  the  overcounting  factor  generated  by  re¬ 
arranging  the  labels  among  lines  1,  2  and  3  in  Fig.  4a.  Tie  factor  D(Q^,  C^) 

23  13 

represents  the  Green's  function  for  the  intermediate  phonon  line“  ’  in  Fig.  4  j  , 

By  using  the  form  (2.  14)  for  the  anharmonic  coefficient  and  keeping  onh  one 
term  in  the  sum  over  nearest  neighbors  for  both  the  numerator  and  the  denomi¬ 
nator  of  the  ratio  (5.  1) ,  we  obtain  the  ratio  as 


0<2>  2 
2  '  , 
b4 


|ux(Q4)|2D(Q4>?4) 


<5. 2) 


The  factorial  that  represents  the  number  of  ways  the  states  in  each  vertex  are 
labelled  cancels  neatly  the  factorial  in  the  anharmonic  coefficient  (2.14),  leaving 
the  counting  factor  in  front  of  (5,2).  This  factor  is  just  the  ratio  of  the  number 
of  ways  of  rearranging  the  labels  on  the  equivalent  outgoing  phonon  lines  of 
Fig.  4  a  to  the  corresponding  number  for  Fig.  4  b . 

Tie  factor  0^  comes  from  the  fact  that 


!#(3)|7*<4).  /2> 


(5.3) 


by  virtue  of  Eq.  (2. 18).  The  momentum  and  frequency  of  the  intermediate-phonon 
Green's  function  are  given  by 


34  =  3.1  +  32 


(5.  4) 
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1  hen,  the  ratio  of  the  two  vertices,  (5.2),  becomes  3d2  £,  and  the  vertex 
renormalization  factor  is 

A3  =  1  +  3^  £  .  (5-9) 

Estimates  of  d2  and  £  are  provided  in  the  next  section.  It  is  easy  to  verify 
that  (5.  9)  correctly  accounts  for  all  the  three-phonon  absorption  processes 
shown  in  Fig.  4  c. 

The  reasoning  used  in  this  simple  example  can  be  applied  to  the  general  case 
to  deduce  the  rules  for  writing  down  the  renormalization  factor  for  n -phonon  ab¬ 
sorption.  The  simplest  term  in  the  total  vertex  is  the  simple  n  +  1  phonon  vertex 
with  one  of  the  phonons  being  the  fundamental  mode  driven  at  the  optical  frequency, 
as  depicted  by  Fig.  5a.  A  typical  reducible  vertex  is  formed  by  joining  a  number 
of  4 r reducible  vertices  of  lower  order  such  that  there  is  only  one  phonon  line  con¬ 
necting  any  pair  of  irreducible  vertices.  Some  examples  are  shown  in  Fig.  5b  -5 d. 
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If  a  vertex  contains  m  internal  lines,  then  its  ratio  to  the  simple  vertex 
(Fig.  5  a)  contains  a  factor  £  m,  with  £  defined  by  Eq.  (5.  8).  Thus, 


l  =  Vs(ra)(m  , 

n  ^  n  s  ’ 


(5. 10) 


with  the  coefficient  obtained  as  follows.  Draw  all  topologically  distinct 

reducible  vertices  with  m  internal  lines  and  n+1  external  lines,  one  of  which 
is  the  fundamental  phonon  driven  at  frequency  n  U3  .  Each  diagram  contri¬ 
butes  to  a  term  of  the  form 


_  (m) 

c  d.  d4 
n  /j  i2 


(5. 11) 


(jjj) 

where  is  the  ratio  of  n  1  to  the  number  of  ways  of  rearranging  the  states 

of  the  n  outgoing  phonon  lines  that  do  not  change  the  reducible  vertex.  The 
factors  of  d^  come  from  the  intermediate  phonon  lines,  SL  being  determined  by 
energy  conservation,  assuming  that  all  outgoing  phonon  lines  have  frequency  to 

ITU 

For  example,  the  vertices  with  one  internal  line,  as  in  Fig.  5b,  give 


s (1)  =  "i:  ( n  \  d  , 

n  m/  m 

m=2  '  ' 


(5. 12) 


and  Fig.  5  e  contributes  to  S 


(n-2) 


the  term 


( n  I  /  2  1 )  d2d3  •••  dn-1 


(5. 13) 
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Armed  with  the  general  rules,  we  can  calculate  the  contribution  of  any  vertex. 


Figs.  6-8  show  the  relevant  vertices  for  four-  to  six-phonon  absorption,  respec¬ 
tively,  and  the  corresponding  contributions  to  S 


From  the  considerations  in  the  next  section,  the  factor  v  in  Eq.  (5.  6)  is 


0.5  or  less;  thus  (5.7)  gives 


Sec.  C 


Figure  6.  The  4-phonon  absorption  vertices.  The  only  external  line  shown 


is  the  fundamental  mode. 
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Figure  8.  The  6-phonon  absorption  vertices. 
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VI.  FREQUENCY  AND  TEMPERATURE  DEPENDENCE 
OF  THE  ABSORPTION  COEFFICIENT 

It  is  straightforward  to  evaluate  0^  given  by  Eq.  (3.7)  with  the  Ot 's  given 
by  the  Brillouin-zone  sums.  We  shall  confine  ourselves  to  two  rough  estimates 
of  the  ct's . 

For  a  linear-chain  model  with  two  atoms  per  unit  cell,  all  with  equal  masses, 

it  is  possible  to  evaluate  explicitly  the  a's  in  the  low-  and  high -temperature 

limits.  Table  I  shows  the  results  for  the  a's  with  us  chosen  to  be  the  top 

mx 

of  the  phonon  spectrum.  This  simple  model  illustrates  nicely  all  the  important 

features  that  follow  from  Eq.  (3.7).  As  a  function  of  frequency,  the  absorption 

coefficient  0^  due  to  the  n-phonon  summation  process  peaks  near  naJ+  ^mx  ’ 

which  is  about  (3/4)  n  us  .  Thus,  the  total  absorption  coefficient,  which  is 

mx 

the  sum  of  all  0^  with  n  S  2,  is  dominated  at  a  particular  frequency  by  the 
nearest  n-phonon  summation  process.  The  frequency  dependence  of  0  in  the 
range  2  u:^^  ~  3  us^^  is  ,  therefore,  approximately  exponential  since  the 
strength  of  the  peak  in  0n  as  a  function  of  n  is  approximately  exponential.  The 
small  values  of  c^_  confirm  the  validity  of  neglecting  En_  in  Eq.  (3.7).  As  the 
temperature  is  raised,  the  strength  of  the  peak  in  0^  increases,  the  position  of 
the  peak  is  shifted  toward  the  lower  frequency,  and  the  width  is  either  narrowed 
or  broadened,  depending  on  the  temperature  dependence  of  the  phonon  frequencies. 
Thus,  8  increases  with  temperature,  but  less  rapidly  than  Tn  *  at  high  tem¬ 
peratures. 

Now  we  give  a  more  realistic  estimate  for  XaCl-structure  crystals.  There 

i  ~  i2 

are  several  wavevectors  for  which  the  explicit  expression  of  |  U  (Q)  |  can  be 


71 


Sec.  C 


written  down.  For  the  acoustic  branch  at  the  zone  boundary  in  the  (1, 1, 1) 
direction,  the  light-mass  ions  stand  still,  and  the  heavy  ions  move  in  the 


direction  w^q,  say.  Then, 


ux(Q) 


(x  •  w>0  )  m<  /m> 


(6.1) 


Similarly,  for  the  optical  branch  at  the  same  wavevector. 


I  UX(Q)  \2  =  (  x  •  w^  )2  . 


(6.2) 


For  the  optical  modes  near  the  zone  center,  Ux(Q)  is  given  by  Eq.  (2. 17).  For 
the  acoustical  modes  near  the  zone  center,  |  ( Q  )  |  ^  is  nearly  zero. 

i  i2 

As  a  rough  approximation,  |  U  (Q )  |  will  be  set  equal  to  zero  for 

A 

2 

co  <  f  to  where  f  <  1 ,  and  where  it  is  not  negligible,  j  U  ( Q )  |  will  be 

UlX  X 

approximated  by  an  average  of  the  three  known  expressions  (6.  1),  (6.2),  and 
(2.  17) ;  thus. 


U  (Q)|' 


(1  +  m</m>  )  6(u>0  -  fco^  )  , 


(6.  3) 


.  A  /v 

where  6  is  the  unit  step  function.  In  the  average,  we  have  replaced  (x  •  w<q) 

by  1/3,  which  is  the  value  for  the  ( 1,1,1 )  zone -boundary  mode  and  is  also  the 

2 


angular  average  of  cos  0  .  The  remaining  factor  is  the  average  of  1 ,  n^  /  m> 

and  /m  . 

C  '  r 


The  estimates  of  the  a's  are  then,  from  (3. 1),  (3. 2),  (3.  5),  (3.  6),  and  (6.  3), 

CO 

with  the  usual  approximation  of  the  sum  over  Q  by  6  (  d  to  g(co): 

*  n 


=  y  (1  +  m</m>  )  <(n+l)6/co  >  co^ 


(6.4) 
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a1+  =  <(n+l )  0  >  /  <  (n+1)  0  /co  >  co^  , 


(6.  5) 


a9,2  -  [  <(n+l)0w  >/<(n+l)  0/oj)  co^2  ]  -  a1+ 


(6.6) 


where  n  is  the  Bose-Einstein  distribution  factor , 


( A  0  )  =  \  d  co  g(co)  A  (to)  0  ( to  -  f  co  )  , 


(6.7) 


and  g (co)  is  the  phonon  density  of  states  normalized  to  unity.  Similar  estimates 


ocn_  <  a«.  /( 1  +  m</m>  )  . 


(6.8) 


The  density  of  states,  shown  as  the  solid  curve  in  Fig.  9,  is  approximated 


by  the  Debye  model , 


K<“>  = 


(6.9) 


sketched  as  the  dashed  curve  in  Fig.  9.  The  value  of  co  is  taken  as  the  Debye 


cut-off  frecuencv  in  (6.  9).  In  the  high -temperature  limit. 


n(co)  +  l  co rp  / co  + 


(6. 10) 


Then,  the  averages  in  Eqs.  (6.  4)  to  (6.  6)  are  easily  evaluated. 


Tiie  value  of  f  will  be  chosen  as  f  corresponding  to  the  assumption 


that,  for  1/8  of  the  modes  (1/4  of  the  acoustical  modes),  W  (Q)  is  negligible. 


In  Table  II ,  we  list  the  data  of  NaCl  along  with  the  values  of  the  a's  at  room 


temperature  corresponding  to  cOy/co^  =  1.03  for  NaCl. 


V 
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To  estimate  the  magnitude  of  the  vertex  correction,  we  need  to  know  the 
contribution  of  the  intermediate  phonon  in  the  form  of  £  given  by  Eq.  (5.  8). 

In  the  process  depicted  by  Fig.  4b,  the  intermediate  phonon  splits  into  two 
phonons  Q.  and  Q„  ,  which  were  taken  at  frequency  to  ,  i.  e. ,  in  the  optical 
branches.  This  is  reasonable  since  the  high-frequency  side  of  the  Gaussian 
P n curves  contribute  to  $  =  as  seen  in  Fig.  iO.  By  a  quasi-selection 

rule,10  Q4  must  be  an  acoustic  mode,  the  largest  contribution  of  which  will  be 
at  the  edge  of  the  Brillouin  zone.  Thus,  we  take 

2 

(6.11) 


U  (Q.  ) 
x '  M  7 


2 


1  m<  /  “Q4  \ 

T  mT  —  0<^ 

>  '  1  mx 


mx  '  WQ4 ) 


using  Eq.  (6.  1) .  The  frequency  of  the  highest  acoustic  mode  is  taken  to  be 

3  2 

^mx  wlt^  ^  =0.5,  and  the  factor  (tOq^  /iju^  )  approximately  simu¬ 

lates  the  effect  of  the  polarization  for  the  long-wavelength  acoustic  inodes. 
Substituting  Eq.  (6.11)  into  Eq.  (5.8),  averaging  over  the  possible  modes  of 
Q4 ,  and  summing  over  three  branches,  we  obtain  an  estimate  of  £ : 


€ 


,  A2)  2 

3  0  /5  m^,  to 

>  mx 


0. 18 


(6. 12) 


Substituting  this  value  into  Eq.  (5. 15)  yields  the  following  estimates  for  the 
vertex- renormalization  factors: 


A0  =  1  ; 


2  2 
A,  =  (1  +  0. 142 )  =  1.30  ; 


=  (1  +  0.388)  =  1.93  ; 


A 


5  =  ( 1  +  0.  844  T  =  3.  40 


A,  =  (1  +  1.72  f  =  7.38  . 


(6.  13) 
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We  note  that,  from  Eqs.  (5. 1)  and  (5.2),  the  vertex  ratio  can  be  shown  to 

be  equal  to  the  ratio  of  the  real  parts  of  the  self-energy  terms  given  by  Figs,  lib 

and  1 1  a  at  zero  temperature  and  frequency  2  to  .  From  R.  A.  Cowley's 
13 

calculation  for  KBr  ,  our  estimate  of  £  appears  to  be  somewhat  too  large. 

The  multiphonon  absorption  calculated  from  (3.  7)  using  the  values 
of  parameters  from  Eq.  (6. 13)  and  Table  II  is  shown  in  Fig.  10, 
where  the  individual  are  plotted  as  light  curves  and  the  sum  of  the  is 
plotted  as  the  heavy  curve.  The  agreement  is  rather  good  in  view  of  the  crude 
approximations  used  to  estimate  the  a's  .  It  should  be  noted  that  no  parameters 
have  beer  adjusted  in  the  theoretical  result. 

By  adjusting  two  parameters  in  Eq.  (3.7),  such  as  K  and  De  (keeping  the 
a's  at  f  =  4- ),  the  experimental  data  can  befitted  to  within  the  scatter  of  the 
data.  In  fact,  by  changing  only  the  value  of  the  single  interaction  strength 
parameter  1  /p  from  9.0  to  12  ,  the  dashed  curve  in  Fig.  10  is  obtained. 

This  larger  value  could  be  partly  explained  by  the  fact  that  the  higher -order 
anharmonic  coefficients  are  much  more  sensitive  to  the  shape  of  the  potential 
curve  than  the  quadratic  terms  from  which  the  value  of  p  is  determined. 

Errors  introduced  by  the  approximations  and  uncertainties  in  the  values  of 
the  parameters  used  also  could  account  for  the  difference,  of  course. 

The  near -exponential  frequency  dependence  is  evident  in  Fig.  10.  The 
vertex  correction,  which  is  included  in  Fig,  10,  slightly  ifnproves  the  agree¬ 
ment  with  the  experimental  result.  Without  this  correction,  the  ^ "jSg  curves 
would  be  shifted  down  by  factors  of  1 ,  1.  3,  1.  9,  and  3.  4,  respectively. 
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Figure  11.  Phonon  self-energy  terms  of  order  e 


Table  II.  Values  of  parameters  for  NaCl  at  room  temperature 


dynes /cm' 


1.50  (formally  for  all  u>) 
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The  n-phonon  regions,  marked  on  Fig.  10,  do  not  correspond  to 


n  tcf  <  co  <  ( n+ 1 )  oof  ,  to  n  <  in  <  ( n+  1 )  10^ 


or  to  n  to  j^q  <  to 


<  (n  +  1 )  to  j^q  ,  as  is  often  assumed  in  the  literature.  In  fact,  the  n-phonon 
regions  shift  as  the  temperature  changes,  as  discussed  below. 


The  n  =  2  central-limit  curve  is  included  in  Fig.  10  even  though  its  ac¬ 
curacy  is  not  expected  to  be  good.  The  two -phonon  structure  is  lost,  of 
course,  in  approximating  jS2  by  a  Gaussian,  and  the  peak  does  not  occur 


at  CO  =  to- . 


The  temperature  dependence  of  0  at  a  given  frequency  in  the  nearly  expo- 

25 


nential  region  is  considerably  weaker  than  that  of  the  simple  expression 

-n 


/  —  to / to t"1  \  /  -to/ntoT  \ 

jSn(T)/0n(O)  =  (  1  -e  |^1  -e  ~T 


n-1 


(6. 14) 


obtained  formally  from  the  occupation -number  factor  (2.4)  by  setting  all 

n  “  1 

OOq  =  oo/n.  The  approximation  0n  ~  T  in  (6. 14)  is  valid  in  the  high- 
temperature  limit,  and  n  has  been  assumed  to  be  independent  of  temperature 
in  the  past.  The  T  dependence  of  /3  results  from  the  temperature  dependence 
of  the  parameters  a,e*,  and,  particularly,  the  phonon  frequencies  C0q  and 
from  the  explicit  temperature  dependence  of  ,  a1+,  and  0i2+. 


The  following  example  of  NaCl  at  300K  and  10.  6fim  illustrates  the  strong 

,n-l 


devj.ai.ior  from  the  frequently  quoted  result  /3  T  .  The  value  of  the  slope 
(T//3)d|3/dT  of  /3  as  a  function  T  on  a  log- log  plot  can  be  estimated  from 


Eq.  (3.7).  Using  n  =  5.  5  from  Fig.  10  and  the  following  approximate  expres- 


26,27 

sions  for  the  temperature  dependence  of  the  parameters,  cc 


Q 


=  co^(l-3.8x  10'4T),  a  =  aQ(l  -4.4x  10_5T),  ande*  =  e^  ( 1  -  1.06x  10‘4T), 


'QO 
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we  find 

(T/j3)dj3/dT  =  2.5  (6.15) 

which  is  considerably  smaller  than  n  -  1  =  4.5.  The  Bom -Mayer -potential 
parameters  C  and  pK  =  pa  in  (2.11)  are  essentially  temperature  independ¬ 
ent,  being  electronic  in  nature.  In  particular,  (2. 12)  should  not  be  used  to 
ascribe  a  temperature  dependence  to  p^  from  measured  values  of  B(T)  and  a(T). 
The  temperature  dependence  of  B  arises  from  anharmonic  and  volume  effects, 
not  from  a  T  dependence  of  p^  . 

A  weakening  of  the  temperature  dependence  such  as  that  in  (6. 15)  is  ap- 

25  8 

parent  in  the  data  of  Harrington  and  Hass,  Barker,  of  Kaiser  and  co- 
8  8 

workers,  and  of  Denham  and  coworkers.  Finally,  it  is  mentioned  that  in 
a  material,  possibly  a  zinc -blende -structure  crystal,  in  which  the  position 
of  a  given  multiphonon  peak  can  be  traced  as  a  function  of  temperature,  the 
temperature  dependence  should  be  quite  different  from  that  of  |3  at  a  given 
frequency.  A  detailed  presentation  of  the  temperature  dependence  of  £  will 
be  given  elsewhere. 

The  /3-  cc  curves  of  the  alkali  halides  and  alkaline  earths  show  less 
structure  than  those  of  the  semiconductor  materials.  It  is  plausible  that  the 
greater  anharmonicity  of  the  NaCl -structure  crystals  could  give  rise  to  such 
short  lifetimes  of  the  zone -boundary  phonons  that  the  fine  structure  in  the 
density  of  states  is  essentially  eliminated. 

13  28 

The  lifetime  of  the  fundamental  phonon  is  short,  '  1  and  the  lifetimes 
of  the  zone-boundary  phonons  should  be  even  shorter  since  the  selection  rules 
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and  quasi-selection  rules  do  not  apply  to  the  zone-boundary  modes  (with  nonzero 
wavevectors).  A  value  of  relative  linewidth  2T/u :  of  the  order  of  0.  3  for  the 


zone-boundary  phonons  at  resonance  should  be  sufficient,  and  this  value  is  reason- 

28 


able  in  view  of  the  value  of  2170.'  =  0.07  for  the  fundamental  mode  in  NaCl  and 
the  fact  that  2r/o  is  expected  to  be  larger  at  the  zone  boundaries.  Furthermore, 
as  n  becomes  larger,  more  convolutions  are  involved  |  seeEq.  (3.3)],  and  each 
convolution  tends  to  smooth  out  any  fine  structure  in  the  density  of  states. 


This  explanation  is  consistent  with  the  experimental  results  which  show  that 


the  two -phonon  peaks  are  wider  in  the  alkali  halides  than  in  the  semiconductor 


materials,  that  T~T  at  the  fundamental  resonance  in  NaCl  (implying  that  the 


two-phonon  contribution  is  small  at  resonance),  and  that  the  two-phonon  peaks 

,2 


have  been  observed  in  NaCl  even  though  T~T  at  resonance.  A  careful  study  of 
the  temperature  dependence  of  the  two -phonon  summation  peaks  could  show  an 
increase  in  the  widths  of  the  peaks  as  T  is  raised  from  77  K  to  the  highest  prac¬ 


tical  temperature  of  the  solid.  Such  increases  are  apparent  in  the  small  amount 
29 


of  existing  data.  As  the  temperature  is  reduced  below  room  temperature,  addi¬ 
tional  multiphonon  peaks  could  appear  in  higher-n  regions  where  /3 (co )  is  relatively 

smooth  at  room  temperature.  Of  the  three  existing  known  cases  (for  CaF2  ,  BaF  , 

29 

and  SrF^  at  77  and  300  K),  two  show  a  small  additional  peak  at  77  K  that  is  ab¬ 
sent  at  300  K . 


It  should  be  emphasized  that  the  two -phonon  peaks  are  associated  with  peaks 
in  the  appropriate  density  of  states  and  are  not  resonance  lines.  Thus,  an  extrapo¬ 
lation  of  /3(eo)  from  the  reststrahl  region  should  not  be  subtracted  from  0  at  higher 
frequencies  to  obtain  the  multiphonon  contribution,  as  is  sometimes  done  in  the 
literature.  An  alternate,  though  unlikely,  explanation  of  the  lack  of  structure  is 
that  the  raw-phonon  density  of  states  shows  little  structure. 
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VII.  ASSUMPTIONS  AND  APPROXIMATIONS 

The  assumptions  and  approximations  made  in  the  previous  sections  are  now 

summarized:  (1)  The  photon-phonon  coupling  is  given  by  the  leading  dipole  term, 

9 

and  the  Lax-Burstein-Bom  mechanism  is  neglected.  (2)  For  the  anharmonic 

forces,  only  the  nearest-neighbor  Born- Mayer  repulsion  term  is  included  and  is 

further  approximated.  (3)  The  lifetimes  of  the  intermediate  and  final-state  phonons 

are  assumed  to  be  infinite.  (4)  The  central-limit  theorem  is  used  to  reduce  the 

n-fold  multiple  sum  in  (2.3)  to  a  Gaussian  whose  parameters  a  are  given  by  single 

20 

sums,  although  it  is  possible  to  improve  the  asymptotic  approximation.  (5)  Rough 
estimates  were  given  for  the  various  Brillouin-zone  sums  over  the  phonon  coordi¬ 
nates  for  the  coefficients,  a's  .  All  of  these  approximations  except  (1)  were  shown 
to  be  reasonable.  The  perturbation  approach  used  is  justified  by  showing  that  all 
diagrams  not  included  in  the  results  are  negligible.  Concerning  (1),  the  long¬ 
standing  question  of  the  importance  of  the  Lax-Burstein-Bom  mechanism  in  NaCl- 
structure  materials  remains  unanswered.  The  mechanism  is  quite  simple  to  in¬ 
clude  formally;  estimating  the  strengths  of  vertices  has  been  the  problem. 

Our  calculation  gives  good  agreement  with  experimental  results  for  the  fre¬ 
quency  dependence  of  the  optical  absorption  and  demonstrates  the  general  nature 
of  this  dependence  for  crystals  with  tetrahedral  symmetry.  The  estimates  listed 
in  (5)  above  enable  us  to  see  explicitly  the  nature  of  our  results.  Some  of  the  esti¬ 
mates  must  be  regarded  as  tentative.  However,  these  approximations  are  not  es¬ 
sential  to  our  theory.  We  plan  to  perform  both  the  multiple  sums  in  (2.21)  fr»r 
n  =  2  -  6  and  the  single-phonon  sums  in  Eqs.  (3.2),  (3. 5),  (3.  6),  and  (5.  8)  by  com¬ 
puter.  This  will  enable  us  to  examine  more  rigorously  the  validity  of  the  other 
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approximations,  especially  (4).  The  computer  results  for  the  multiple  sums  in 


(2.21)  should  provide  greater  accuracy  in  the  small-n  regions,  say  n  =  2  and  3, 


where  the  central-limit  results  are  less  accurate,  and  should  afford  a  good  test 


of  the  approximations  in  the  region  of  n  =  4  -  6 .  The  temperature  and  frequency 


dependence  of  ft  for  a  number  of  crystals  will  be  included  in  the  computer  pro¬ 


gram,  which  is  being  performed  in  collaboration  with  A.  Karo  of  the  Lawrence 


Livermore  Laboratory. 
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D.  TEMPERATURE  DEPENDENCE  OF  MULTIPHONON  INFRARED  ABSORPTION* 

M.  Sparks 

Xonics,  Incorporated,  Van  Nuys,  California  91406 

and 

L.  J.  Sham 

University  of  California,  San  Diego,  La  Jolla,  California  92037,  and 
Xonics,  Incorporated,  Van  Nuys,  California  91406 

Measurements  of  Harrington  and  Hass  and  of  Barker  indicate  that 
the  temperature  dependence  of  the  infrared  absorption  coefficient 
in  the  n-phonon  region  is  considerably  weaker  than  ~  Tn  \  which 
had  been  predicted  for  the  high-temperature  limit  of  multiphonon  ab¬ 
sorption.  This  discrepancy  is  resolved  by  taking  into  account  the 
temperature  dependence  of  the  phonon  frequencies  and  the  lattice 
constant.  The  agreement  between  the  experimental  and  theoretical 
results  with  no  adjustable  parameters  is  good.  A  new  evaluation  of 
the  multiphono .t  sums  yields  0  ~  exp  (-to  T  )  directly,  rather  than  as 
a  sum  on  n . 

*  1 

The  nearly  exponential  frequency  dependence  of  infrared  absorption  in  the  region 

of  low  absorption  has  been  explained  recently  by  a  simple  multiphonon-absorption 

2  “5  6 

theory  °  and  by  independent-molecule  models.  The  frequency  and  temperature 

dependence  of  the  optical  absorption  coefficient  |9  are  of  fundamental  interest,  and 

both  should  be  useful  in  identifying  intrinsic  absorption  and  in  distinguishing  between 

intrinsic  and  extrinsic  absorption.^’  ^  Harrington  and  llass'^  have  shown  that  the 

2-4 

temperature  dependence  of  /3  is  considerably  weaker  than  the  expected  dependence 
jS  [1  -expJ-to/oCrpH  1  -exp(-iL'/nu:T)J  ns:Tn  1 


for  n-phonon  absorption. 
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The  approximate  equality  is  the  high-temperature  limit,  to  is  the  lasei  frequency, 

9 

and  tOrp  =  k^T/fi.  Reexamination  of  earlier  data  reveals  similar  discrepances, 
which  constitute  the  most  serious  problem  in  the  recent  developments  in  the  theory 
of  multiphonon  absorption. 

These  discrepancies  are  explained  by  including  the  temperature  dependence  of 
the  phonon  frequencies  U!q  and  lattice  constant  in  our  previous  theory.  A  simple 
estimate  indicated  previously  that  the  resulting  deviations  would  be  quite  large.  111 
The  previous  expression^  4’  111  for  |3  is 

0  =  f(a.')u:'4  E  (An2vn/n!)  £•••£  6  fcc-T  )  n  an  (1) 

n=2  Qj  Qn  \  j-1  Qj  /  £=1  Q£ 


where  f  (w)  =  con.  [  1  -  exp  ( -u:/wT )  ] ,  v  =  Ti /2p„2  m,  an  =  W„  (n„  +1)/Na*n  , 

1  K  Q£  Q£  Q£  Q£ 
Wq  is  of  order  unity  for  large  u.'q  and  is  very  small  for  small  u)~  ,  nn  are 

£  W£  h  *£ 

phonon  occupation  numbers,  the  Q's  denote  wavevectors  and  branches,  m  is  the 

reduced  mass,  pK  is  the  Born- Mayer  repulsive-potential  parameter  (cn),  2N  is 

the  number  of  ions  in  the  crystal,  and  the  higher-order  terms  in  the  pe  turbation 

expansion  give  rise  to  the  vertex-correction  factors4, 10  A  =1  +  A  f+Hf2), 

n  n  ’  ’ 

where  A4  1.  94,  A,  =  3.  92,  A^  =  7.  15,  and  ^  =  9BQaQ  /5  ( 1  -  2  pQ  )  m  co  Here 
B  is  the  bulk  modulus,  p  =  pK/a,  the  subscript  0  denotes  T=  0,  and  to  is  a 
frequency  near  the  top  of  the  phonon  spectrum.  Kq.  (1)  can  be  written  dowr  immedi¬ 
ately  from  the  well  known  expression11,  10,2*4  for  p  with  ^  »  (reststrahl  fre- 

2  1  O  n  j 

fluency,  cc f  )  and  standard  perturbation -theory  results,  apart  from  the  details  ' 
of  An  and  cCq  which  are  not  needed  here. 
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See  1) 


10  ? - 4  v 

It  is  simple  to  show  tli.it  for  n>2  the  central-limit  approximation  to  the 

multiple  sums  on  Qj*--  Q  is  satisfied,  thus  reducing  (1)  to  0  ~  &n  vV*t'1 


/3n  =  f  ( tc )  tc  f  (uvC) 


ex)  (-  a  /p^  )  G/a 


3 


where  av  denotes  the  weighted  (with  \V())  average  o  'er  phonons,  a  is  the  near 

neiglibor  distance,  and  G/v'n"'  is  the  Gaussian  from  the  central-limit  approximation.  ’ 

The  small  temperature  dependence  of  the  effective  charge  was  neglected,  and 

A  ^/n  1  /fT  was  approximated  by  C11,  which  is  accurate  over  the  range  of  n's  used 

here  (typically  n  =  3-6),  as  seen  by  the  straight  line  obtained  by  plotting 
2 

in  [  An  /  n  I  vTT  ]  vs  n.  The  two  parameters  in  the  Born- Mayer  potential  are  essen¬ 
tially  independent  of  temperature,  heing  electronic  in  nature.  'Hie  standard  relation 
2  -1 

B  =  con.  a  Pj,  -2,  with  experimental  values  for  the  temperature  dependence  of  B  and  a, 

should  not  be  used  to  determine  the  T  dependence  p^  since  the  T  dependence  of  B  and  a 

arises  from  anharmonic  and  volume  effects,  not  from  the  temperature  dependence  of  p,, . 

K 


The  average  in  (2)  is  easily  evaluated  by  using  the  Taylor's  expansion 
Hq  +  1  ^  (ui^/oCq  )  +  j  +  (tCq  /  12  uirj,  ),  which  is  well  satisfied  for  all  cases  con¬ 
sidered,  and  a  Debye  spectrum  (with  cut-off  frequency  u.’nu( )  truncated  at  ^  cc 

to  account  for  the  angle  factor  in  the  vertex  [see  in  (1)].^  ^  This  gives 

2  *  3 

u  =  9fi(n  +l)/16p,,  xni,  with  x  =  3a.'  /4.  Tlie  factor  of  -r  is  obtained  when 

x  rK  mx  4 

the  three  terms  in  the  expansion  are  recombined  after  integration. 


A  simple,  accurate  expression  for  the  slope  (T//3)dj3/dT  of  in  ^  vs  £nT, 
which  is  a  convenient  and  sensitive  measure  of  the  T  dependence  of  0 ,  is  obtained 
by  approximating  j9  by  the  locus  of  the  inflection  points  on  the  high -a;  sides  of  the 
peaks  of  the  Gaussians.  The  analysis  leading  to  (4)  below  or  inspection  of  the 


n  i 
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curves  and  their  sum2"4*  10  show  the  validity  of  this  approximation.  In  ,hc 
absence  of  detailed  information  on  the  temperature  dependence  of  the  toQ .  it  is 
assumed  that  (T/^ld^/d  r  -  (T/c,f )  d  *f/dT  .  Then  (!)  gives 


T  dfi 

7T  77  =  +  D  +  D  , 

P  d  T  1  to  a 


D  = 

tO 


D  = 


;( 


/  •*  to 

\  [  Xn  C  v  u  + - !H*.  n 

\  4  a.'  3 

'  T 

4  mx 


1  + 


1  3 

4  k'  mx 

tO 

n 

“)  ’ 

A„€-i  " 

T 

du 

l 

HA  £ 
n 

_ 

to 

f 

d  T  ’ 

3  +  -L  \  I  liL 
P0  /  a  dT  * 


The  results  fro,,,  (.3)  for  NaP,  UP.  NttCl,  and  KBr,  for  which  experimental 
ues  arc  available,  ar.  -sted  rn  Table  I  along  with  the  values  of  input  parameters 
and  experimental  values  of  Harrington  and  Hass8  and  Barker.9  The  agreement  is 
surprisingly  good,  probably  fortuitously  so.  tn  view  of  the  uncertainties  in  the  ex¬ 
perimental  values  and  the  fact  that  the  ccQ  (T)  are  no,  well  known.  Note  that  there 
are  no  adjustable  parameters  in  Ihe  theoretical  result,  and  that  the  value  of  D  is 
sensitive  to  small  errors  In  the  values  of  the  input  parameters.  Fig.  1  shows  thl  excel¬ 
lent  agreement  bet..een  theoretical  curves  and  the  experimental  points  of  Harrington 
and  Hass  and  of  Barker  and  illustrates  the  greater  sensitivity  of  the  parameter 
<T/0)d0/dT  than  of  such  a  plotted  comparison,  the  10  percent  difference  between 
experiment  and  theory  fot  KBr  a,  1000  K  appearing  smaller  than  the  corresponding 
40  percent  difference  in  the  values  of(T//3)d/3/dT. 

In  order  to  obtain  the  experimental  intrinsic  value  of  (T/0)  d0 /dT  for 
NaCl  at  943  cm  ,  which  is  obscured  by  the  extrinsic  contribution,  the  ex¬ 
trinsic  contribution  was  assumed  to  be  independent  of  temperature,  and  the 
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Table  Caption 

Table  I.  Theoretical  and  experimental  values  of  ( T//3)dj3 /dT.  Experimental  values 
marked  ~  and  «  are  uncertain  (~50%)  and  highly  uncertain  (factor  of  ~  3). 

a.  A.  M.  Karo  and  J.  R,  Hardy,  Lawrence  Radiation  Laboratory  Report 
UCRL-14822,  April  (1966). 

b.  American  Institute  of  Physics  Handbook,  2nd  Ed. 

c.  J.  P.  Jasperse,  A.  Kahan,  J.  N.  Plendl,  and  S.  S.  Mitra,  Phys.  Rev.  H6, 

526  (1966). 

d.  I.  F.  Chang  and  S.  S.  Mitra,  Phys.  Rev.  B5^,  4094  (1972). 

e.  J.  E.  Mooij ,  W.  B.  VanDeBunt,  and  J.  E.  Schrijvers,  Phys.  Letters  28  A, 

573  (1969). 


Table  I.  Theoretical  and  experimental  values  of  O/0)d/3/dT.  Experimental  values 
marked  ~  and  **  are  uncertain  (~50%)  and  highly  uncertain  (factor  of  ~3). 
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Figure  Caption 

Figure  1.  Comparison  of  experimental  points  of  Harrington  and  Hass 
(NaF,  943  cm"1 )  and  Barker  (KBr,  418  cm  1 )  with  theoretical  curves  fit 
to  the  data  at  300  K  . 
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intrinsic  value  of  0  at  300  K  was  assumed  to  be  equal  to  7x  10"4cm_1.  *  These 
]uestio  utble  assumptions  surely  introduce  considerable  errors.  Indeed,  the  large 
experimental  value  of  (T/ 0)  d/J/ dT  =  6  at  300 K  and  the  equal  values  of  2.4  at 
700  and  943  cm  1  appear  to  be  unreasonable. 

In  the  classical,  high-temperature  limit  of  »  to,  all  occupation  numbers 
nu)Q  are  rePlaced  by  wT/tOg.  For  LiF  at  1800cm"1  and  300K  the  classical  ap¬ 
proximation  to  (T/0)(d0/dT)  from  (3)  is  approximately  three  times  greater  than 
the  unapproximated  value.  Even  at  1000  K,  the  classical  limit  gives  3.5,  compared 
with  2.4  from  theory  and  2.5  for  experiment.  Thus,  considerable  care  must  be  ex¬ 
ercised  in  applying  classical  theories  to  practical  cases,  as  already  discussed  by 
Maradudin  and  Mills  in  the  accompanying  Letter. 

In  a  material,  such  as  AlSb  for  example,  in  which  the  position  of  a  given  multi¬ 
phonon  pe.  tc  can  be  traced  as  a  function  of  temperature,1'1  the  temperature  dependence 
should.be  quite  different  from  that  of  at  a  given  frequency,  since  following  the  peak 
eliminates  the  contribution  from  dn/dT. 

The  exponential  frequency  dependence  of  /3  can  be  obtained  directly ,  rather  than 

as  a  sum  over  n ,  as  follows:  Representing  the  delta  function  by  an  integral  over  t 

reduces  (1)  to  0  =  f(w)w'4  j(dt/2ir)  exp(iwt)  £  n4  [ g ( t ) ] ",  where  g(t) 

=  Dq  CTQ  cxp("  1  tore  An2  vn/n  n!  was  approximated  by  Dn,  which  is 

quite  accui  ate  for  n  -  3-<S.  The  sum  on  n  can  be  written  as  a  linear  combination 
"1  n 

of  ( 1  -  g)  =  Eg  and  its  first  four  derivatives.  The  contour  integrals,  which  have 
poles  at  t  =  i  T  ,  arc  easily  evaluated  by  residues,  giving 

-OCT 

&  ~  Pq  C  (4) 
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where14  ~  f(a',)[dg(iT)/dT  ]  ^  and  T  is  the  solution  to  l-g(ir)=0,  i.e. , 

1  -DN  1  £  Wn(nn+1)  con  1  exp  (to n  r)  =  0,  which  is  easily  solved  numerically. 15 
q  W  x  v  v 

The  following  approximation  illustrates  the  general  dependence  of  r  on  T,  An,  etc. , 
although  it  is  too  crude  to  afford  accurate  values  of  T .  Neglecting  the  angle  depend¬ 
ence  of  Wq  and  app,  ixtmating  the  density  of  states  by  6(o  -  co^)  gives 
T  =  -co 1  Xn[6 C Wb(nj3  + 1 ) /co^  ]  and  dg(iT  )/dT  =  co^ . 


The  temperature  dependence  o'  0  from  (4)  is  formally  the  same  as  that  from 
(3).^  The  numerical  values  will  differ  slightly,  corresponding  to  the  slight  numeri¬ 
cal  difference  between  r  in  (4)  and  the  effective  T  from  (2). 


We  would  like  to  thank  Dr.  J.  A.  Harrington  and  Dr.  M.  Hass  for  sending  their 
data  prior  to  publication.  Discussions  with  Dr.  A.  A.  Maradudin  and  Dr.  D.  L.  Mills 
are  gratefully  acknowledged.  Dr.  C.  J.  Duthler  kindly  assisted  with  the  calculations. 
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Abs  true t 

have*  calculated  the  frequency  dependence  of  infrared 
absorption  in  the  classical  limit  for  an  exactly  soluble 
model  of  a  lattice  of  noninteracting  diatomic  molecules, 
each  bound  internally  by  a  potential  for  which  the  classical 
equation  of  motion  can  be  solved  in  closed  form.  Four  poten¬ 
tials  have  been  used:  a  Morse  potential,  a  potential  of  the 
form  V(x)  =■-  (a/x  )  +  bx  ,  an  infinite  square  well  potential, 
and  a  triangular  well  potential.  The  analytic  results  we  ob¬ 
tain  show  that  the  absorption  coefficient  for  large  frequen¬ 
cies  associated  with  potentials  which  admit  an  harmonic 
approximation  decreases  nearly  exponentially  over  the  fre¬ 
quency  region  covered  by  recent  experiments,  with  significant 
deviations  from  exponential  behavior  at  higher  frequencies. 
For  the  square  and  triangular  well  potentials,  the  absorption 
decreases  as  uj  for  frequencies  large  compared  to  a  char- 


actexistic  frequency. 
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I  .  I  ntrochiction 

The  absorption  of  electromagnetic  radiation  by  { he  lattice 
vibrations  in  anhannonie  crystals  has  received  considerable 
attention  from  both  theorists  and  experimentalists  for  many 
years.  However,  most  studies  have  focused  attention  on  frequen¬ 
cies  either  in  the  near  vicinity  of  the  fundamental  roststrahl 
absorption  bands,  or  at  frequencies  sufficiently  low  that  the 
dominant,  portion  of  the  absorption  may  bo  accounted  for  by  pro¬ 
cesses  which  involve  at  most  two  phonons.  ^ 

Recently,  interest  has  been  aroused  in  the  behavior  of  the 
absorption  coefficient  at  frequencies  several  times  (say  2  to  10 
times)  the  maximum  vibrational  frequency  of  the  crystal,  but  still 
small  compared  to  the  electronic  band  gap.  In  this  frequency 
region,  the  principal  contribution  to  the  absorption  coefficient 
in  a  pure  crystal  presumably  comes  from  multi  phonon  processes, 
where  the  number  of  phonons  involved  may  be  quite  large.  The 
behavior  of  the  absorption  coefficient  in  this  frequency  regime 
is  clearly  important  to  understand  for  fundamental  physical 
reasons.  There  is  also  a  great  deal  of  practical  interest  in 
this  region,  since  high  power  CC>2  lasers  produce  intense  beams 
of  radiation  at  10.6jj  .  This  corresponds  to  a  frequency  several 
times  that  of  the  maximum  vibrational  frequency  of  many  materials 
ti  at  may  prove  useful  for  the  fabrication  of  windows  and  lenses 
for  use  with  these  devices.  Because  the  radiation  from  these 
lasers  is  very  intense,  even  a  small  amount  of  absorption  can 
lead  to  appreciable  heating  of  any  window  through  which  the*  beam 
passes.  It  is  therefore  of  interest  to  understand  the  nature  of 
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the  Intrinsic  .absorption  processes,  as  well  as  impurity  and 
surface  induced  absorption  at  frequencies  high  compared  to  the 
characteristic  vibrational  frequencies  of  the  crystal. 

One  may  readily  come  to  appreciate  the  difficulty  of  carry¬ 
ing  out  a  first  principles  calculation  of  the  frequency  depend¬ 
ence  of  the  absorption  coefficient  in  the  mul  tiplionon  regime, 
for  a  realistic  model  of  an  anharmonic  crystal  lattice.  What 
is  quite  intriguing  is  that  experimental  studies  of  the  frequency 

dependence  of  the  absorption  coefficient  in  several  alkali  halide 
(2) 

crystals  have  revealed  that  in  all  the  cases  studied,  for  fre¬ 

quencies  in  the  region  of  200  cm"1  to  800  cm"1  the  absorption 
coefficient  at  room  temperature  may  be  fitted  quite  accurately 
by  the  empirical  formula 


P(uu)  =  A  exp(-Biu)  , 


(1-1) 


whore  p(oj)  is  the  absorption  coefficient  at  frequency  u),  and  A 
and  B  are  constants  charr iteristic  of  the  particular  crystal.  It 
is  extremely  important  to  know  whether  Eq.(l)  can  be  derived  from 
a  theoretical  model  of  some  generality,  and,  if  so,  j t  is  import¬ 
ant  to  know  if  it  holds  to  frequencies  as  high  as  10. 6p  (which 
lies  outside  the  range  accessible  to  the  experimental  studies), 
and  also  if  it  holds  at  temperatures  higher  than  room  temperature. 

In  this  paper,  we  wish  to  address  ourselves  to  these  ques¬ 
tions.  Because  of  the  difficulty  of  carrying  out  calculations  of 
the  absorption  coefficient  in  the  multiphonon  regime  that  are 


both  i  o.ilis  tic  and  that  lead  to  conclusions  of  a  general  nature, 
wo  have  chosen  to  exploit-  the  properties  of  a  model  of  a  solid 
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that  is  highly  schematic,  but  v/hicli  allows  simple  analytic  expres¬ 
sions  to  be  obtained  for  the  absorption  coefficient  for  a  variety 
of  interatomic  potentials.  We  replace  a  diatomic  solid  which  con¬ 
sists  of  N  unit  cells  by  an  array  of  N  electric  dipole  active, 
but  anharmonic  oscillators.  While  such  a  model  is  rather  over¬ 
simplified  if  we  choose  to  represent  a  real  solid  by  it,  by  an 
examination  of  the  model  we  can  gain  insight  into  the  question  of 
whether  the  form  in  Eq.(I-l)  is  valid  quite  generally.  If  it  is 
valid  quite  generally,  it  should  also  be  valid  for  our  model.  If 
a  realistic  potential  is  chosen  for  the  anharmonic  oscillator,  we 
feel  the  model  also  provides  a  reliable  semi-quantitative  estimate 
for  the  magnitude  of  the  absorption  in  the  multi-pbonoi.  regime. 

On  the  basis  of  our  model,  we  will  also  be  led  to  the  conjecture 
that  at  high  frequencies,  the  magnitude  of  the  absorption  coeffi¬ 
cient  might  be  quite  sensitive  to  the  presence  of  certain  impuri¬ 
ties  . 

Since  the  region  of  experimental  interest  to  date  is  room 
temperature  and  above,  we  have  used  the  methods  oi  classical 
phjsics  to  compute  the  absorption  coefficient.  We  obtain  a 
general  expression  for  the  absorption  coefficient  lor  the  model 
described  above,  and  then  apply  the  expression  to  the  study  of 
the  frequency  dependence  oJ  the  absorption  coefficient  for  four 
potential  functions.  We  consider  absorption  by  anharmonic 
oscillators  described  by  the  Morse  potential,  a  second  potential 
which  posse  ses  a  hard  core  and  admits  a  harmonic  approximation 
(\  'x)  b  x"  i  a  x  tlit'  square  well,  and  a  potential  of 

trj angular  shape. 
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The  outline  of  the  papei  is  as  follows.  In  section  II,  we 
obtain  a  general  expression  for  the  absorption  coefficient  of  the 
oscillator  array  by  the  use  of  the  methods  of  classical  statis¬ 
tical  mechanics.  In  section  III,  we  apply  this  expression  to  the 
four  examples  mentioned  in  the  preceding  paragraph.  In  section  TV, 
we  present  a  discussion  of  some  implications  of  the  results  ob¬ 
tained  in  section  III. 

(3) 

Quite  recently,  McGill,  Hellwarth,  Mangir  and  'Yinston 
have  also  presented  a  theoretical  discussion  of  multiphonon  ab¬ 
sorption  by  an  array  of  uncoupled  oscillators.  In  the  body  of 
this  paper,  these  authors  present  a  diagrammatic  calculation  of 
the  absorption  coefficient  which  they  argue  leads  to  an  exponen¬ 
tial  form  identical  to  that  displayed  in  our  Eq.(I-l),  for  a 
specific  model  of  the  interatomic  potential.  Their  model  presumes 
that  in  the  crystal  Hamiltonian  the  term  proportional  to  the  n 
power  of  the  atomic  displacements  is  proportional  to  only  the 
quantity  gn  V n 1  ,  where  g  is  independent  of  n,  and  the  factor  n ! 
apparently  comes  from  the  Taylor  series  expansion  of  the  crystal 

potential.  This  model  is  quite  special,  since  one  may  find  a 
large  variety  of  realistic  potentials  which  admit  an  harmonic 
approximation,  and  ior  which  their  factor  gn  will  be  replaced 
b}  a  quan 1  i  t  y  that  exhibits  a  fundamentally  different  dependence 
on  n.  (Consider  the  Lennard-Joncs  G-12  potential,  or  any  poten¬ 
tial  which  contains  a  term  which  varies  inversely  with  a  power 
of  the  interatomic  separation.)  Their  conclusion  that  the  ab¬ 
sorption  eoeifici'i’t  varies  exponential  ly  witli  frequency  follows 
upon  counting  the  number  of  important  diagrams  in  the  first  lew 


If 
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orders  of  perturbation  theory.  In  view  of  the  discussion  in  the 
following  oaragraph  of  the  present  paper  and  also  that  n  cur 
section  IV,  we  are  led  to  question  the  conclusion  that  lor  a 
potential  of  general  form,  the  theory  produces  an  analytic  expres¬ 
sion  for  the  absorption  coefficient  exponential  in  character  in 
the  multi  oiionon  regime. 

(3) 

In  Appendix  A  of  their  paper,  McGill  et  al .  consider  an 
array  of  non-interacting  oscillators,  each  of  which  ir.  described 
by  a  Morse  potential.  They  insert  the  expression  exhibited  by 
Heaps  and  Hcrzberg^^  for  the  appropriate  electric  dipole  moment 
matrix  element  into  the  quantum  mechanical  form  for  the  absorp¬ 
tion  coefficient.  A  simple  analytic  expression  for  the  absorp¬ 
tion  coefficient  of  the  model  follows  from  this  procedure.  If  u> 
is  tho  maximum  vibration  frequency  of  the  crystal,  then  when 
hpT  >  the  correspondence  principle  applies,  and  their  expres¬ 
sion  may  he  compared  with  the  result  we  obtain  below.  The  two 
results  agree  in  this  regime.  However,  while  both  results  provide 
a  rather  good  lit  to  the  loom  temperature  data  in  the  regime  of 
frequencies  explored  by  Dcutsch,  and  they  thus  appear  qualita¬ 
tively  consistent  with  the  form  in  Eq.(l-l),  at  higher  frequencies 
significant  deviations  are  predicted  by  both  expressions.  As  the 
temperature  is  increased,  these  deviations  are  expected  to  set  in 
at  progressively  lover  frequemics.  Nci  thoi  our  calculation  nor 
that  presented  in  the  Appendix  of  the  paper  of  McGill  <  l  al., 
produces  an  analytic,  expression  for  thr  absorption  coef  1  i  c  i  en  t 
of  an  array  of  independent,  oscillators  which  exhibits  the  exponen¬ 
tial  behavior  suggested  l»\  l.q.(l-j),  although  as  remarked  above, 
the  quantitative  di  lioivners  arc  small  in  l  ’»« •  frequency  rep.  i  me 


17 


Sec.  E 


explored  by  the  experiments  so  far.  Thus,  while  one  may  con¬ 
struct  a  particular  potential  that  leads  to  something  close  in 
form  to  an  exponential  law,  we  feel  that  quite  generally,  th( 
exponential  law  is  not  valid,  and  one  may  fit  the  data  quite 
well  by  the  forms  we  obtain  below. 
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II.  General  Theory 

S^nce  we  consider  the  crystal  to  be  an  array  of  non-interacting 
Molecules,  infrared  absorption  by  the  collection  of  oscillators 
will  bo  N  times  that  of  a  single  molecule,  whore  N  is  the  number 
of  molecules  in  the  crystal.  Thus,  in  what  lollows,  we  consider 
only  the  absorption  by  a  single  molecule.  For  this  the  motion  of 
the  center  of  mass  of  the  molecule  is  irrelevant,  since  it  makes  no 
contribution  to  its  dipolo  moment.  Consequently,  the  equation  of 
motion  which  provides  the  starting  point  for  our  treatment  is 


mit  -- p  *= 


-  dVJiO  +  dM(x) 

dx  dx  v  ; 


(II-l) 


where  m  is  he  reduced  mass  of  the  molecule,  x  is  the  relative 
coordinate  of  the  two  atoms  comprising  the  molecule  with  p  the 


momentum  canonically  conjugate  to  x,  V(x)  is  the  interaction  poten¬ 
tial  energy  between  these  two  atoms,  M(x)  is  the  dipole  moment  of 
the  molecule  and  E(t)  is  the  electric  field  of  the  incident  infra¬ 
red  rrdintxon. 

To  obtain  the  rate  at  which  energy  is  absorbed  by  the  molecule 
from  the  electromagnetic  field,  we  multiply  both  sides  of  Eq.(Il-l) 
by  x  and  rewrite  the  result  in  the  form 


d 

dt 


•  ^ 


mx 


+  V  (x)  J  =  E  ( t) 


dM(x) 

dt 


(I  I  -2) 


The  left  hand  side  of  this  equation  is  the  ins-t  antaneous  time  rale 
of  change  of  the  energy  of  the  molecule;  we  will  denote  it  by  d£  4lt  . 
It  is  not  the  instantaneous  time  rate  of  change  of  the  energy  in 
the  molecule  we  require,  but  rather  its  average  with  respect  to  the 
cannonieal  ensemble  described  by  the  Hamiltonian  lor  the  system 
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II 


2 

=  P_ 

2n 


+  V (x)  -  M (x)E (t) 


and  the  time  average 
The  aveiagc  with 
pressed  in  the  form 


(11-3) 


-  Hq  -  M(x) T (t)  , 

of  the  resulting  expression. 

respect  to  the  canonical  ensemble  can  be  ex- 


+  co  +  co 

4£> "  dp  ^  dx  p(x,p,t> 

—  CO  —  CO 

-f  CO  -f  CO 

~  [  dp  ^  dx  p(x,p,t)E(t)  m  (x) 

—  CO  —  CO 


(I  I  -4) 


where  p(x,p,t)  is  the  canonical  distribution  function  which  obeys 
the  Liouvillc  equation 


dp  +  _dp  _dH  dp  5H  _ 
dt  dx  dp  dp  dx 


In  view  of  Eq.(II-3),  this  equation  can  be  rewritten  as 


ae  ^  ao  8Ho  a£  SHo  _  a£  cim(x)  .  . 

"a  t  +  eT  15  ■§£ - Ip  ~a~  E(t)  • 

We  now  expand  p  in  powers  of  the  driving  electric  field  of 
infrared  radiation, 


(II-5) 


(II-6) 


the 


P 


Pq  +  px  +  •••' 


(I  I  ~  7  ) 


where  the  subscript  denotes  the  order  of  the  corresponding  terra  in 
E(t).  When  we  substitute  the  expansion  of  Eq.(II-7)  into  Eq.(II-6), 
and  equate  terms  of  lihe  order  in  E(t)  on  both  sides  of  the  equa¬ 
tion,  we  obtain  the  system  of  equations 
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Sx 

dpx  dp1 

5Ho  . 
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Sp 

Sx 

“  0 


Sp 


°  4?  E(t) 


Sp  dx 


(11-8) 


(TI-9) 


We  now  use.  tne  results  that 


an 

o  p 

Sp  m  x 


(II -10a) 


SH 

"ax 


o  cf 
dx 


A%r  ~  P 


(1 1  -  10b) 


Equation  (I I  —  8 )  and  Eq.(II-9)  can  then  he  rewritten  as 

It  p0  ‘ 0  (u-ii) 

<L  „  «  2^2  dM  T'/i  \  ,TT  -9x 

dt  P1  Sp  dx  L(t)  »  (II-J.2) 


where  dp  /dt  ana  dpn/dt  are  the  total  time  derivatives  of  p  and 
o  r  i  f  o 

p^,  ,  respectively.  From  a  physical  point  of  view,  dp/dt  is  the 
change  in  the  distribution  function  seen  by  an  observer  moving 
with  a  particle  that  traverses  the  orbit  generated  by  the  Hrunilton- 
ian  H,  and  which  passes  through  the  point  (x,p)  m  phase  space  at 
time  t  . 

For  the  equilibrium  distribution  function  pQ  we  assume  the 
canonical  form 


„  cxp(~  e}I0) 
°  Z 


f 


(11-13) 


where  the  partition  function  Z  is  defined  by 

+  CO  +• 

z  “  ^  tip  ^  dx  exp(-'tJiIo)  .  (11-14) 

—  CD  —  CO 
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If  we  now  use  the  fact  that 


dp  dp  oil 
po  Ko  o 


dp  dHo  dp 


“  PP 


J2  = 


o  m 


-  Bp  x 


(11-15) 


Eq. (11-12)  for  p^  becomes 

dpl  .  •  dM 


dt 


.  •  UiW  n  /  .  \ 

Pl'o*  d5  E(t) 


(11-16) 


We  new  assume  that  the  perturbing  electric  field  was  switched  on 
adiabatically  in  the  infinitely  distant  past,  so  that 

p^(-e=)  =  0  (11-17) 

The  solution  of  Eq. (11-16)  which  obeys  the  initial  condition 
Eq.  (11-17)  is 


■>!<«  "  Po  4(t')['dl]x(t,)E<t')dt' 


(11-18) 


With  the  results  given  by  Eqs .  (11-7),  (11-13)  and  (11-14),  we 
can  rewrite  Eq.(II-4)  in  the  form 

<4f>  -  \  d»  5  d*  <=o  -  dx 


i  E(t) 


—  09  —  00 


+  00  +09  t 


+  s  \  dp  V  dx  Uo  *<t'>[!U  ,  ft  En')E(t)dt' 


—  CO  —  00  —  00 


(11-19) 


If  we  use  Eqs.  (II-3)  and  (11-13),  and  rewrite  the  first  term 


on  the  right  hand  side  of  Eq.  (11-19)  in  the  form 

+  «  +re 

Vdxo'6V(x)fi  , 

m  j  dx  ’ 


E  (t)  ^  dp 


(11-20) 
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we  see  that  it  vanishes  cluo  to  the  vanishing  of  the  integral 


dp 

over  p.  Thus,  the  first  non-vanishing  contribution  to  (-^) 


comes  from  the  second  term .  This  terra  may  be  arranged  to  read 

t 

,de 


<4l)  "  P  l  <M(t  /)M(t)  >o  E(t')E(t)dt'  , 


(11-21) 


—  CO 


where  we  have  introduced  the  notation 


M 


(x(t)j  =  M(t)  , 


(13  -  22) 


and  where  for  any  function  A(x,p,t)  , 

4-  co  4  co 


<A)q  -  dp  ^  dx  po(x,p)A(x,p, t) 


(11-23) 


For  E(t),  we  now  assume  the  form 


E(t)  «  Eqcos  at  e  , 


(11-24) 


where  er,t  is  an  adiabatic  switching  factor  (r)  is  a  positive  infin- 


tesmal) .  With  th^s  choice,  Eq. (11-21)  becomes 

t 


=  eE2  \  <fi(l')ft(t)  )  cos  at '  cos  at  er|^t  +  1  ^dt'.  (11-25) 
d  t  o  j  o 


Let  t 1  *-  t  -  T  ,  and  integrate  over  t  rather  than  t*  .  Then 

r 

“  f.E2  e2r|tcos  at  V  <M(t-  v)M(t)  >  cos  u(t-  Oo"“Td’; 
'dl  f  o  j  c’ 


o 


(11-20) 


Because  H  is  time  independent,  our  system  possesses  time  Iransla- 
o 

tion  invariance,  which  in  the  context  o'  1 1 1  -  ■  prese  r.  I  problem  u; 
expressed  by 
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(M(t-  T)M(t)>0  ~  <M(o)M(t))0 


(11-27) 


The  time  averaRe  of  E,.  (11-26)  thus  becomes  (in  the  limit  n  -  «) 

r“  UT  cos(ui7)c'1t(M(o)M(t))0  •  (I1  28) 

'  ci  t  ^ 

.  _c  m  in leer at  ion  variable, 

In  Eq .  (II-2S)  ,  we  can  replace  T  by  - 

and  use  time  reversal  symmetry,  «hich  loads  to  the  identity 

(M(o)M(-  7  )  >c  *  (M(o)»i(7)>0  •  (I1 

d  ^ 

We  then  obtain  an  alternative  expression  for  «aj» : 

iffi  l  d7  cos(o;T)e+’,lT<M(o)H(7)>0  (II'30) 


Upon  adding  half  of  Eq.<H-*»  t0  *‘-(II-3°)*  ”  “ 

^  <il(o)!l(7))cos  0)7  e  ^  'dT  .  (H-31) 

—  C3 

ave raced  rate  at  which  eneigy 
for  the  time  and  thermodynamically  axel  age 

i noiii ii  np tod  on  by  nn  extern^.-  n  c 
is  absorbed  by  a  diatomic  molecule  acted  o) 

electric  field. 

The  result  in  Eq. (11-31)  may  be  recognised  as  a  classical  vei- 

sron  of  the  well  Known  Kobo  fo.mula  for  the  absorption  coefficient. 

viia t  follow* .  we  shall  confine  our  attention  to  the  case  of 

infrared  absorption  by  a  first  order  dipole  moment.  Thai  is.  1* 

VO  expand  M(x)  in  a  Macla.urin  series 

,  .  .  ^  1  .,(/  ,nw?  ...  (I I -d2) 

y,(j:)  =-•  M'  (e)*  +  ~j  N»  to)x  '  » 
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where  the  primes  denote  differential  ion  with  respect  to  x  (we 
assume  the  equilibrium  configuration  of  the  molecule  has  no  dipole 
moment),  we  retain  only  the  contribution  from  the  leading  term  in 
the  expansion.  McGill  et  al . ,  have  examined  the  effect  of  the 
second  term  on  the  absorption  coefficient  and  find  its  effect 
quantitatively  small. The  coefficient  M'(o)  has  the  dimensions 
of  a  charge,  and  we  rs-e  it  by  q  in  what  follows.  Thus,  the 
starting  point  for  the  investigations  in  this  paper  is  the  follow¬ 
ing  expression  for  the  average  rate  of  energy  absorption  by  a 


diatomic  molecule: 


«3f»- 


2  2  +oo 

pq  E  _  -nl-rl 

- T  \  <P(°)P(T>>0cos<^>e  ^  d' 

4m  ^ 


(11-33) 


We  next  turn  to  the  problem  of  casting  the  momentum  autocorre 


lation  function  <p(o)p(t)>0  into  a  form  convenient  for  computa¬ 
tional  purposes.  This  autocorrelation  function  can  be  written 
explicitly  in  the  form 


+  OT 


+  00 


<p(c»)p(t)  >o  =  \  jj  dp  §  dx 


-ih +  v(x)] 


pp  (t)  .  (11-34) 


Because  the  Hamiltonian  is  time  independent,  we  ha\  c  expressed  it 
in  terms  of  the  values  of  p  and  x  at  time  t  -•  o  in  Lq .  (11—34)  . 
Thus,  here  and  in  what  follows  p  and  x  denote  p(o)  and  x(o), 
respectively .  In  addition,  as  integration  of  the  equations  oi 
motion  si  lows ,  the  value  of  the  momentum  at  time  t  is  a  function 
of  the  initial  values  x  and  p,  and  we  indicate  this  explicitly 
by  writing  ]>(t)  as  p(x,p,t). 

We  now  rewrite  Eq .  (1 1-3 '. )  in  the  term 
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+  00 


+  00  +  CO 


<p(o)p(t) >o  - 


1  ^dEe~pE  ^  dp  ^  dx  6(E-  V  (x)  -  -^-)pp_(x,P, t} 

—  00  —  00  —  00 

x9(E) 

"  or  2  t  6(p-pF)  +  6(p+PF) 

-  |  ^  dEe~pE  ^  dx  ^  dp  - - — PP(x,P,  t) 


E  . 
mxn 


x1(E)  -oo 


(11-35) 


where 

p£  -  y2m(E-  V(x))  (11-36) 

and  x  (E)  and  x2(E)  (chosen  so  x^E)  s  x2(E))  are  the  classical 
turning  points  for  motion  in  the  potential  V  (x)  ,  i.e.,  they  are 
the  solutions  of 

E  =  V  (x)  .  (11-37) 


We  assume  the  potential  V (x)  is  such  that  there  are  two  classical 
turning  points  for  energies  E  *  Emin  ,  where  Emin  is  the  minimum 
value  of  V (x) .  The  physical  significance  of  pk  is  that  it  is  the 
momentum  at  t  -  o  in  a  motion  corresponding  to  the  total  energy  E. 

Upon  carrying  out  the  integration  over  p  in  Eq. (11-35),  we  ob¬ 
tain  the  result  that 


<p(o)p(t))0  -  ~  ^  dE  e"pE  H  dx  p(x,pE,t)  + 


x2(E) 


min 


x1(E) 


X1.(E) 


(11-38) 


^  dx  p(x,  -  pE,  t) 


x2(K) 


Since  1  he  momentum  at  t  ~  o  for  x  in  the  interval  (x2 (}.)  ,x^ (E)  ^ 
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is  the  negative  ol  that  at  the  shmc  point  in  l lie  interval 
^x^(E),  x,^(E)^,  because  the  motion  reverses  itself  at  each 
turning  point,  the  expression  in  brackets  is  the  integral  over 
one  period  of  the  motion  beginning  at  x.j(t),  and  returning  to 
x^(E)  after  one  period.  Thus,  we  m;  v  write  Eq. (11-38)  in  the 
form 

CD 

<p(o)p(t))o  «  ^  dl’  e  dx  p  (x ,  p]; ,  t  )  .  (11-39) 

E  . 
mill 


The  one-dimensional  motion  of  a  particle  in  a  region  bounded 
by  two  turning  points  is  a  periodic  function  of  time  with  a  period 
T(E)  given  by ^ 


*2<e> 


T(E)  -  ^ 

xx(E) 


dx 

Ve-  V  (x) 


(11-40) 


This  result  holds  for  any  initial  position  x  and  momentum  p^  in  a 

motion  corresponding  to  total  energy  E.  Thus,  the  integral  over 

a  period  in  Eq. (11-39)  is  a  periodic  function  of  time  with  the 

same  period  T(E),  and  wc  expand  it  in  a  Fourier  scries: 

+» 

(j^  dx  p(x,pE,t)  £  Pn(E)p  1,lU)(EK  (11-41) 

n  =  -ra 


where 


T  (E) 


Pn(E)  "  t~(K)  \  dt  dx  P(x»EE,t)o' 


*  i  r.  a  (E)  t 


(1  j  -42) 


and  we  have  defined 


a'(E) 


:vr 


(11-43) 


To  obtain  a  simple  expression  for  p  (E) ,  note  that  the 


Si 


soli,.!,,.,  ,.|  the  equation  of  .notion  for  a  panic],.  moving  in  a 
one  dimensional  potential  it,  given  explicitly  by<5) 


/-AS  n  C  dx  ' 

'2  J  IFr^  ’  * 

Xj(E)  '!L  V<x  ) 


(11-44 


~ U  Vlit'illJ>  tllis  ession,  v/c  aro  measuring  time  with  respect  tc 
an  instant  tQ  at  which  the  particle  is  at  the  loft  hand  turning 
point  x;|(E).  It  is  necessary  to  know  x  as  a  function  of  t  only 
foi  (E)  -  x  s»  Xg (I.)  because  as  t  increases  past  t  +T(E)/2,  where 
x  -  x2(F),  the  motion  reverses  itself  (i.e.,  x(t)  is  symmetric 
a !)o lit  t  «  tQ  +  T(F)/2)  until  the  time  tQ  +  T (E)  is  reached,  at 

which  point  the  particle  has  returned  to  x^E),  and  the  motion 
begins  to  repeat  again. 


The  solution  of  Eq . (11-44)  can  be  written 
x(t)  -  x1(E)  +  fE(t  -  to) 


(11-41 


where  the  function  f£(t)  is  an  even  function  of  t,  is  periodic 
m  t  with  period  T(E),  is  even  about  t  -  T(E)/2,  and  vanishes  as 
t  -  o.  The  momentum  Pf,(t)  for  the  orbit  of  energy  E  is 


dl  (t-  t  ) 

PE(t)  -  m  - - »  m  gE(t-  tQ)  , 


(II-4 


where  gj.U)  is  an  odd  function  of  t,  is  periodic  in  t  with  period 
T(L),  and  is  odd  about  t  =  T(E)/2 

With  these  results ,  Eq. (11-42)  becomes 

T(E) 

Pn(l)  YfL)  ^  clxc  ^  ^  L°  ^  m  gI'(t"  t0)cinU)^L)  °"to) 

(11-47) 
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whore  the  integrations  over  x  and  t  have  boon  interchanged. 
Since  g£(t)  is  periodic  with  period  T(E),  and  wo  integrate  over 
a  complete  period,  this  result  becomes 


i. n li; ( E )  t  T(E) 

Pn  “  YXE)  &  d*  e  °  S  dl  m  eE(t)sin(nu)(E)  t)  , 


(11-48) 


where  we  have  used  the  fact  that  g£(t)  is  odd  about  t  =  T(E)/2  . 
We  now  convert  the  first  integral  from  an  integral  over  x  into  an 
integral  over  t  : 

-T(E)  T (E) 

.  r-  anu;(E)t  .  . 


pn(E)  "  tTe)  $1  dto 


^  dt  m  gK(t)sin(nu(E) t)  . 


(11-49) 


The  coordinate  x  is  obtained  as  a  function  of  t  bv  setting  t  =  o 

in  Eq. (11-44).  The  limits  on  the  first  integral  follow  from  the 

fact  that  the  original  integral  on  x  around  a  period  of  the  motion 

corresponds  to  t  -  t  increasing  from  0  to  T(E)  as  x  goes  from 

x-^(K)  to  Xg^),  and  back  to  x^(E).  Setting  t  =  0,  we  see  that  t 

goes  from  o  to  -T(E)  as  x  makes  the  same  circuit.  Making  the 

change  of  variable  t  =  -t',  and  using  the  fact  that  f„(t)  in 

o  i. 

Eq. (11-45)  is  an  even  function  of  t,  we  obtain  for  pfi(E) 

T(K)  ^  T(E) 

pn (E)  “  TCeJ  S  ‘H'  gE(t/)e';,Ua(E)t  ^  dt  m  gE(t)sin(nuj(E)t) 
o  o 

(11-50) 

or,  using  the  oddness  of  eE(t. ')  about  t'  i=  T(E)/2,  we  find  fiv.ally 


pn(E)  "  tTE)  pn  <“>• 


(11-51) 
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where 


T  (E) 

Pn(E)  “  ^  cl t  pE(t)sin  nu,(E)t  , 
o 


(11-52) 


where  the  momentum  pr(t)  "  n  gr(t)  appearing  in  Eq. (11-52)  must  be 

Ij  L 

understood  to  be  obtained  from  Eq. (11-44)  with  t  =  o  (since  that 
is  how  g (t)  is  defined) . 

Upon  combining  Eqs .  (11-39),  (11-41),  (11-43)  and  (11-51), 
we  obtain  for  the  momentum  auto-correlation  function 


+  oo 


<P(o)p(t)>  “ 


1  s 


dEe 


-PE 


2nZ 


0)(E)Pn2  (E)c“inUj(E)L 


(11-53) 


n~  -  t» 


E 


min 


When  this  result  is  substituted  into  Eq. (11-33),  and  the  integra¬ 
tion  over  tine  is  carried  out,  we  obtain  (as  n  -♦  o  +  ) 


WV  +” 


«H»"— 27  I  S  «Ee-^(E)P>x 


n=-  »  E 


min 


x  nuu(E)  ^  +  6^uj  +  nw(E)  . 


(11-54) 

2 

If  we  note  that  (E)  vanishes  for  n  -  o  ,  and  is  an  even 
function  of  n  for  n  f-  o,  wo  finally  obtain  for  the  average  rate  of 
energy  absolution  by  a  diatomic  molecule 


_  2r,  2 
a>pcj  Eo 

4m2Z 


c&  CO 


US 

ii“l  E 


dEe 


-PE02- 


nun 


(11-55) 


where  we  have  assumed  0;.  >  o  . 

The  dynamics;  ol  the  problem  are  seen  to  enter  simply  through 
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the  necessity  of  knowing  pr(t)  for  the  evaluation  of  9  (E) ,  and 

s  n 

uu(E)  .  Moreover,  if  uu(E)  is  a  sufficiently  simple  function  of  E 
so  that  the  equation  uu  =  nuu(E)  can  be  inverted,  the  integration 
over  E  can  be  carried  out  using  properties  of  the  delta  function. 

We  briefly  summarize  in  words  the  procedure  for  computing 
the  quantities  which  enter  Eq. (11-55).  Given  a  potential  function 
V (x)  ,  one  requires  the  period  T(E)  as  a  function  of  energy.  This 
function  may  be  obtained  from  Eq.(TI-40),  and  w(E)  is  defined  by 
Eq. (11-43).  By  solving  the  equations  of  motion,  one  finds  the 
momentum  as  a  function  of  time  p^(t)  for  an  orbit  of  energy  E, 
with  the  origin  of  time  chosen  so  the  parameter  t  in  Eq. (11-44) 
is  set  equal  to  zero.  The  quant  ty  f>n(K)  is  a  measure  of  the 
amplitude  of  the  n^1  harmonic  in  the  function  p„(t),  and  is  ob- 
tained  from  Eq. (11-52).  Finally,  7.  is  the  partition  function. 

We  conclude  the  present  section  by  displaying  a  remarkably 
simple  relation  between  the  partition  function  of  the  oscillator, 
and  the  function  T(E),  that  gives  the  period  as  a  function  of 
energy  . 


Wc  have 
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or  finally 

03 

Z  =  ^  dE  e“pE  X(E)  . 

E  . 
mi  a 

Wo  next  proceed  to  apply  the  results  of  the  present 
to  the  examples  mentioned  in  section  I . 


(11-56) 


section 
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III.  Applications  of  the  Formal ism  to  the  Study  of  Mu3  tiphonon 
Absorption  for  Some  Specific  Potentials 

In  this  section,  we  study  the  behavior  o'1  the  absorption  co¬ 
efficient  as  a  function  of  frequency  foi-  four  specific  forms  of  the 
interatomic  potential  V(x)  .  We  first  derive  the  form  of  the 
absorption  coefficient  for  the  Morse  potential,  a  form  used  fre¬ 
quently  in  molecular  physics.  We  then  consider  the  potential 
o  o 

V(x)  =  ax  +  b/x  ,  the  square  well,  and  the  potential  V(x)  “ 
yx  for  x  s  0,  V(x)=®  for  x  <  0.  The  last  two  potentials  are 
interesting  examples  to  consider,  since  one  cannot  construct  a 
discussion  of  a  perturbative  nature,  because  an  harmonic  approxi¬ 
mation  does  not  exist  for  either  case, 
a)  The  Morse  Potential 

The  Morse  potential  may  be  written  in  the  form 

V(x)  «=  D[  1  -  exp(-a[x  -  xq])]2  .  (III-l) 

The  minimum  value  of  V(x)  is  zero  and  the  minimum  occurs  at  x"xo< 
For  large  values  of  the  interparticle  separation  V(x)  approaches 
the  constant  value  D,  the  dissociation  energy  of  the  molecule. 

For  most  cases  of  interest  here,  D  assume  a  value  the  order  of 
one  electron  volt,  an  energy  very  large  compared  to  k^T,  as  long 
as  we  confine  our  attention  to  temperatures  the  order  of  or  lower 
than  the  melting  temperature  of  the  solids  of  interest  to  us.  We 
shall  make  use  of  the  fact  that  KRT<<D  i r*  the  discussion  below. 

If  we  consider  only  motions  of  small  amplitude,  x  remains 
near  x  ,  and  V(x)  is  well  approximated  by  t  lie  parabolic:  form 

V(x)  -  Da2 ( x  -  xq) “  .  (111-2) 
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Thus,  in  the  limit  where  the  amplitude  of  the  motion  is  small,  the 
molecule  behaves  like  an  harmonic  oscillator  with  frequency  oi0 


given  by 


2a2D 

m 


(III-3) 


We  begin  by  deriving  an  expression  for  the  quantity  Pj/t) 
defined  in  section  II,  following  the  procedure  outlined  there. 
We  begin  with  Eq.  (11-44),  which  with  t^  equal  zero  becomes 
t  f)1/  *"  a(x~“x)  2a  (x„— '1  "2 

x1(E) 


dx’ 


^E-D 


a(x  -x) 

h  2De  0  -Do 


X  -x) 
o 


(III-4a) 


The  integral  can  be  evaluated  in  closed  form.  This  may  be 

done  by  letting  8  =  exp(-ax)  ,  and  01  *=  cxp(-ax^(E) )  .  If  we  then 
define  C  =  D  cxp(2axo) ,  B  -  2D  exp(ax  )  and  A  =  D  -  E,  the 


integral  becomes 


1  -  id//1  [B9'  -A 

e 


(III-4b) 


For  bound  motions  of  the  molecule,  the  only  case  cf  interest  here, 

2  /  g\ 

B0  ;?  A  +  C0  everywhere.  The  integral  is  then  readily  evaluatodv  ' 


to  give 


at  [j^(D-E)  ]?:  -  sin  1 


ax 

(D-E)-De  °  0 


ax 


./DE  e  V 


J  0- 


-1 

sin 


ax 


(D-E)-Do  °6 
ax 

/DE  c  °e 


TT 


( 1 1 1-5) 
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If  thiL  relation  is  solved  for  c  a  function  of  t,  and  x(t) 
is  obtained  front  this  result,  one  finds 


x(t)  = 


xo^ 


1 

a 


tn 


I) 

D-E  + 


1  ,  n 

-  til[  1  -  (jr) 

a  D 


cos  i  (F.)  tj 


(II 1-6) 


where  wo  have  intro Juced  the  Quantity  a (E)  ,  deiined  by 

i 

*(E)  -  x  [1  -  £l2  .  (II1-7) 

o  I)  • 

From  Eq.  (II1-6),  it  is  evident  that  the  period  T(E)  of  the 
motion  is 

'  <IU‘8) 

For  small  energies,  T(E)  assumes  a  value  independent  of  energy 
and  equal  to  2 tt/j  ,  and  the  period  lengthens  as  E  increases. 

The  momentuu  PE(t)  is  found  from  r'.x(t),  with  x(  t)  given  by 
Eq.  (Ill -6) .  If  this  differentiation  is  carried  out  and  the 
result  substituted  into  the  expression  for  f>n(E)  given  in  Eq .  (11-52), 
one  obtains 

j  T(E) 

(E)  (2mE) 2  i  (E)  r  sin  nr (E)  t  siin  (E)t  dt 

Pii  '  J  y  2 

°  O  1  -  (^)  COSJL'(E)  t 


2(2mE)  2 
w 


rr 


/ 


sm  n  ^  s  i  n  o 

i  i  2 

1  ~  tp)  cos;-, 


d? 


(II 1-9) 


The  integral  in  Eq .  (II-9)  may  be  evaluated  exactly,  to  give 

l  (n  'JL)  l 

<i>  ~2  o-u-ji,2]  . 


(6) 


( 1 11-10) 
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The  discussion  shows  that  for  the  Morse  potential,  exact 
results  for  all  quantities  which  enter  the  expression  for  the 
absorption  coefficient  are  readily  obtained,  for  the  bound  orbits. 
We  shall  compute  an  approximate  form  for  the  absorption  rate  in 
the  linn  t  i^T«l)f  where  only  orbits  with  energy  E«D  contribute 
significantly  to  the  rate  of  absorption.  In  this  regime,  we 
leplace  ^(F)  by  its  leading,  contribution  when  E«D: 


p  (E)  =  — -  ( 2mJ))  "  (  E  \ 

n  %  \~w  I 


n 

2 


(111-11) 


In  this  limit,  the  function  u,(E>  may  also  be  replaced  by  the 
approximate  form  wo  independent  of  energy,  for  E«D.  From  Eq. 

(II-jG),  one  readily  sees  that  the  expression  for  the  partition 
function  becomes 


Z  *= 

V 

If  these  approximations  ar 


(HI-12) 


e  inserted  into  the  expression  for 


,de 


«dt»  displayed  in  Eq.  (11-55),  one  finds  that 


where 


^dt^  ^ 

n=l  n  ° 


q2K2 


%  “  " 


n  **  “o  ,  /kBT\  n_1 
4  m  n*  (  4D  ) 


Of, 


( III-13a) 


( II I-13b) 


For  the  case  n  -  1 ,  we  have 
tt  q2l‘o 

]  4  “ Jfj  >  independent  of  temperature, 


( 1 1 1  ~l3c  ) 


126 


tlu*  well  known  result  for  the  integrated  strength  of  the  funda¬ 
mental  absorption  line  of  a  simple  harmonic  oscillator  of  mass  m, 
charge  q  and  frequency  jlq.  (Recall  that  E  is  the  peak  value  of 
the  field.) 

The  integrated  strength  of  the  absorption  peak  at  the  fre¬ 
quency  a)  =  njj  (the  n  phonon  absorption  peak  for  this  model)  is 
° 

related  to  that  of  the  fundamental  absorption  peak  by  the  simple 
relation 


O'. 


n 


c*i 


n  ! 


n-1 


(II I- 14) 


b)  The  Potentia 1  V(x)  =  (a/x2)  +  bx2 

The  second  potential  function  we  consider  is 

V(x)  =  -^  +  bx2  x  2t  0.  ( 1 1 1  -15) 

Setting  tQ  =  0  in  Eq.  (11-44),  we  obtain  for  the  equation  deter¬ 
mining  x(t) 

x 

t  =  §  f  -  -xdx  =.  .  ( 1 11-16) 

6  J  / o 4 

Xi(E)  W-a  H-  Ex  -  bx 


The  turning  points  x.^  are  the  solutions  of  the  equation 

E  =  ~  +  bx2, 

X 


( 1 1 1 -17) 


and  are  found  to  be 


Xj_(E) 

x2(U) 


L 

[ 


E  -  [  E2 
2b 


E  +  [  K2  ~  4  a 
2b 


(III-lSa) 

(Ill-lSb) 
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The  energy  Kmin  equals  the  minimum  value  of  the  potential 
enei  CV  (HI-15).  Equivalently,  it  is  the  energy  for  which 
x,  (K)  *=■  x0(L)  ,  and  hence  is  given  by 

E  .  -  2,/aK  .  (III-19) 


The  change  of  variable  x2  =  y  transforms  Eq.  CHI -16)  into 


It  should  be  noted  that  the  period  T(E)  and  the  corresponding 
frequency  «(E)  for  the  potential  (III-15)  are  independent,  of 
the  energy  E. 

It  follows  immediately  from  Eqs.  (11-56)  and  (III-21)  that 
the  partition  function  Z  is  given  by 


z  =  ln_  e-2^ 


(111-23) 
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Ihc1  expression  (11-55)  for  the  rate  at  which  energy  is  absorbed 
by  the  molecule  described  by.  the  potential  function  (111-15)  now 
becomes 


«5T» 


U.JU 


8*7 


,2 


m 


Eo  J,  n  j  e2^  /  }  6<»-  n»o: 

E  . 
mm 


(111-24) 


If  we  make  the  change  of  variable  E  =  2/ab(z +  1),  we  obtain  the 
convenient  form 

<(at»  =  ~&r  P2  \  Eo  ^  n  j  /  dze~2d'/*15*  p2(2y^F(z+l)  )  j 

"7  n=l  (  ^  “  I 

o 


x  6(uj  -  nxQ) 


(II 1-25) 


Sxncc  2v/ab  -  Emin,  we  see  that  for  temperatures  such  that 

Lmin  ~  kBT’  only  volues  of  z  ^  1  contribute  signif icantlv  to  the 
integral.  fte  now  turn  to  a  determination  of  iJ  (F.)  . 


Inverting  Eq.  (II 1-20)  we  obtain  for  x(t) 


x(t)  -  - 


i  E  VP2  -  4; 


4a  b 


I  2b  cos  ,J,ot  \ 

from  which  we  immediately  obtain  p(t): 


i  i 


P(t)  =  ^ 


V F2  -  4a b 


s.in  x  t 
o 


]j 


cos  I  t 

o 


(II I -2d) 


(IIJ-27) 


The  integral  i n  this  case  is  given  by 
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,  „  /— -  6  f  sm  x  y 

3n(k)  “  2,v/mL  /  - 

11  xo  •L  yi  -  6cosx 


sin  nx 


(III  -28) 


where  v,e  lmve  set 


f  o 

Yr  -  4ab 


( I 11-29) 


The  change  of  variable  x  —  tt  —  2cp  yields  the  following  expression 


for  'Pn(E)  , 


P„(E)  -  (-1) 


n-1  ,/roE  26  C  cos  2(n-l)o-cos2(n+l)cp 


iK  26  C  cos  2 

’°  yus  J  XT 


26  2 
U6  Sin  * 


(I 11-30) 


The  integral 


£* 

Kn(k)  ’  / 


cos  2nx 


o  Vl  ~  k2sin2x 


has  the  following  expansion  in  powers  of  k 


0  <  k  <  1 
n  -  0,1,2,. 

2  (7) 


( 1 1 1-31) 


Kn(k)  ,  |  E  ^ 


(v-n)  !  (v+n)  I 


,  -t  Itn •],)  k2"  \  (2ml)  k2 

-  (-1)  -2  ]  T2T,)T|1+—  4  k 


( I II-32a) 


(2ii+1)(2ih3)2  ,.4  ,  I 

- K  +  •  •  •  i 

32 (2n ! 2)  ' 


(III  -32b) 


Consequently,  since  in  the  present  case 


2  ° ; 

k  t-t—  <  1  • 


(111-83) 
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the  first  two  terms  in  the  expansion  of  f>n(E)  are  found  to  bo 


o  (E) 

n  2  % 


retail 

1 

_ (_2j0_ 

i'W 

(2n-2TT 

(li6)n 

(2n-l) 

26  ' 
1+6  +  °V 

(  26  \2\ 

4 

[1+6/  ) 

V 


Consequently,  v/c  have  that 


n2,r>  rr  inE 
Pn(E)  ~  T  ~2 
“o 


Hn-il 

r(£) 


[(2n-2)!]2  (1+6)2n-l 


( III-34) 


j  i  .  ( 2n-l)  26  „//  26  \2  \  ) 

|1+~ -  1+3  +  °(i(t+t)  )  •  (111-35) 


Since  6  <1,  we  can  simplify  this  expression  to 


<V(E> 


2  mE  [  (2n-2)  !  1  ,2nr,  _  ,  2.  , 

n  “2  *~Gn-6 - ~~4  6  +  0(6  )]  . 

2bn  b[(n-l)I]4 


(1 11-36) 


If  v/e  make  the  replacement  E  =  2  %/aF  (z  +  1)  in  Ei .  (1 1 1-29)  , 
we  find  that  as  a  function  of  z  5  is  given  by 

,  _  ,/z2  +  27, 

1  7  +  1 

-  Jto,  (1  ••  |  z  h  %?+...)  o  s  z  <  1  (IIT-37) 

Thus,  the  leading  term  in  the  expansion  of  (2  N/aF(l  A  z)')  for 
small  z  is 
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F  2  (2 (l+z))  ~  n2  -Eg-  2^ir-T^1-~g>-’32  Zn  +  0(Zn+1)  . 
V  '  ui„  2',n  6[(n-l)  !]4 


(II 1-38) 


When  this  result  is  substituted  into  Kq .  (III-25)  ,  and  the  integral 
over  z  is  carried  out,  we  obtain  finally 


«af»" 


_n  (j2  2r  C(2n-2);12n!  /  kBT 


L  n5n-5r,  ,  N  ,  -,4  \  K 

n=1  2  C(n-l) !  ]  \  mm 


6  (urnoio) 


(III  -39) 


We  see  from  this  result  that  the  integrated  strength  of  absorption 
at  the  frequency  uu  «  gu^,  (nq  /4m) Eq,  is  independent  of  temperature, 
and  has  the  same  value  as  in  the  case  of  the  Morse  potential.  The 
integrated  strength  of  the  absorption  at  the  frequency  u)  =  nuu  is 
related  to  that  of  the  fundamental  absorption  peak  (n  =  1)  by 


[ (2n  -  2) l]2 
!5n"5[(n-l)*.)4 


/kRT 

nI  (e~“ 

'  min 


(II 1-40) 


which  is  very  close  in  form  to  the  corresponding  relation,  (III-14) , 
obtained  for  the  Morse  potential, 
c)  The  Square  Well 

We  next  apply  the  method  to  compute  the  shape  of  the 
absorption  spectrum  for  a  particle  trapped  in  a  square  well  poten¬ 
tial,  with  inJinitely  steep  sides.  This  example  is  an  interesting 
application  of  the  formalism  developed  in  section  II  of  the  present 
paper,  since  the  potential  does  not  admit  an  harmonic  approxima¬ 
tion.  In  section  IV,  wc  shall  argue  that  there  are  certain  impuri¬ 
ties  in  a?. Pali  halides  that  may  bo  regarded  as  moving  in  a  very 
steep  sided  potential  that  can  he  approximated  by  a  square  well, 
at  least  for  qualitative  purpose's. 


We  suppose  the  square  well  lias  width  L,  with  infinitely  steep 
sides.  If  we  consider  motion  of  a  particle  of  mass  m  in  the  well, 
with  snergy  L\  then  the  velocity  v  of  the  particle  is  (21J/m)  '  . 

The  period  T(E)  is  2L/v,  or 


T(E)  -  L 


<f,5 


2n 

TTeT 


For  the  momentum  p  (t.)  ,  one  has 


(111-41) 


A 


short 


pe(d  - 

calculation 


+  ( 2m  E)  ^ 

-  (2inE)2 
gives  for  this 


0  <  t  <  T(E)/2 
T(E)/2  <  t  <  T(E) 
potential 


(III-42) 


Pn(E)  “  TT  [1  " 


It  is  a  straightforward  matter  to  insert  this  expression 
into  Eq.  ( 11-55)  ,  and  obtain  the  form  of  <<^»  .  The  computation 
of  the  partition  function  is  also  quite  elementary.  The  final 
result  is  best  expressed  in  terms  of  a  characteristic  frequency 
u>T,  given  by 

(IJI-43) 


In  terms  of  u»T,  the  expression  for  the  rate  of  energy  absorption 
is 


Q2p2  g 

/«>  -  Lfo  ajf. 

' dt ;  3/2  m  3 

it  aT  n=o 


exp 


v 

iU 


2 


(2ml) 
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The  absorption  coefficient  is  thus  a  superposition  of  a 
sequence  of  Gaussians.  The  function  on  the  right  hand  side  of 
Eq.  (Ill— 14)  has  a  preunent  peak  very  near  .r  «  xT,  with  weak 
subsidiary  maxima  at  higher  frequencies.  01  particular  interest 
is  the  behavior  of  the  absorption  coefficient  for  frequencies 
larged  compared  to  xT-  I).  this  region,  one  may  find  the  asymptotic 
form  o l  the  absorption  coefficient  by  replacing  the  sum  in  Eq.(III-44) 
by  an  integration.  For  yj  »  ,  by  this  means  one  finds 


«ST»  “ 


4 

TT 


2  .,2 
q  E 
_ _ o 
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U’ 
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UJ_ 
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UL’-j, 
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3/2 
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2_2 

q  Eq  xT 
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( 1 1 1-45) 


Thus,  for  frequencies  large  compared  to  the  characteristic 

frequency  x-pi  the  square  well  potential  gives  rise  to  an  absorption 

-2 

coefficient  which  falls  off  as  uu  ,  a  result  qualitatively  different 
from  the  empirical  form  displayed  in  Eq.  (1-1). 
d)  The  'J  riangulor  IVc  1 1 

IVe  next  display  the  form  of  the  absorption  coefficient 
for  a  potential  V(x)  of  the  form 

/  ®  x  <  o 

V(x)  »  (111-46; 

(  Yx  x  >  o  , 
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where  the  constant  y  is  presumed  positive.  This  is  a  second 
example  of  a  potential  for  which  a  harmonic  approximation  fails 
to  exist. 

The  classical  equations  of  motion  for  this  example  are  quite 
elementary.  For  an  orbit  with  energy  E,  the  period  T(E)  is 


T(E)  =  2  (2mE)  ^ 


( III-47) 


and  the  momentum  of  the  particle  as  a  function  of  time  is 


P 


E 


(t) 


o  <  t  < 


T(E) 

2 


T(E) 

2 

t  <  T(E) 


A  short  and  straightforward  calculation  then  gives 


p  (E)  .  IS* 

n  nyn  ’ 


(III-48) 


and  for  the  partition  function 

Z  *=  (k^T)372 


(III-49) 


When  these  results  are  inserted  in  Eq.  (11-55),  one  finds 
that  tlie  results  assume  the  form 


«at» " 


n 


5/2  m 


qV  ”  3  2  fiT 

q  S  n  exp  (-ii  ) 

oj  n=l  uu 


( I 11-50) 


whore  is  a  characteristic  frequency  given  by 


135 


Sec.  E 


fiT  -  rrY(2mkBT)"2  (111-51) 

When  10  »  n,p,  the  form  of  the  absorption  coefficient  may  be 
deduced  by  replacing  the  sum  over  n  by  an  integration,  as  we  did 
in  the  preceding  example.  In  the  high  frequency  regime,  one  finds 


that 


''dt '' 


m 


(II 1-52) 


a  result  remarkably  similar  to  that  obtained  for  the  square  well. 

In  both  examples,  for  uu  large  compared  to  a  characteristic  frequency, 

_2 

the  absorption  falls  off  as  ju 


IV  .  General  Discussion 

The  purpose  of  this  section  is  to  examine  some  implications 
of  the  results  in  section  III. 

We  first  examine  the  question  of  the  validity  of  the  phenom¬ 
enological  form  displayed  in  Eq.(I-l),  which  is  suggested  by  the 
absorption  data  reported  to  date.  If  this  form  is  in  fact  a 
general  result  which  holds  in  the  limit  of  high  frequencies,  then 

we  should  expect  it  to  emerge  from  our  analysis.  If  we  consider 

2  2 

the  Morse  potential,  and  also  the  potential  bx  +  a/x  ,  our 
independent  oscillator  model  predicts  a  series  of  absorption 
peaks  at  the  frequencies  -  naj^  ,  where  u)Q  is  the  fundamental 

vibration  frequency  of  the  anharmonie  oscillator.  If  the  exponen 
tial  law  is  obeyed  for  these  models,  then  we  should  expect  the 
integrated  strength  of  the  n  phonon  peak  (the  absorption  peak  at 
nuu^)  to  vary  with  n  as  gn  ,  where  £  is  a  parameter  independent 
of  n  that  depends  on  the  details  of  the  intcratonic  potential 
and  the  temperature. 

If  we  examine  the  results  obtained  in  section  III,  we  see 
that  a  relation  of  this  form  does  not  hold  for  any  of  the  four 
potentials  we  have  examined.  The  two  potentials  most  directly 

applicable  to  real  physical  systems  are  the  Morse  potential  and 

9  2  -  . 

the  potential  bx~  h-  n/x  ,  since  both  of  these  potentials  aomit 

an  barmen  ie  approximation.  It  vc  demote  the  integrated  strength 

of  the  n- phonon  peak  by  ct,  .  then  for  the  Morse  potential  v.e  find 
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where  is  the  integrated  strength  of  the  fundamental  reststrahl 

band  at  u'  .  A  very  similar  relation  obtains  for  the  potential 
2  2 

bx  +  a/x  ,  as  we  have  seen.  Because  of  the  factor  of  n!  that 
appears  in  Eq.(IV-l),  a  plot  of  c*n  vs .  n  does  not  give  an  exponen¬ 
tial  lav  of  the  form  displayed  in  Eq.(I-l).  It  is  for  this  reason 
that  we  feel  that  the  proof  offered  by  McGill  et  al . ,  leads  to  con¬ 
clusions  that  arc  not  correct.  Both  examples  we  have  investigated 
lead  to  results  which  contradict  this  conclusion.  As  we  remarked 
in  section  J,  these  authors  have  also  calculated  the  intensity  of 
the  multiphonon  absorption  peaks  for  a  set  of  independent  oscilla¬ 
tors,  each  ot  which  is  described  by  the  Morse  potential.  Their 
calculation  is  carried  out  by  a  quantum  mechanical  means,  and 

yields  a  result  in  agreement  with  ours  in  the  limit  ftiw  <  k„T  , 

o  B  ’ 

where  the  correspondence  principle  dictates  that  the  classical  and 
quantum  mechanical  results  must  concur. 

One  must  then  inquire  into  the  reason  why  the  data  are  so 
well  fitted  by  the  exponential  form  displayed  in  Eq.(I-l).  Of 
course,  it  may  be  that  our  independent  oscillator  model  is  so  highly 
simplified  that  conclusions  based  on  it  are  unreliable.  We  feel  that 
that  the  problem  does  not  lie  here,  but  in  the  fact  that  the  data 
obtained  to  date  extend  only  over  a  rather  small  range  of  frequency, 
from  2xq  to  6.»  at  best.  The  frequency  only  varies  by  a  factor 
oi  at  most  three  through  this  range.  While  the  absorption  coeffi¬ 
cient  changes  by  many  decades  as  one  passes  through  this  frequency 
range,  measurements  over  a  wider  range  of  frequencies,  or  at 
higher  temperatures  will  he  required  to  del ermine  whether  the 
phenomenological  form  in  Eq.(l-l)  is  valid. 
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We  illustrate  this  point  in  Figure  (1) ,  where  we  present  the 
multiphonon  absorption  data  obtained  by  Deutsch^ ,  and  compare  it 
with  the  prediction  of  Eq.(Ivr-l)  at  integral  values  of  nujQ  .  The 
solid  line  in  Figure  (1)  is  a  straight  line,  chosen  with  slope 
such  that  it  passes  through  the  data  obtained  by  Deutsch.  Wo 
have  omitted  the  data  points  simply  to  avoid  cluttering  the  figure 
but,  as  noted  by  Deutsch,  the  data  fall  on  the  straight  line.  It 
will  be  important  for  our  purposes  to  note  that  the  data  points 
all  lie  below  the  frequency  of  750  cm  \  to  the  left  of  the  square 
bracket  which  has  been  placed  on  the  straight  line.  Thus,  the  por¬ 
tion  of  the  straight  line  to  the  right  of  this  bracket  represents 
an  extrapolation  of  the  data  to  higher  frequencies,  in  particular 
to  the  frequency  of  the  C02  laser.  This  extrapolation  assumes  the 
phenomenological  form  given  in  Eq.(I-l)  is  valid  for  all  frequen¬ 
cies,  since  it  is  a  linear  extrapolation  on  a  semi  log  plot. 

In  the  figure,  the  barred  circles  represent  the  prediction  of 
Eq.(IV-l).  To  place  the  points  on  the  figure,  we  have  adjusted  ^ 
and  D  so  that  Eq.(IV-l)  fits  the  data  at  the  two  points  n  “  2  and 
n  "  5  .  When  we  then  calculate  the  strength  of  the  absorption  at 
n  =  3  and  n  -  4,  we  obtain  results  that  agree  very  well  with  the 
data,  to  within  the  accuracy  of  the  graph.  Thus,  we  can  also  fit 
the  function  in  Eq.(JV-l)  to  the  data,  and  it.  is  clear  that  the 
measurements  do  not  extend  over  a  range  oi  frequencies  large  enough 
to  warrant  the  conclusion  that  the  exponential  lav  provides  a 
unique  description  of  the  results. 

It  is  not  hard  to  sec  why  wo  obtain  a  good  Jit  to  the  data.  The 
ratio  k  T/-II.)  is  quite  small,  of  the  order  of  10  2  typical]},  as  we 
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Figure  Caption 

Figure  (1)  Comparison  between  the  multiphonon  data  in  KC£  ,  and 
the  non-interacting  oscillator  model,  for  the  case  where  the 
motion  of  the  oscillator  is  governed  by  the  Morse  potential.  The 
solid  line  passes  through  the  room  temperature  data  points  of 
Deutsch,  and  all  the  data  lie  below  the  square  bracket  placed  on 
the  curve  near  /50  cm  ^  .  The  barred  circles  are  computed  from 
the  theoretical  model,  with  parameters  adjusted  to  fit  the  data 
at  11  =  2  and  n  =  5.  The  crosses  give  the  theoretical  prediction 
when  T  -  600  K,  and  the  dashed  line  is  a  straight  line  drawn  as 
an  aid  to  the  eye. 
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shall  sec  shortly.  II  we  plot  the  data  on  a  semi-log  plot,  then 
we  are  concerned  with  the  behavior  of  lop  ^  as  a  function  of  n. 
Eq.(lV-l)  then  gives 


i  r  41)  1  “1  —  4D(y**  *n 

ln  ""  '■  ■  "t1"  V  -  n  1"  »0  -  in  tjfi]  ■ 

If  n  is  large  enough  for  Stirling’s  approximation  to  be  used, 
log  n’.  n  log  n  -  n,  and 

ln  “n  “  -  "t1  +  111  S  -  1"  »]  ♦  m  [~£]  . 

13  -1  LixB*  J 

Since  4D/kBT  >>  1  ,  the  ln  (41)/k^,T)  contribution  to  the  quantity  in 
square  brackets  is  quite  large.  Furthermore,  if  we  plot  log  &  as  a 
function  of  n,  deviations  from  a  straight  line  arise  only  because 
of  the  log  n  term.  Since  log  n  is  a  slowly  varying  function  of  n, 
a  plot  of  In  an  vs.  n  can  give  a  result  approximated  very  well 
by  a  straight  line,  if  only  a  small  range  of  n  is  examined. 

II  we  accept  the  results  of  our  independent  oscillator  model 
as  realistic  for  the  moment,  then  from  Figure  (1)  one  can  see  that 
extrapolation  of  the  exponential  law  of  Eq.(I-l)  to  higher  frequen¬ 
ce  can  lead  to  serious  errors  in  estimates  of  the  absorption 
coefficient.  For  example,  by  the  time  n  =  7,  the  expression  in 
Eq.(IY-l)  gives  a  value  of  tho  absorption  coefficient  larger  by  a 
factor  of  5  than  that  obtained  from  extrapolation  ol  the  exponen¬ 
tial  lav.  This  example  suggests  that  to  estimate  the  absorption 
coefficient  at  10. 6„  by  the  use  of  Eq.(I-l)  to  extrapolate  data 
Irom  lower  frequencies  may  lead  to  a  serious  underestimate  of  the 
ab.soi  pt  ion  coefficient  at  10. Gp. 
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If  the  temperature  of  the  crystal  is  raised,  the  discussion 
above  suggests  that  deviations  from  the  near  exponential  behavior 
should  be  expected  to  be  more  severe,  and  to  set  in  at  lower 
values  of  n.  To  illustrate  this  point,  on  Figure  (])  we  have 
placed  a  series  of  crosses  to  represent  the  prediction  of  Eq.(IV-l) 
at  T  =  600CK,  once  and  D  have  been  adjusted  to  lit  the  room 
temperature  data  at  n  ~  2  and  n  ~  5,  as  described  above.  The 
dashed  line  is  a  straight  line  place  on  the  graph  as  an  aid  to 
the  eye.  The  deviations  ax'e  indeed  more  pronounced,  although  on 
the  semi-log  plot  they  do  not  look  large  if  one  adjusts  the  slope 
of  the  straight  line  to  give  the  best  fit. 

In  order  to  see  if  the  independent  oscillator  model  provides 
a  fit  to  the  data  with  realistic  parameters,  we  have  done  the  fol¬ 
lowing,  for  the  case  where  the  molecular  potential  is  assumed  to 
be  the  Morse  potential.  We  have  determined  the  parameters  of  the 
Morse  potential  by  the  use  of  the  value  of  the  TO  frequency,  and 

the  multiphonon  data  on  the  four  alkali  halide  crystals  studied 

(2) 

by  Deutsch  ‘  .  When  this  information  is  combined  with  the  tabu¬ 
lated  value  of  the  interatomic  spacing,  the  coefficient  of 
(linear)  thermal  expansion  may  be  calculated  for  the  model.  We 
shall  describe  the  details  of  the  analysis  below.  The  results 
of  the  investigation  are  summarized  in  Table  1.  For  NaC-t  ,  KCt 
and  KBr,  this  procedure  gives  results  in  remarkable  accord  with 
measured  values  of  the  thermal  expansion  coefficient.  For  LiF, 
the  agreement  is  poorer,  but  the  estimated  and  measured  value  of 
the  thermal  expansion  coefficient  still  differs  by  little  more 
than  a  factor  of  2.  On  the  basis  of  this  analysis,  we  conclude 
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that  multiphoncn  absorption  processes  of  an  intrinsic  character 
are  indeed  responsible  for  the  absorption  measured  by  Deutsch. 

To  obtain  the  numbers  displayed  in  Table  1,  wo  have  employed 
the  following  procedure.  From  Eq. (111-14),  one  sees  that  the 
slope  of  a  plot  of  log  an  vs.  n  is  controlled  only  by  the  param¬ 
eter  D,  for  the  Morse  potentiil.  tfe  have  obtained  the  value  of  D 
given  in  Table  1  by  fitting  the  ratio  c^/q^  to  the  data  ol  Deutsch. 
The  value  of  a  is  obtained  by  identifying  the  frequency  ujq 
^Eq. (1II-3))  with  the  TO  phonon  frequency  of  the  crystal.  The 
parameter  x  in  the  Morse  potential  is  chosen  to  be  equal  to  the 
nearest  neighbor  interatomic  spacing  in  the  crystal. 

In  the  quasi-harmonic  region,  where  kgT<<D,  it  is  an  elementary 
matter  to  calculate  the  me.'n  value  <x>  of  the  interatomic  separa¬ 
tion  of  the  two  atoms  in  the  molecule.  One  finds 

3kBT 

<x>  ~  Xo  +  4aD  * 


where  the  second  term  is  small  compared  to  x  . 
coefficient  of  (linear)  thermal  expansions  with 
where  6(x)  is  the  change  in  the  moan  separation 
duced  by  the  temperature  change  ‘T .  If  we  call 


We  identify  the 
{  ( x ) 

the  ratio  , 

xo  01 

of  the  atoms  pro- 
the  expansion 


coefficient,  then 


Pt 


1u0.\ 

o 


(IV- 2) 


The'  figures  in  the  third  column  of  Table  1  have  been  obtained  from 
Lq  .  (IV- 2)  ,  and  in  the  right  hand  colui.  n  we  give  the  experimental 


See.  E 


the  strength  the  n-phonon  absorption  peak  for  the  Morse  poten- 
tial,  and  for  the  potential  bx  +  a/x  exhibits  nearly  the  same 
functional  dependence  on  n.  This  leads  one  to  expect  that  the 
result  in  Eq.(IV-l)  may  be  rather  general  in  the  limit  of  large  n, 
valid  in  the  classical  regime  for  any  potential  which  admits  an 
harmonic  approximation,  as  long  as  the  anharmonic  corrections  to 


the  particle  motion  are  small.  At  this  time,  we  have  not  suc¬ 
ceeded  in  providing  a  general  proof  of  this  result,  however. 

We  conclude  this  section  with  some  comments  about  the  possible 
role  impurities  may  play  in  affecting  the  behavior  of  the  absorp¬ 


tion  coefficient  at  high  frequencies.  If  the  impurities  are 
coupled  to  the  ions  of  the  host  lattice  by  means  of  an  interac¬ 
tion  that  may  be  crudely  represented  by  either  the  Morse  potential 

2  2 

or  the  potential  bx  +  a/x  ,  then  within  the  framework  of  our 
model,  the  impurities  will  not  affect  the  qualitative  behavior 
of  the  absorption  coefficient,  although  they  will  affect  it  in  a 
quantitative  sense,  since  a  certain  fraction  of  the  lattice  oscil¬ 
lators  will  then  bo  described  by  parameters  which  differ  from  the 
oscillators  which  describe  the  host  lattice. 

However,  there  are  certain  impurity  ions  which  behave  in  an 
anomalous  manner  when  present  as  substitutional  impurities  in 
alkali  halides.  An  example  is  the  Li  +  ion,  which  frequently 
givc.s  rise  to  a  very  low  frequency  resonance  phonon  mode,  even 
though  it  is  a  very  light  ion  which  produces  a  high  frequency 
local  phonon  mode  if  it  is  coupled  to  the  host  ions  by  interac¬ 
tions  characterised  by  harmonic  force  constants  comparable  to 


those  in 


the  host  matrix.  Evidently  in  the  ease  oi  Li,  in  the 
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harmonic  approximation  the  force  constants  are  very  much  smaller 
than  those  which  characterize  the  host,  small  enough  to  offset 
the  tendency  of  the  light  mass  to  create  a  high  frequency  local 
mode.  In  fact,  in  KC*,  the  Li  ion  sits  off  the  substitutional 
latLicc  site,  in  the  (111)  direction  while  it  sits  on  site  in 
KBr .  These  facts  suggest  that  the  Li  ion  moves  in  a  potential 
well  with  a  rather  flat  bottom,  while  terms  higher  order  than 
quadratic  in  the  displacement  of  the  Li+  ion  from  the  substitu¬ 
tional  site  play  an  important  role  in  the  lattice  potential  energy. 
This  notion  is  supported  by  theoretical  studies^,  and  by  experi¬ 
mental  observations  which  show  very  large  electric  field  induced 
shifts  in  the  frequency  of  the  resonance  modc^10l  It  is  also 
true  that  a  number  of  other  ions  have  been  observed  to  produce 
resonance  modes  with  frequencies  very  much  lower  than  the  frequency 
expected  on  the  basis  of  mass  defect  considerations  alone . 


The  remarks  of  the  preceding  paragraph  indicate  that  there 
are  a  certain  number  of  impurities  which  when  placed  in  alkali 
halides  may  be  crudely  described  as  moving  in  the  cage  formed 
by  their  nearest  neighbors,  with  the  floor  of  the  cage  quite  flat 


in  character.  The  calculation  presented  in  section 
contribution  to  the  absorption  coefficient  from  the 
the  square  well  suggests  that  these  impurities  may  g 
bution  to  the  absorption  coefficient  which  falls  oil 
frequencies  large  compared  to  the  r  sonaacc  frequene 
the  peak  in  the  impurity-induced  absorption  occurs, 
presence  of  a  significant  number  of  impurities  which 
very  low  frequency  resonances  may  have  a  significant 
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effect  on  the  behavior  of  the  high  frequency  absorption  coeffi¬ 
cient. 

Before  we  proceed  to  an  estimate  of  the  quantitative  effect 
of  these  impurities  on  the  absorption  coefficient  near  10. 6p,  we 
first  note  that  for  typical  values  of  parameters,  the  character¬ 
istic  frequency  tUp  in  the  discussion  of  section  III(c)  indeed 

_  8 

lies  in  the  proper  spectral  region.  For  example,  if  L =  3  x 10  cm, 
and  if  m  is  chosen  to  be  fifty  atomic  units,  then  ulu,  »  15  cm  \  a 
frequency  in  the  spectral  region  where  the  low  lying  impurity  in¬ 
duced  lattice  resonances  are  observed.  Also,  the  resonance  mode 
frequency  in  some  cases  is  observed  to  increase  significantly 
with  temperatures.  Of  course,  we  cannot  expect  our  very  crude 
model  to  account  for  the  features  of  these  resonance  modes  in  a 
systematic  and  complete  manner,  but  the  overall  qualitative 
features  seem  to  be  reasonable. 

It  is  a  straightforward  matter  to  find  an  expression  for  the 
contribution  of  Nj  impurities  to  the  absorption  coefficient,  in 
the  limit  uj  >>  (Up.  This  may  be  done  by  multiplying  Eq.  (II 1-45) 
by  Nj  ,  the  number  of  impurities  in  the  sample,  dividing  the 

result  by  t lie  (time  averaged)  energy  stored  in  the  electromagnetic 
2 

field  ( C  coV  F,o  /P  rr ,  where  is  the  optical  dielectric  constant 

i 

and  V  the  crystal  volume),  then  multiplying  this  by  Cm2/c,  to 
obtain  1/L,  where  J,  is  the  distance  required  for  the  energy 
density  o J  the  wave  to  decay  to  1/e  of  its  initial  value.  One  finds 
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where  n  r-  NA’  is  the  n.nnber  of  unit  cclls/unit  volume  of  the  host 
crystal  and  f  is  the  impurity  concentration.  In  Kq.(lV-3),  q  and 
m  are  the  effective  charge  and  mass  of  the  impurity.  Let  qQ  and 
in^  be  the  effective  charge  and  reduced  mass  of  the  unit  cell  in 
the  host  crystal,  and  let 

2  in 

R  =  -2-rp  .  (IV-4) 

£  m 

% 

Then  Eq.(IV-3)  may  be  written 


(IV-5) 


2  2 

where  “  Cw  +  4nnqo  /  moa0  is  the  static  dielectric  constant  of 
the  host. 

For  the  purposes  of  providing  a  crude  estimate  of  the  sensi¬ 
tivity  of  the  absorption  coefficient  at  10. 6p  to  the  presence  of 

i 

these  anomalous  impurities,  we  set  R  and  (£s  -  £j)/  €M2  equal  to 
unity,  and  suppose  uu/idq  ~  7  and  u,p  =  20  cm  *  .  We  then  find 


5  f  (cm"1)  (IV- 6) 

where  f  is  the  impurity  concentration. 

The  quantitative  estimate  of  L  displayed  in  Eq.(IV-G)  must 
be  regarded  as  extremely  crude,  because  the  model  is  highly  over¬ 
simplified.  It  does  suggest  that  near  10. Gp,  the  absorption 
coefficient  of  the  crystal  may  be  quite  sensitive  to  small  concen¬ 
trations  of  Li,  Ag,  Cu  or  other  impurities  which  give  rise  to 
impurity  induced  resonance  modes  with  frequency  very  much  lower 
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than  expected  oil  ma* ;;  defect  considerations.  It  would  be 
extremely  i  n  l  ores  t  inr;  to  measure  the  efiect  on  the  absorption 
coeificient  of  doping  KC {  with  Li,  with  concentrations  in  the 
ran^e  of  0.1‘c  to  tost  this  conjecture. 
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ABSTRACT 

The  theory  of  infrared  absorption  by  an  array  of  inde¬ 
pendent,  anharmonic  oscillators  is  discussed.  When 
the  oscillator  potential  is  the  Morse  potential,  the 
theory  provides  an  excellent  description  of  the  tem¬ 
perature  dependence  of  the  absorption  coefficient  at 
10. 6p.  in  NaCt  and  NaF  reported  by  Harrington  and  Hass. 
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There  has  recently  been  interest  in  the  mechanisms  for  absorp¬ 
tion  of  infrared  radiation  by  insulating  crystals  at  frequencies 
high  compared  with  the  Reststrahl  frequency.  In  the  vicinity  of 
the  C02  laser  line  at  10. 6m,  and  in  alkali  halide  crystals  where 
the  electronic  band  gap  is  very  large,  the  principal  ct ntr ibution 
to  the  absorption  coefficient  from  the  bulk  of  the  crystal  comes 
from  multiphonon  processes  in  which  five,  six  or  perhaps  a  larger 
number  of  phonons  are  created  in  the  absorption  process. 

Quite  recently,  Deutsch^  has  completed  a  detailed  experimental 
study  of  the  frequency  dependence  of  the  absorption  coefficient  in 
several  alkali  halide  crystals  at  room  temperature.  Several  groups 

have  addressed  the  theory  of  the  absorption  process  by  models  that 

1-3 

differ  significantly  in  physical  content,  "  and  there  seems  to  be 
general  agreement  that  the  frequency  dependence  of  the  intrinsic 
contribution  to  the  absorption  coefficient  may  be  understood  if 
the  absorption  has  its  origin  in  multiphonon  processes.  All  of  the 
theories  which  have  been  applied  to  the  analysis  of  the  data  pre¬ 
sume  that  anharmonic  effects  on  the  lattice  motion  may  be  treated 
by  perturbation  methods.  One  then  predicts  that  for  temperatures 
T  large  compared  to  the  Debye  temperature  Bp,  the  contribution  to 
the  absorption  coefficient  from  processes  which  involve  n  phonons 

n_  j 

should  vary  with  temperature  like  T 

A  recent  experimental  study  of  the  temperature  dependence  of 
the  absorption  coefficient  or  at  10. 6u  in  several  alkali  halides 

4 

has  been  reported  by  Harrington  and  Hass.  These  authors  find  that 

1 

a  varies  with  T  more  slowly  than  Tn  in  the  crystals  examined  by 
them.  The  purpose  of  this  paper  is  to  apply  the  simple  theoretical 


154 


See.  F 


model  developed  in  our  previous  paper  to  an  analysis  of  this  data. 

We  obtain  an  excellent  quantitative  description  of  the  observed 
temperature  dependence,  and  the  model  also  predicts  absolute  values 
of  a  at  high  temperatures  close  to  the  observed  values.  We  analyze 
the  data  on  NaC£  and  NaF  in  detail.  In  the  case  of  NaCt ,  a  crystal 
considered  in  our  previous  paper,  an  excellent  description  of  the 
data  may  be  obtained  without  introducing  any  parameters  not  found 
in  our  earlier  work  and  for  NaF,  a  crystal  not  examined  earlier  by 
us,  one  parameter  is  required.  The  value  of  this  parameter  may  be 
checked  by  comparing  the  thermal  expansion  coefficient  predicted 
by  our  simple  model  with  the  experimental  result,  and  the  two  agree 
well.  Thus,  we  conclude  that  our  model  provides  an  excellent 
account  of  the  data  reported  by  Harrington  and  Hass.  Furthermore, 
the  physical  content  of  the  model  suggests  that  at  high  temperatures 
anharmonic  effects  have  a  strong  influence  on  multiphonon  excita- 
tations  in  the  alkali  halides. 

In  our  previous  paper,  we  replaced  the  crystal  by  a  set  of  N 
classical,  non-interacting  by  anharmonic  oscillators,  where  N  is 
the  number  of  unit  cells  in  the  crystal.  Each  oscillator  has  a  re¬ 
duced  mass  m,  and  a  transverse  effective  charge  q.  If  the  crystal 
is  illuminated  by  radiation  with  the  elective  field  E(t)  =  Eq  cos  u;t, 

then  the  time  and  ensemble  averaged  rate  at  which  energy  is  absorbed 

3 

by  a  single  oscillator  is  given  by 

2  2 

, .  uu8q  E  ^  i  p00  np  o 

= - _2  2  1  f  dEePEP^(E)6(uj-nw(E))  ,  (1) 

dt  4tn  Z  n=l  b 

where  8  =  l/k„T,  Z  is  the  partition  function  for  a  single  oscillator, 

the  period  T(E)  of  the  bound  orbit  of  energy  E  is  written 
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t  h 

T(E)  -  2n/ou(E),  and  P^CE)  is  a  measure  of  the  amplitude  of  the  n 
harmonic  in  the  orbit  of  energy  E.  If  p^(t)  gives  the  time  depen¬ 
dence  of  the  momentum  in  the  orbit  of  energy  E,  with  the  origin  of 
time  chosen  so  that  the  particle  is  at  a  classical  turning  point 
at  t  -  0,  then 


P  (E) 
n 


2 

tTeT 


.T(E) 


dtp£(t)  sin  nuu(E)t 


(2) 


The  virtue  of  the  model  is  that  exact  expressions  for  a  may  be 

obtained  with  it.,  even  when  the  oscillator  motion  is  very  anharmonic, 

2 

For  the  Morse  potential  V(x)  =  D[  1-exp (a[x~xQ]  )]  we  previously  ob¬ 
tained  exact  expressions  for  uj(E)  and  Pn(E),  although  we  found  a 
only  when  k^T  was  small  enough  for  the  oscillator  motion  to  be 
treated  as  nearly  harmonic.  In  this  paper,  we  work  with  the  full 
form  of  a,  without  this  last  assumption.  In  our  earlier  paper,  we 
found 


u) (E)  -  uuq ( 1-E/D) 


\ 


(3) 


(n+1) 

P  (r,  =  2rr  ( 2mE )  2  /  D\  2  f,  uu(E) 

Pn(E>  „0 -  li) 

where  =  a(2D/m)2  is  the  frequency  of  the  oscillatory  motion,  in 
the  hax'monic  approximation.  Equations  (3)  and  (4)  may  be  inserted 
into  Eq.  (1),  and  the  integration  carried  out.  We  make  one  simpli¬ 
fication  in  the  exact  result.  For  k^T  <<  D,  a  limit  that  applies 
to  our  discussions  here,  the  partition  function  is  well  approximated 
by  the  result  Z  =  2^kgT/oio  obtained  from  the  harmonic  approximation.' 
In  the  results  that  follow,  we  use  this  form  for  Z.  Then  if  we  let 
§  -  ul)/uj0  ,  and  nm  is  the  first  integer  larger  than  l,  we  find 


n 

,  (4) 
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At  fixed  frequency,  the  temperature  dependence  of  the  absorp¬ 
tion  coefficient  is  controlled  by  the  single  parameter  D/k^T,  which 

o 

for  the  case  of  NaCt  may  be  obtained  from  our  earlier  work. 


We  first  recover  from  Eq.  (5)  the  result  valid  in  the  quasi¬ 
harmonic  regime,  where  for  the  nth  term,  this  requires  nkBT  «  D. 

t  h 

In  this  limit,  the  n  term  of  Eq.  (5)  peaks  sharply  just  below 
5  =  n.  Thus,  let  §  =  n  -  €,  with  €  small.  If  the  nth  term  in 
Eq.  (5)  is  denoted  by  (<“jf)>n, 


-■  2TTq  E°  /  D  \2  €n9(€) 

‘dt;/n  muu  \k„T/  „n  n+1 

o  '  B  '  Z  n 


2D 

nkBT 


(6) 


where  0(€'  =  1  if  €  >  0  and  is  zero  for  €  <  0.  The  integrated 
strength  of  the  absorption  line  in  Eq.  (6)  is 


This  expression  is  identical  to  that  produced  by  the  quasi-harmonic 
approximation  utilized  earlier  by  us.^ 

In  Figure  1,  we  present  a  calculation  of  the  frequency  depen¬ 
dence  of  the  absorption  coefficient  for  the  model  for  T  =  300° K  and 
T  *=  900°K,  fron  Eq.  (5).  We  have  used  the  value  of  D  obtained  for 


157 


ABSORPTION  COEFF.  (CM"1) 


Sec.  F 

NaCt,  from  our  earlier  work,  and  with  the  absolute  value  of  the  a 
adjusted  to  fit  the  10. 6U  data  of  Harrington  and  Hass  at  900°  K. 

Several  features  of  these  results  deserve  comment.  For  T=300  K, 
the  n=3  term  in  Eq.(5)  produces  a  very  sharp  peak  centered  a  bit 
below  ^  =  3.  This  peak  is  quite  narrow  and  well  defined,  as  one  would 

I 

expect  in  the  quasi-harmonic  approximation.  For  £  near  5  or  6,  the 
effect  of  anharmonicity  is  severe  enough  that  a  varies  smoothly  with 
frequency,  displaying  only  gentle  shoulders  as  a  reminder  of  the 
sharp  structure  present  in  the  quasi-harmonic  approximation.  But  the 
time  T=900°K,  the  theory  produces  a  very  smooth  dependence  of  a 
with  frequency. 

These  calculations  suggest  that  for  large  values  of  u>/<V  even 
at  room  temperature,  the  lattice  motion  cannot  be  regarded  as  nearly 
harmonic,  so  the  absorption  coefficient  cannot  be  calculated  by  per¬ 
turbation  theoretic  methods  which  tre^t  the  anharmonic  terms  in  the 
crystal  Hamiltonian  as  small.  This  is  also  clear  from  the  data  of 
Harrington  and  Hass,  which  we  shall  see  is  well  fitted  by  our  model, 
since  the  data  show  very  large  departures  from  the  Tn  1  behavior 

cited  earlier. 

It  must  be  emphasized  that  at  fixed  T,  the  relative  importance 
of  the  anharmonicity  increases  as  u)Adq  increases.  If  we  examine  the 
absorption  coefficient  for  the  model  near  the  Reststrahl  region  (the 
Reststrahl  absorption  is  described  by  the  term  n  =  1  in  Eq.  (5)), 
then  the  parameters  we  employ  produce  rather  modest  anharmonic 
effects.  For  example,  at  temperature  T,  the  Reststrahl  peak  occurs 
very  close  to  uJq  -  uJo(l-kBT/2D)  ,  where  u)Q  is  the  Reststrahl  frequency 
at  T  -  0.  For  our  parameters,  at  room  temperature,  kgT/2D  ==0.03 
for  NaCt.  The  half  width  at  half  maximum  is  also  close  to  kgT/2D 
in  magnitude,  and  increases  linearly  with  T.  As  the  parameter  w/% 


160 


Sec.  F 


increases,  each  successive  peak  shifts  to  lower  frequency  by  a  frac¬ 
tional  amount  that  increases  with  increasing  order,  and  the  width  of 
each  feature  described  by  the  terms  in  Eq.  (5)  increases  rapidly 
enough  so  that  by  the  time  f,  =  6 ,  the  absorption  coefficient  varies 
smoothly  with  frequency  at  room  temperature. 

In  Figure  2a,  we  compare  the  temperature  dependence  of  the  ab¬ 
sorption  coefficient  at  10. 6^  obtained  from  our  calculations  with 
the  data  of  Harrington  and  Hass.  The  theory  has  been  fitted  to  the 
data  at  900rK,  and  gives  a  very  good  account  of  the  observed  temper¬ 
ature  dependence  for  all  but  the  lowest  temperature.  In  his  previous 
work,  Deutsch  has  reported  a  value  of  0.001  cm’1  for  the  absorption 
coefficient  at  room  temperature,  while  the  value  reported  by  Harring¬ 
ton  and  Hass  is  larger  by  a  factor  of  2.7.  Thus,  near  room  tempera- 
ture,  the  absorption  coefficient  measured  by  Harrington  and  Hass 
presumably  contains  an  extrinsic  contribution  which  is  dominated  by 
the  intrinsic  contribution  at  higher  temperatures.  If  we  accept 
Deutsch  s  value  as  the  correct  one  at  room  temperature  (see  the 
diamond  in  Fig.  2a),  then  we  may  fit  a  to  within  a  factor  of  two 
over  the  entire  temperature  range. 

Our  model  also  predicts  the  absolute  magnitude  of  the  absorp¬ 
tion  coefficient.  The  rate  at  which  energy  is  absorbed  from  the 
field  is  found  by  multiplying  Eq.  (5)  by  the  number  of  molecules  N 
(equal  to  the  number  of  unit  cells)  in  the  crystal.  Then  upon 
dividing  by  the  time  averaged  energy  V€qE^/8tt  stored  in  the  field, 
where  €q  is  the  high  frequency  dielectric  constant,  and  dividing  by 
the  propagation  velocity  c/€^,  one  obtains  the  inverse  of  the  length 
d  required  for  the  energy  density  in  the  beam  to  decay  to  1/e  of  its 
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initial  value.  This  quantity  is  the  absorption  coefficient  a  meas¬ 
ured  in  the  experiments.  The  magnitude  of  a  may  be  estimated  in 

the  multiphonon  regime  from  dielectric  constant  data  once  D  is  known 

2 

since  €  -  €  =  4nnq  /m  for  ionic  crystals.  For*NaCt,  the  theory 

3  U 

predicts  that  at  T  =  900°K,  the  absorption  length  should  be  0.011  Cm" 
at  10. 6u,  while  the  data  shows  it  to  be  0.020  cm-1.  Our  model  thus 
gives  an  absolute  value  for  a  in  good  accord  with  the  data. 


In  Figure  2b,  we  compared  the  temperature  dependence  predicted 

for  the  absorption  coefficient  at  10. 6U  with  the  data  on  NaF.  Again 

we  obtain  an  excellent  fit  for  500° K  <  T  <  1200° K.  The  theory 

does  poorly  for  T  <  500° K,  presumably  because  8^  «  490°  in  NaF, 

and  our  classical  model  is  inapplicable  for  T  s  0^.  In  Figure  2b, 

we  have  chosen  f>NaF  =  DNaC^  >  and  adjusted  the  magnitude  of  the 

absorption  length  to  fit  the  data  at  700° K.  For  NaF,  §  s  3.96  at 

10. 6^.  The  value  selected  for  D  predicts  the  coefficient  of  linear 

—  6  ®  —  1 

thermal  expansion  for  the  crystal  to  be  38  x  10  K  when  T  >  0^ , 

if  we  use  our  earlier  procedure  to  make  this  estimated  This  value 

is  in  excellent  accord  with  the  measured  value  of  the  thermal 

7 

expansion  coefficient.  From  the  dielectric  constant  data,  the 
theory  predicts  the  absorption  length  to  be  1.2  cm  ^  at  900° K, 
while  the  measured  value  is  3.4  cm-1.  Thus,  while  the  absolute 
value  for  ot  in  NaF  agrees  less  well  with  the  data  than  in  the  case 
of  NaCt,  the  theory  gives  reasonable  semiquantitative  agreement. 


One  can  inquire  about  sensitivity  of  these  results  to  the 

v 

details  of  the  potential.  To  test  this,  the  dashed  curve  in 

Fig.  (2a)  gives  the  dependence  on  T  for  a  at  6uiq  for  the  potential 

2  2  (3) 

bx  +  a/x  examined  previously  ,  with  the  fit  to  the  room 
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temperature  data  carried  out  as  in  our  earlier  paper.  This  potential 
also  provides  a  reasonable  fit  to  the  data,  so  the  detailed  form 
of  the  potential  does  not  seem  to  be  critical. 

Thus  our  model  gives  a  good  account  of  the  dependence  of  or 
on  T  in  alkali  halides  through  the  use  of  only  a  single  parameter 
whose  value  may  be  checked  through  use  of  data  on  the  thermal  ex¬ 
pansion  coefficient.  Quite  recently,  Sparks  and  Sham8  have  extended 
their  earlier  work2  to  provide  an  account  of  the  dependence  of  & 
on  T.  Their  theory  includes  the  effect  of  phonon  dispersion  in  an 
approximate  way.  They  also  treat  the  problem  by  quantum  mechanical 
methods  while  our  theory  is  purely  classical.  Thus,  their  theory 
may  be  applied  to  the  region  T  <  eD  while  ours  may  not  be.  The  two 
pieces  of  work  are  complementary  in  a  certain  sense.  While  their 
model  is  more  detailed  than  ours,  they  introduce  a  sequence  of 
approximations,  such  as  the  inclusion  of  +erms  only  through  a  finite 
order  of  perturbation  theory,  approximate  treatment  of  phonon  dis¬ 
persion,  and  finally  a  phenomenological  inclusion  of  certain  anharmonic 
effects  by  allowing  the  phonon  frequencies  to  be  temperature  dependent. 
Our  model  is  more  schematic  in  nature,  but  its  virtue  is  the  use  of 
only  a  small  number  of  parameters  to  characterize  a  given  material, 
and  we  obtain  exact  results  for  or.  The  two  theories  concur  in  one 
important  regard.  At  high  temperature,  and  in  the  multiphonon  regime, 
in  the  alkali  halides  the  effect  of  anharmonic ity  is  very  large,  so 
large  deviations  from  the  Tn_1  law  occur. 
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G.  TEMPERATURE  AND  FREQUENCY  DEPENDENCE  OF 
INFRARED  ABSORPTION  AS  A  DIAGNOSTIC  TOOL* 

M.  Sparks 

Xonics,  Incorporated,  Van  Nuys,  California  91406 

Recent  developments  render  untenable  a  proposed  method 
of  distinguisi  ing  between  intrinsic  and  extrinsic  infrared  ab¬ 
sorption  on  the  basis  of  the  proposed  temperature  dependence. 

However,  when  the  proper  temperature  dependence  of  multi - 
phonon  absorption  is  accounted  for  and  the  possibility  of  other 
intrinsic  processes  is  taken  into  account,  the  temperature  and 
frequency  dependence  of  the  absorption  of  both  the  best  avail¬ 
able  and  intentionally  imperfected  crystals  should  be  useful 
in  studying  extrinsic  processes. 

It  has  been  suggested  that  the  temperature  dependence  of  the  optical  absorption 
coefficient  )3  could  be  used  to  distinguish  between  intrinsic  (characteristic  of  a 
perfect  crystal)  and  extrinsic  (caused  by  imperfections)  infrared  absorption. 
Several  recent  developments  bear  on  this  suggestion.  First,  a  c<  mbination  of 
frequency  o o  and  temperature  T  dependence  should  be  far  more  useful  than  the 
T  dependence  alone. ^  Second,  the  T  dependence  of  $  for  mult,  phonon  absoiption 
deviates  considerably^’^  from  the  expected^'  (j’  *  result.  Third,  other  extrinsic 

processes  in  addition  to  those  considered  in  Ref.  1  should  be  included. 

Consider  the  first  development.  The  near  exponential  frequency  dependence 
of  0  observed'*8  in  a  number  of  materials  including  LiF,  NaF,8  NaCl,  KC1,  KBr, 
MgF2,  CaF2,  BaF2,  SrF2 ,  MgO,  A12C>3,  SiC>2  ,  TiC>2  ,  SrTiOg ,  and  BaTiOg , 
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when  extrapolated  to  10.  6 Jim,  gives  a  fair  estimate  of  the  intrinsic  value  of  fl 
at  10.  6jim.  Thus,  it  is  already  known  with  fair  accuracy  whether  the  10.  6jim 
values  of  j3  are  intrinsic  or  extrinsic  for  a  number  of  materials.  Nevertheless, 
additional  to  and  1  measurements  would  be  of  great  interest  in  studying  the 
sources  of  extrinsic  absorption  and  possibly  new  intrinsic  absorption  mechan¬ 
isms.  Intentionally  introducing  various  types  of  imperfections  into  crystals  and 

9 

on  their  surfaces  should  be  invaluable  in  such  studies.  Since  values  of  /3  of 
-4  -1 

10  cm  and  even  lower  are  of  current  interest,  standard  transmission  measure¬ 
ments  are  inadequate.  Emissivity  techniques  should  be  capable  of  yielding  values 

-4-1  9 

of  well  below  10  cm  in  a  carefully  designed  simple  instrument.  Calori¬ 
metric  measurements  of  j3’s  approaching  10_4cm'1  at  the  single  wavelength 
10.  6jim  are  of  course  now  common. 

Concerning  the  second  development,  the  temperature  dependence  of  the 

intrinsic  n-phonon  absorption  (annihilation  of  one  photon  and  creation  of  n  phonons) 

S  6  2  1 

is  considerably  weaker  than  the  well  known'  ’  '  ' 1  explicit  temperature  dependence 

„a)30  (l-e'“/WT)  (  1  -  e'“/n“T  (1) 

where  to  is  the  laser  frequency  and  tOy  =  k^T /1i .  Eq.  (1),  which  can  be  written 
as  p  ~  (nq  + 1 )  -  nQ  ,  where  ng  is  the  Bose-Einstein  occupation  number  for 
a  phonon  of  frequency  to/  n ,  or  in  several  other  simple  forms  by  using  energy 
conservation,  reduces  to  /3  ~  Tn  1  in  the  high-temperature  limit  to^  »  to. 
Experimentally,  this  deviation  from(l)  was  observed  by  Harrington  and  Hass,3 
and  is  apparent  in  the  data  of  Barker.10  Sparks  and  Sham4  have  explained  this 
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deviation  by  including  the  temperature  dependence  of  tile  phonon  frequencies 
and  lattice  constant  in  their  previous  theory  of  the  nearlv  exponential  fre¬ 
quency  dependence  of  j9  discussed  below.  In  addition,  the  particular  extrinsic 
processes  suggested  bv  Hardy  and  Agrawal  ^  and  bv  Rosenstock  ^  will  show 
similar  deviations  from  their  predicted  explicit  T  dependence  when  the  T 
dependence  of  the  parameters  is  included,  lentil  more  experimental  informa¬ 
tion  is  available  to  establish  faith  in  the  ability  to  predict  the  T  dependence  of 
the  multiphonon  and  band-phonon  plus  localized-impurit\  processes,  distin¬ 
guishing  between  the  two  on  the  basis  of  the  1'  dependence  alone  probablv 
would  be  difficult. 

The  third  development  concerns  a  further  more  serious  problem  in  dis¬ 
tinguishing  between  intrinsic  and  extrinsic  absorption  on  the  basis  of  the  T  de¬ 
pendence  in  Ref.  1.  That  is,  there  are  many  possible  extrinsic  processes  in 
addition  to  the  band-phonon  plus  localized-impurity  process  suggested  in  Refs.  1 
and  6.  For  example,  small  amounts  of  macroscopic  inclusions  can  give  rue 
to  a  temperature  dependence  ranging  from  £  ~  T  to  T  in  typical  cases, 
depending  on  the  type  of  impurity.^  There  is  also  some  experimental  evi¬ 
dence,  though  it  is  not  conclusive  at  present,  for  temperature  independent 
12  1 

absorption.  ’ 

The  study  of  intrinsic  and  especially  extrinsic  processes,  by  measuring 
0(oc,T)  in  ultrapure  and  intentionally  imperfected  crystals,  for  example, 
indubitably  will  become  important  if  the  current  interest  in  obtaining  low- 
absorption  materials  is  unabated.  Already  there  are  a  number  of  cases  in 
which  multiphonon  absorption  is  obviously  the  source  of  0  and  others  in  which 
it  obviously  is  not.  For  example,  there  is  now  little  doubt  that  the  nearlv 
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exponential  frequency  dependence  ’  of  0  mentioned  above  is  the  result  of 

multiphonon  absorption.  The  contributions  of  the  individual  n-phoncn  processes 

have  been  calculated  and  summed  to  obtain  the  nearly  exponential  frequency 

dependence.  '’  More  recently,  a  direct  derivation  of  the  exponential  in  closed 

form,  rather  than  as  a  sum  of  the  individual  n-phonon  contributions,  has  been 
13 

given. 

The  fact'  that  the  10.6  pm  values  of  0  for  a  number  of  crystals,  such  as 
KBr,  CdTe,  and  KC1  as  examples,  lie  well  above  the  values  obtained  by  extra¬ 
polating  longer -wavelength  values  to  10.6  (im  surely  indicates  that  these  meas¬ 
ured  values  of  0  arise  from  processes  other  than  multiphonon  absorption.  A  clear 

example  of  an  extrinsic  experimental  value  of  0  is  that  of  a  NaCl  sample  meas- 

3 

ured  at  10.6pm  from  300  K  to  near  the  melting  point  by  Harrington  and  Hass. 

The  room  temperature  value  of  0  =  0.003  cm  1  for  this  sample  is  considerably 

- 1  14 

greater  than  both  the  lowest  value  of  0  =  0.0015  cm  observed  to  date  and  the 

-3  - 1 

estimated  intrinsic  value  of  slightly  less  than  10  cm  obtained  by  extrapolating  the 
lower  frequency  data  to  10.6pm.  Furthermore,  the  temperature  dependence  deviates 
strongly  from  the  multiphonon  value  evjn  when  the  T  dependence  of  phonon  fre¬ 
quencies  and  lattice  constant  are  taken  into  account.  In  fact,  0  decreases  as 
the  temperature  increases  near  room  temperature. 
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H.  SHORT-PULSE  OPERATION  OF  INFRARED  WINDOWS 
WITHOUT  THERMAL  DE FOCUSING* 

M.  Sparks 

Xonics,  Incorporated,  Van  Nuys,  California  91406 

The  possibility  of  transmitting  short  infrared  pulses 

through  materials  with  little  thermally  induced  optical 

distortion  is  shown  to  exist.  For  sufficiently  short  pulses, 

-8  -9 

of  the  order  of  10  —  10  sec,  the  absorbed  energy'  does  not 

have  time  to  thermalize,  thus  avoiding  heating  effects  until 
after  the  pulse  has  been  transmitted. 

In  high-power  laser  systems,  heating  of  a  window  or  other  transparent 
component  by  the  laser  beam  causes  changes  in  the  index  of  refraction  n  and 
in  the  thickness  of  the  window.  ’  ’  The  resulting  defocusing  is  one  of  the 
most  serious  problems  in  such  systems.  It  has  been  pointed  out  that  much 
higher  intensities  can  be  transmitteo  if  the  laser  can  be  operated  for  a  period 
of  time,  say  of  the  order  of  a  second,  after  which  the  window  is  cooled  before 
the  next  pulse  is  transmitted.  ’  *  Hloembergen  suggested  that  if  the  pulse 
duration  is  much  shorter  than  a  characteristic  time,  of  the  order  of  an  acousti¬ 
cal  velocity  times  a  linear  dimension  of  the  window,  the  change  in  window 
thickness  is  negligible,  and  only  changes  in  n  (at  constant  strain)  should  be 
considered. 

The  purpose  of  the  present  note  is  to  point  out  that  in  principle  still 
shorter  pulse  durations  allow  the  window -to-transmit-thtrpplsc  before  the 
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temperature  rises,  thus  avoiding  the  thermal  defocusing.  The  requirement 
for  this  effect  is 


'puls  *  "r<w> 


where  tpUjg  is  the  pulse  duration  and  r(to)  is  the  relaxation  frequency  of  the 
fundamental -phonon  mode  (transverse  optical  mode  with  wavevector  0). 
When  (1)  is  satisfied,  the  energy  absorbed  from  the  laser  field  by  the  funda¬ 
mental  mode  has  not  had  sufficient  time  to  relax  out  of  the  fundamental  mode 
and  thermalize.  Thus  the  temperature,  which  corresponds  to  the  thermal- 
equilibrium  values  of  the  phonon  occupation  numbers,  remains  at  its  initial 
value. 

Practical  values  of  t  jg  depend  critically  on  the  deviation  of  r(oo)  from 
the  resonant  value  =  T  (ocf),  where  tof  is  the  frequency  of  the  fundamental 
mode.  For  example,  from  the  value  of  the  linewidth  in  thin  NaCl  samples,  ^ 


Tf  2=  2.5  x  1012  sec"1 


r>  -13 

The  corresponding  value  of  t  <  4  x  10  sec  from  (1)  is  smaller  than 
values  of  current  interest.  The  value  of  T  (to)  decreases  rapidly  as  to 
increases  above  tOj- .  On  resonance  ,  is  large  in  general  because  there 
are  many  phonon  states  into  which  the  fundamental  mode  can  decay.  As  to 
increases,  the  lower-order  processes,  in  which  the  fundamental  phonon  splits 
into  only  a  few  phonons,  cannot  conserve  energy,  and  T  (to)  decreases  ac¬ 
cordingly. 
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In  order  to  obtain  an  estimate  of  the  size  of  T  (to),  consider  the  optical 

6,7 

absorption  coefficient  /3(to) 


/3(co)  =  A 


_ to  r<to) _ 

/  2  2 .2  2  .2 

(to  -  tOj  )  +  tOj  1  (to) 


(2) 


where  A  is  *  constant.  On  resonance,  (2)  gives 

0(cof)  =  A/cof  rf  .  (3) 

At  the  laser  frequency  to,  assumed  to  be  high  [to  »  tof,  ]?(to)],  (2)  gives 
0  (to)  =  AT(to)/to3  .  (4) 

Tiie  ratio  of  /3(tof)  to  |3 (to)  from  (3)  and  (4)  is 

0(tof)/0(to)  =  to3/tof  rfT(to) 

which  gives,  with  (1),  the  central  result 

^  j  $(tOj )  tOj 

lpuls  <  T(to)  "  J(u)  ^3 


(5) 


(6) 


Consider  the  materials  in  which  distortion -free  transmission  is  likely  to  be  ob¬ 
servable.  In  order  to  make  t  lg  large,  T(to)  should  be  small,  according  to  (6). 

3 

Since  0 (to/ tof)  decreases  more  rapidly  than  (co/tof)  with  increasing  to  in 

7 

the  multiphonon  absorption  region,  a  large  value  of  to/tof  is  desirable  in 
order  to  make  the  denominator  in  (6)  small.  However,  if  to/tof  is  too  large, 
other  absorption  processes  dominate  the  multiphonon  absorption.  Thus,  the 
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smallest  value  of  j3(oj)  in  the  multiphonon  region  is  desired.  The  smallest 

_4  -1  8 

values  of  p(oo)  measured  to  date  at  10.  6  Jim  are  of  the  order  of  10  cm  . 


The  value  of  0  at  10.  6 Jim  for  KC1  obtained  by  extrapolating  from  mea 

8  o  i  4  1 

sured  values  between  0  =20  and  5x  10  cm  is  0  =  10  "cm"  .  In  view  of 

-4  - 1 

current  interest  in  KC1  and  the  fact  that  0  =  4  x  10  cm  has  already  been 

-4  -1 

demonstrated,  it  is  likely  that  10  cm  will  be  attained.  Thus,  KC1  is  a 
good  candidate  in  which  to  study  the  effect.  At  other  wavelengths,  oth^r  ma¬ 
terials  would  be  more  appropriate;  for  example,  the  extrapolated  value  of  j3 
equals  10  4cm  *  at  ~3.  5iim  for  LiF.  If  values  of  /3  smaller  than  10  4cm  ^ 
are  attained,  the  choice  of  materials  will  change. 

For  KC1,  with8  to 2.7  x  10*3  sec  *,  2:  0.  08  oof,  uo  =  1.  8  X  10*"* 

sec  1  (i.e. ,  10.6ym),  /3 (cof )  s-  5  x  I04cm‘*,  and  0(u  )  ^  10_4cm’*, 

(6)  gives 


1 

T(u» 


5  x  10 


sec 


(7) 


Such  short  pulses  are  not  presently  available  at  lOfirn.  However,  they 
are  currently  available  at  shorter  wavelengths.  For  example,  6TW  with  a 
pulse  duration  of  10  **  sec  is  available  at  1.06Mm.  Picosecond  pulse  dura¬ 
tions  are  theoretically  possible  at  longer  wavelengths,  of  course.  Values  of 
l/r(to)  at  wavelengths  shorter  than  10jim  can  be  as  large  as  that  in  (7).  For 
example,  at  3.5jim  for  LiF  with  u:f  =  5.  8  x  10*3  sec’*,  =  0.08  , 

co  =  S.5  x  10*4  sec  *,  (u)j)  =  2  x  lo'’  cm  *,  and  /3(ce)  =  10  4  cm  *,  (6)  gives 

_9  - — - 

l/r(to)  =  3  x  10  sec.  Thus,  it  is  possible  that  the  no-heating  short-pulse 
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effect  could  be  observed  over  a  substantial  part  of  the  near  infrared  region. 
Unfortunately,  at  1.06fim,  where  10  ^  sec-pulse-duration  lasers  are  avail¬ 
able,  there  is  no  material  with  /3  determined  by  multiphonon  absorption.  I'lius, 
tests  and  applications  of  the  theory  will  have  to  await  the  availability  of  short - 
pulse  sources  at  longer  wavelengths. 

Hie  analysis  above  applies  to  the  case  of  multiphonon  absorption  bv  the 
anharmonic  potential  mechanism.  The  possibility  of  ?.  similar  no-heating 
effect  at  other  frequencies  and  for  other  absorption  mechanisms  also  exists. 
However,  the  situation  is  likely  to  be  considerably  more  complicated.  For 
example,  for  absorption  in  the  visible  and  very  near  infrared,  the  absorption 
is  likely  to  involve  electronic  excitations,  such  as  in  electric-field-induced 
absorption  or  excitation  of  impurity  -type  or  surface-state  levels.  The  time 
required  for  the  energy'  to  be  transferred  to  the  lattice  varies  greatly  with  the 
number  of  impurities,  type  of  material,  and  the  temperature  in  general. 
Furthermore,  the  direct  effect  on  the  index  of  refraction  of  the  electronic 
processes  involved  in  the  absorption  and  the  possibility  of  induced  transparency 
from  emptying  the  impurity  levels  should  be  considered. 

In  addition  to  being  of  fundamental  interest,  the  short-pulse  effect  could 
have  practical  applications.  It  could  be  useful  for  diagnosis,  as  in  distinguish¬ 
ing  between  multiphonon  absorption  and  electronic -type  absorption,  for  example. 
It  could  possibly  afford  a  method  of  measuring  the  relaxation  frequency  of 
the  fundamental  phonon  mode  directly.  Its  usefulness  in  high-power  applications 
is  limited  by  material  breakdown  at  the  high  intensities  required  to  obtain  high 
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power  with  short  pulse  times.  A  typical  value  of  breakdown  intensity  for  alkali 
10  2  M  -9 

halides  is  2  x  10  W/cm  .  For  t  jg  =  5  x  10  sec,  this  gives  100  joules 

2 

per  pulse,  or  100  W/cm  average  for  a  duty  cycle  of  1  sec.  This  value,  which 

possibly  could  be  raised  by  operating  at  lower  than  room  temperature  or  further 

4  2 

improving  materials,  can  be  compared  with  values  of  ~0. 2  to  ~2  x  10  W/cm 
estimated  to  be  obtainable  for  candidate  window  materials  (ranging  from  Si  to 
KBr)  for  a  one  second  pulse,  two  second  duty  cycle,  beam  truncated  at  ^  its 
maximum  intensity,  and  halving  the  intensity  at  tl  e  target.*  If  greater  than  half 
the  intensity  must  be  available  at  the  target  or  if  the  beam  is  truncated  at  a  lower 
intensity  (as  it  usually  is),  these  one-second-pulse  values  are  reduced,  by  fac¬ 
tors  of  a  thousand  or  even  greater  in  extreme  cases.  In  such  cases,  the  no¬ 
heating  effect  would  be  useful. 

It  should  be  mentioned  that  local  heating  of  imperfections,  which  can  lead  to 
material  failure,  is  a  greater  problem  for  short  pulses  than  long  ones  for  a  given 

9 

value  of  the  energy  per  pulse.  The  temperature  rise  tends  to  be  lower  for  long 
pulses  since  the  heat  has  time  to  diffuse  away  from  the  generation  site.  In  this 
regard,  notice  that  when  (1)  is  satisfied,  the  time  constant  for  the  energy  to 
leave  the  fundamental  mode  (after  the  pulse  is  transmitted)  is  1/1^,  which  is 
considerably  shorter  than  l/r(u.')»  since  the  fundamental  mode  oscillates  at  its 
resonant  frequency  after  the  drive  field  is  turned  off.  A  further  limit  to  the 
high-power  use  of  the  short  pulse  effect  is  the  air  breakdown  at  the  high  in¬ 
tensities  encountered  in  high-energy  short-pulse  systems. 

Discussions  with  Dr.  C.  J.  Duthler  are  gratefully  acknowledged. 
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